Changeset 11908 for NEMO/branches
- Timestamp:
- 2019-11-14T20:20:07+01:00 (5 years ago)
- Location:
- NEMO/branches/2019/ENHANCE-02_ISF_nemo
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo/cfgs/SHARED/field_def_nemo-oce.xml
r11835 r11908 264 264 265 265 <!-- * variable related to ice shelf forcing * --> 266 <field id="isfdebug3d2" grid_ref="grid_T_3D" /> 267 <field id="isfdebug3d1" grid_ref="grid_T_3D" /> 268 <field id="isfdebug2d" /> 266 <field id="isftfrz_par" long_name="fzp temperature at ocean/isf interface" unit="degC" /> 267 <field id="isftfrz_cav" long_name="fzp temperature at ocean/isf interface" unit="degC" /> 269 268 <field id="fwfisf_cav" long_name="Ice shelf melting" unit="kg/m2/s" /> 270 269 <field id="fwfisf_par" long_name="Ice shelf melting" unit="kg/m2/s" /> -
NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcpl.F90
r11902 r11908 61 61 ! check presence of variable needed for coupling 62 62 ! iom_varid return 0 if not found 63 id = 063 id = 1 64 64 id = id * iom_varid(numror, 'ssmask', ldstop = .false.) 65 65 id = id * iom_varid(numror, 'tmask' , ldstop = .false.) … … 145 145 !!---------------------------------------------------------------------- 146 146 !! 147 INTEGER :: ji, jj, jd 147 INTEGER :: ji, jj, jd, jk !! loop index 148 148 INTEGER :: jip1, jim1, jjp1, jjm1 149 149 !! … … 195 195 ! 196 196 ! recompute the vertical scale factor, depth and water thickness 197 IF(lwp) write(numout,*) 'isfcpl_ssh : recompute scale factor from sshn (new wet cell)' 198 IF(lwp) write(numout,*) '~~~~~~~~~~~' 199 DO jk = 1, jpk 200 e3t_n(:,:,jk) = e3t_0(:,:,jk) * ( ht_0(:,:) + sshn(:,:) ) & 201 & / ( ht_0(:,:) + 1._wp - ssmask(:,:) ) * tmask(:,:,jk) & 202 & + e3t_0(:,:,jk) * (1._wp -tmask(:,:,jk)) 203 END DO 204 e3t_b(:,:,:) = e3t_n(:,:,:) 197 205 CALL dom_vvl_zgr() 198 206 ! -
NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/TRA/traisf.F90
r11902 r11908 68 68 ! ensure 0 trend due to unconservation of the ice shelf coupling 69 69 IF ( ln_isfcpl_cons ) CALL tra_isf_cpl(risfcpl_cons_tsc, tsa) 70 71 IF ( ln_isfdebug ) THEN 72 CALL debug('tra_isf: risfcpl_tsc T',risfcpl_tsc(:,:,1)) 73 CALL debug('tra_isf: risfcpl_tsc S',risfcpl_tsc(:,:,2)) 74 END IF 70 75 ! 71 76 END IF 72 77 ! 73 78 IF ( ln_isfdebug ) THEN 74 CALL debug('tra_isf: risfcpl_tsc T',risfcpl_tsc(:,:,1))75 CALL debug('tra_isf: risfcpl_tsc S',risfcpl_tsc(:,:,2))76 79 CALL debug('tra_isf: tsa T' ,tsa(:,:,:,1)) 77 80 CALL debug('tra_isf: tsa S' ,tsa(:,:,:,2)) -
NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/isfstp.F90
r11889 r11908 145 145 !--------------------------------------------------------------------------------------------------------------------- 146 146 ! initialisation ice sheet coupling 147 IF ( ln_isfcpl ) CALL isfcpl_init()147 IF ( ln_isfcpl ) CALL isfcpl_init() 148 148 ! 149 149 END IF -
NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/sbcfwb.F90
r11896 r11908 84 84 IF( kn_fwb == 2 ) WRITE(numout,*) ' adjusted from previous year budget' 85 85 IF( kn_fwb == 3 ) WRITE(numout,*) ' fwf set to zero and spread out over erp area' 86 IF( kn_fwb == 4 ) WRITE(numout,*) ' instantaneously set to zero with heat and salt flux correction (ISOMIP+)' 86 87 ENDIF 87 88 ! … … 120 121 zcoef = z_fwf * rcp 121 122 emp(:,:) = emp(:,:) - z_fwf * tmask(:,:,1) ! (Eq. 34 AD2015) 122 qns(:,:) = qns(:,:) + zcoef * ( -1.9 ) * tmask(:,:,1) ! (Eq. 35 AD2015) ! could be sst_m if we don't want any bouyancy fluxes 123 sfx(:,:) = sfx(:,:) + z_fwf * ( 33.8 ) * tmask(:,:,1) ! (Eq. 36 AD2015) ! could be sss_m if we don't want any bouyancy fluxes 123 qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) * tmask(:,:,1) ! (Eq. 35 AD2015) ! could be sst_m if we don't want any bouyancy fluxes 124 sfx(:,:) = sfx(:,:) + z_fwf * sss_m(:,:) * tmask(:,:,1) ! (Eq. 36 AD2015) ! could be sss_m if we don't want any bouyancy fluxes 125 !qns(:,:) = qns(:,:) + zcoef * ( -1.9 ) * tmask(:,:,1) ! (Eq. 35 AD2015) ! could be sst_m if we don't want any bouyancy fluxes 126 !sfx(:,:) = sfx(:,:) + z_fwf * ( 33.8 ) * tmask(:,:,1) ! (Eq. 36 AD2015) ! could be sss_m if we don't want any bouyancy fluxes 124 127 ENDIF 125 128 !
Note: See TracChangeset
for help on using the changeset viewer.