Changeset 14273 for NEMO/trunk/tests
- Timestamp:
- 2021-01-06T11:57:45+01:00 (3 years ago)
- Location:
- NEMO/trunk/tests
- Files:
-
- 1 added
- 7 edited
-
BENCH/MY_SRC/usrdef_istate.F90 (modified) (1 diff)
-
BENCH/MY_SRC/usrdef_sbc.F90 (modified) (9 diffs)
-
DONUT/EXPREF/file_def_nemo-ice.xml (modified) (3 diffs)
-
DONUT/MY_SRC/usrdef_sbc.F90 (added)
-
ICE_ADV1D/MY_SRC/usrdef_sbc.F90 (modified) (1 diff)
-
ICE_ADV2D/MY_SRC/usrdef_sbc.F90 (modified) (1 diff)
-
ICE_AGRIF/MY_SRC/usrdef_sbc.F90 (modified) (1 diff)
-
ICE_RHEO/MY_SRC/usrdef_sbc.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_istate.F90
r14223 r14273 72 72 zfact = REAL(jk-1,wp) / REAL(jpk-1,wp) ! 0 to 1 to add a basic stratification 73 73 ! temperature choosen to lead to ~50% ice at the beginning if rn_thres_sst = 0.5 74 pts( :,:,jk,jp_tem) = 20._wp*z2d(:,:) - 1._wp - 0.5_wp * zfact ! -1 to -1.5 +/- 1.0 degG74 pts(ji,jj,jk,jp_tem) = 20._wp*z2d(ji,jj) - 1._wp - 0.5_wp * zfact ! -1 to -1.5 +/- 1.0 degG 75 75 ! salinity: 76 pts( :,:,jk,jp_sal) = 30._wp + 1._wp * zfact + z2d(:,:) ! 30 to 31 +/- 0.05 psu76 pts(ji,jj,jk,jp_sal) = 30._wp + 1._wp * zfact + z2d(ji,jj) ! 30 to 31 +/- 0.05 psu 77 77 ! velocities: 78 pu( :,:,jk) = z2d(:,:) * 0.1_wp * umask(:,:,jk)! +/- 0.005 m/s79 pv( :,:,jk) = z2d(:,:) * 0.01_wp * vmask(:,:,jk)! +/- 0.0005 m/s78 pu(ji,jj,jk) = z2d(ji,jj) * 0.1_wp * umask(ji,jj,jk) ! +/- 0.005 m/s 79 pv(ji,jj,jk) = z2d(ji,jj) * 0.01_wp * vmask(ji,jj,jk) ! +/- 0.0005 m/s 80 80 END_3D 81 81 pts(:,:,jpk,:) = 0._wp -
NEMO/trunk/tests/BENCH/MY_SRC/usrdef_sbc.F90
r14223 r14273 30 30 PRIVATE 31 31 32 PUBLIC usrdef_sbc_oce ! routine called in sbcmod module33 PUBLIC usrdef_sbc_ice_tau ! routine called by sbcice_lim.F90 for ice dynamics34 PUBLIC usrdef_sbc_ice_flx ! routine called by sbcice_lim.F90 for ice thermo32 PUBLIC usrdef_sbc_oce ! routine called by sbcmod.F90 for sbc ocean 33 PUBLIC usrdef_sbc_ice_tau ! routine called by icestp.F90 for ice dynamics 34 PUBLIC usrdef_sbc_ice_flx ! routine called by icestp.F90 for ice thermo 35 35 36 36 !! * Substitutions … … 83 83 ! 84 84 ENDIF 85 86 85 ! 87 86 END SUBROUTINE usrdef_sbc_oce … … 97 96 INTEGER, INTENT(in) :: kt ! ocean time step 98 97 ! 99 REAL(wp) , DIMENSION(jpi,jpj) :: z2d ! 2D workspace98 REAL(wp) :: zztmp 100 99 INTEGER :: ji, jj 101 100 !!--------------------------------------------------------------------- … … 106 105 ! 107 106 DO_2D( 0, 0, 0, 0 ) 108 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 107 zztmp = 0.1 * ( 0.5 - REAL( mig0(ji) + (mjg0(jj)-1) * Ni0glo, wp ) / REAL( Ni0glo * Nj0glo, wp ) ) 108 utau_ice(ji,jj) = 0.1_wp + zztmp 109 vtau_ice(ji,jj) = 0.1_wp + zztmp 109 110 END_2D 110 utau_ice(:,:) = 0.1_wp + z2d(:,:)111 vtau_ice(:,:) = 0.1_wp + z2d(:,:)112 111 113 112 CALL lbc_lnk_multi( 'usrdef_sbc', utau_ice, 'U', -1., vtau_ice, 'V', -1. ) … … 121 120 !! *** ROUTINE usrdef_sbc_ice_flx *** 122 121 !! 123 !! ** Purpose : provide the surface boundary (flux) condition over 124 !sea-ice 122 !! ** Purpose : provide the surface boundary (flux) condition over sea-ice 125 123 !!--------------------------------------------------------------------- 126 124 INTEGER, INTENT(in) :: kt ! ocean time step … … 130 128 REAL(wp), DIMENSION(jpi,jpj) :: zsnw ! snw distribution after wind blowing 131 129 !!--------------------------------------------------------------------- 130 #if defined key_si3 132 131 ! 133 #if defined key_si3134 132 IF( kt==nit000 .AND. lwp) WRITE(numout,*)' usrdef_sbc_ice : BENCH case: NO flux forcing' 135 133 ! … … 137 135 emp_oce (:,:) = 0._wp ! uniform value for freshwater budget (E-P) 138 136 qsr_oce (:,:) = 0._wp ! uniform value for solar radiation 139 qns_oce (:,:) = 0._wp ! uniform value for non-solar radiation137 qns_oce (:,:) = 0._wp ! uniform value for non-solar heat flux 140 138 141 139 ! ice variables 142 140 alb_ice (:,:,:) = 0.7_wp ! useless 143 141 qsr_ice (:,:,:) = 0._wp ! uniform value for solar radiation 144 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar radiation 142 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar heat flux 143 dqns_ice(:,:,:) = 0._wp ! uniform value for non solar heat flux sensitivity for ice 145 144 sprecip (:,:) = 0._wp ! uniform value for snow precip 146 145 evap_ice(:,:,:) = 0._wp ! uniform value for sublimation … … 148 147 ! ice fields deduced from above 149 148 zsnw(:,:) = 1._wp 150 !!CALL lim_thd_snwblow( at_i_b, zsnw ) ! snow distribution over ice after 151 !wind blowing 149 !!CALL lim_thd_snwblow( at_i_b, zsnw ) ! snow distribution over ice after wind blowing 152 150 emp_ice (:,:) = SUM( a_i_b(:,:,:) * evap_ice(:,:,:), dim=3 ) - sprecip(:,:) * zsnw(:,:) 153 151 emp_oce (:,:) = emp_oce(:,:) - sprecip(:,:) * (1._wp - zsnw(:,:) ) … … 164 162 ! --- shortwave radiation transmitted thru the surface scattering layer (W/m2) --- ! 165 163 qtr_ice_top(:,:,:) = 0._wp 166 167 164 #endif 168 165 -
NEMO/trunk/tests/DONUT/EXPREF/file_def_nemo-ice.xml
r14226 r14273 53 53 <field field_ref="normstr" name="normstr" /> 54 54 <field field_ref="sheastr" name="sheastr" /> 55 <field field_ref="isig1" name="isig1" /> 56 <field field_ref="isig2" name="isig2" /> 57 <field field_ref="isig3" name="isig3" /> 55 <field field_ref="sig1_pnorm" name="sig1_pnorm"/> 56 <field field_ref="sig2_pnorm" name="sig2_pnorm"/> 58 57 59 58 <!-- heat fluxes --> … … 68 67 <field field_ref="albedo" name="albedo" /> 69 68 70 <!-- heat fluxes --> 71 <field id="hfxbog" long_name="heat flux used for bottom ice growth (neg.)" unit="W/m2" /> 72 <field id="hfxbom" long_name="heat flux used for bottom ice melt (pos.)" unit="W/m2" /> 73 <field id="hfxsum" long_name="heat flux used for surface ice melt" unit="W/m2" /> 74 <field id="hfxopw" long_name="heat flux used for open water ice formation" unit="W/m2" /> 75 <field id="hfxdif" long_name="heat flux used for ice temperature change" unit="W/m2" /> 76 <field id="hfxsnw" long_name="heat flux used for snow melt" unit="W/m2" /> 77 <field id="hfxerr" long_name="heat flux error after heat diffusion" unit="W/m2" /> 78 <!-- heat fluxes associated with mass exchange --> 79 <field id="hfxthd" long_name="heat fluxes from ice-ocean mass exchange during thermo" unit="W/m2" /> 80 <field id="hfxdyn" long_name="heat fluxes from ice-ocean mass exchange during dynamic" unit="W/m2" /> 81 <field id="hfxres" long_name="heat fluxes from undiagnosed processes" unit="W/m2" /> 82 <field id="hfxsub" long_name="heat fluxes from ice-atm. mass exchange during sublimation" unit="W/m2" /> 83 <field id="hfxspr" long_name="heat fluxes from ice-atm. mass exchange during snow precip" unit="W/m2" /> 84 <field id="hfxdhc" long_name="Heat content variation in snow and ice (neg = ice cooling)" unit="W/m2" /> 85 86 <field field_ref="hfxcndtop" name="hfxcndtop" /> 69 <field field_ref="hfxcndtop" name="hfxcndtop" /> 87 70 <field field_ref="hfxcndbot" name="hfxcndbot" /> 88 71 <field field_ref="hfxsensib" name="hfxsensib" /> … … 94 77 <field field_ref="vfxice" name="vfxice" /> 95 78 <field field_ref="vfxsnw" name="vfxsnw" /> 96 79 97 80 <!-- categories --> 98 81 <field field_ref="icemask_cat" name="simskcat"/> -
NEMO/trunk/tests/ICE_ADV1D/MY_SRC/usrdef_sbc.F90
r14072 r14273 120 120 emp_oce (:,:) = 0._wp ! uniform value for freshwater budget (E-P) 121 121 qsr_oce (:,:) = 0._wp ! uniform value for solar radiation 122 qns_oce (:,:) = 0._wp ! uniform value for non-solar radiation122 qns_oce (:,:) = 0._wp ! uniform value for non-solar heat flux 123 123 124 124 ! ice variables 125 125 alb_ice (:,:,:) = 0.7_wp ! useless 126 126 qsr_ice (:,:,:) = 0._wp ! uniform value for solar radiation 127 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar radiation 127 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar heat flux 128 dqns_ice(:,:,:) = 0._wp ! uniform value for non solar heat flux sensitivity for ice 128 129 sprecip (:,:) = 0._wp ! uniform value for snow precip 129 130 evap_ice(:,:,:) = 0._wp ! uniform value for sublimation -
NEMO/trunk/tests/ICE_ADV2D/MY_SRC/usrdef_sbc.F90
r13999 r14273 118 118 emp_oce (:,:) = 0._wp ! uniform value for freshwater budget (E-P) 119 119 qsr_oce (:,:) = 0._wp ! uniform value for solar radiation 120 qns_oce (:,:) = 0._wp ! uniform value for non-solar radiation120 qns_oce (:,:) = 0._wp ! uniform value for non-solar heat flux 121 121 122 122 ! ice variables 123 123 alb_ice (:,:,:) = 0.7_wp ! useless 124 124 qsr_ice (:,:,:) = 0._wp ! uniform value for solar radiation 125 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar radiation 125 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar heat flux 126 dqns_ice(:,:,:) = 0._wp ! uniform value for non solar heat flux sensitivity for ice 126 127 sprecip (:,:) = 0._wp ! uniform value for snow precip 127 128 evap_ice(:,:,:) = 0._wp ! uniform value for sublimation -
NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_sbc.F90
r14072 r14273 119 119 emp_oce (:,:) = 0._wp ! uniform value for freshwater budget (E-P) 120 120 qsr_oce (:,:) = 0._wp ! uniform value for solar radiation 121 qns_oce (:,:) = 0._wp ! uniform value for non-solar radiation121 qns_oce (:,:) = 0._wp ! uniform value for non-solar heat flux 122 122 123 123 ! ice variables 124 124 alb_ice (:,:,:) = 0.7_wp ! useless 125 125 qsr_ice (:,:,:) = 0._wp ! uniform value for solar radiation 126 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar radiation 126 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar heat flux 127 dqns_ice(:,:,:) = 0._wp ! uniform value for non solar heat flux sensitivity for ice 127 128 sprecip (:,:) = 0._wp ! uniform value for snow precip 128 129 evap_ice(:,:,:) = 0._wp ! uniform value for sublimation -
NEMO/trunk/tests/ICE_RHEO/MY_SRC/usrdef_sbc.F90
r14063 r14273 179 179 emp_oce (:,:) = 0._wp ! uniform value for freshwater budget (E-P) 180 180 qsr_oce (:,:) = 0._wp ! uniform value for solar radiation 181 qns_oce (:,:) = 0._wp ! uniform value for non-solar radiation181 qns_oce (:,:) = 0._wp ! uniform value for non-solar heat flux 182 182 183 183 ! ice variables 184 184 alb_ice (:,:,:) = 0.7_wp ! useless 185 185 qsr_ice (:,:,:) = 0._wp ! uniform value for solar radiation 186 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar radiation 186 qns_ice (:,:,:) = 0._wp ! uniform value for non-solar heat flux 187 dqns_ice(:,:,:) = 0._wp ! uniform value for non solar heat flux sensitivity for ice 187 188 sprecip (:,:) = 0._wp ! uniform value for snow precip 188 189 evap_ice(:,:,:) = 0._wp ! uniform value for sublimation
Note: See TracChangeset
for help on using the changeset viewer.
