Skip to content

Commit

Permalink
fix nuopc cap for radiation stress components
Browse files Browse the repository at this point in the history
  • Loading branch information
uturuncoglu committed Jul 15, 2024
1 parent d9b3172 commit e8b3ef6
Showing 1 changed file with 9 additions and 6 deletions.
15 changes: 9 additions & 6 deletions model/src/wav_import_export.F90
Original file line number Diff line number Diff line change
Expand Up @@ -145,6 +145,9 @@ subroutine advertise_fields(importState, ExportState, flds_scalar_name, rc)
call fldlist_add(fldsFrWav_num, fldsFrWav, 'Sw_vstokes')
else
call fldlist_add(fldsFrWav_num, fldsFrWav, 'Sw_z0')
call fldlist_add(fldsFrWav_num, fldsFrWav, 'Sw_wavsuu')
call fldlist_add(fldsFrWav_num, fldsFrWav, 'Sw_wavsuv')
call fldlist_add(fldsFrWav_num, fldsFrWav, 'Sw_wavsvv')
end if
call fldlist_add(fldsFrWav_num, fldsFrWav, 'Sw_pstokes_x', ungridded_lbound=1, ungridded_ubound=3)
call fldlist_add(fldsFrWav_num, fldsFrWav, 'Sw_pstokes_y', ungridded_lbound=1, ungridded_ubound=3)
Expand Down Expand Up @@ -753,14 +756,14 @@ subroutine export_fields (gcomp, rc)
call CalcBotcur( va, wbcuru, wbcurv, wbcurp)
end if

if ( state_fldchk(exportState, 'wavsuu') .and. &
state_fldchk(exportState, 'wavsuv') .and. &
state_fldchk(exportState, 'wavsvv')) then
call state_getfldptr(exportState, 'sxxn', sxxn, rc=rc)
if ( state_fldchk(exportState, 'Sw_wavsuu') .and. &
state_fldchk(exportState, 'Sw_wavsuv') .and. &
state_fldchk(exportState, 'Sw_wavsvv')) then
call state_getfldptr(exportState, 'Sw_wavsuu', sxxn, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
call state_getfldptr(exportState, 'sxyn', sxyn, rc=rc)
call state_getfldptr(exportState, 'Sw_wavsuv', sxyn, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
call state_getfldptr(exportState, 'syyn', syyn, rc=rc)
call state_getfldptr(exportState, 'Sw_wavsvv', syyn, rc=rc)
if (ChkErr(rc,__LINE__,u_FILE_u)) return
call CalcRadstr2D( va, sxxn, sxyn, syyn)
end if
Expand Down

0 comments on commit e8b3ef6

Please sign in to comment.