Changeset 4907 for branches/2014/dev_CNRS_2014
- Timestamp:
- 2014-11-27T22:02:56+01:00 (10 years ago)
- Location:
- branches/2014/dev_CNRS_2014/NEMOGCM
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_CNRS_2014/NEMOGCM/CONFIG/SHARED/field_def.xml
r4906 r4907 30 30 <field id="heatc" long_name="Heat content vertically integrated" unit="J/m2" /> 31 31 <field id="saltc" long_name="Salt content vertically integrated" unit="PSU*kg/m2" /> 32 <!-- EOS -- !>32 <!-- EOS --> 33 33 <field id="alpha" long_name="thermal expansion" unit="1/degC" grid_ref="grid_T_3D"/> 34 34 <field id="beta" long_name="haline contraction" unit="1/psu" grid_ref="grid_T_3D"/> 35 35 <field id="bn2" long_name="squared Brunt-Vaisala frequency" unit="1/s" grid_ref="grid_T_3D"/> 36 36 <field id="rhop" long_name="potential density (sigma0)" unit="kg/m3" grid_ref="grid_T_3D"/> 37 <!-- Energy - horizontal divergence -- !>37 <!-- Energy - horizontal divergence --> 38 38 <field id="eken" long_name="kinetic energy" unit="m2/s2" grid_ref="grid_T_3D"/> 39 39 <field id="hdiv" long_name="horizontal divergence" unit="s-1" grid_ref="grid_T_3D"/> -
branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/LIM_SRC_3/limdiahsb.F90
r4902 r4907 75 75 76 76 ! 1/area 77 z1_area = 1. d0 / MAX( glob_sum( area(:,:) * tms(:,:) ), epsi06 )78 79 zinda = MAX( 0. d0 , SIGN( 1.d0 , glob_sum( area(:,:) * tms(:,:) ) - epsi06 ) )77 z1_area = 1.0 / MAX( glob_sum( area(:,:) * tms(:,:) ), epsi06 ) 78 79 zinda = MAX( 0.0 , SIGN( 1.0 , glob_sum( area(:,:) * tms(:,:) ) - epsi06 ) ) 80 80 ! ----------------------- ! 81 81 ! 1 - Content variations ! … … 153 153 ! 3 - Diagnostics writing ! 154 154 ! ----------------------- ! 155 zindb = MAX( 0. d0 , SIGN( 1.d0 , zbg_ivo - epsi06 ) )155 zindb = MAX( 0.0 , SIGN( 1.0 , zbg_ivo - epsi06 ) ) 156 156 ! 157 157 CALL iom_put( 'ibgvoltot' , zbg_ivo * rhoic * r1_rau0 * 1.e-9 ) ! ice volume (km3 equivalent liquid) -
branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r4902 r4907 222 222 ! Second step in limthd_dh : heat remaining if total melt (zq_rema) 223 223 ! Third step in limsbc : heat from ice-ocean mass exchange (zf_mass) + solar 224 hfx_out(ji,jj) = hfx_out(ji,jj) 224 hfx_out(ji,jj) = hfx_out(ji,jj) & 225 225 ! Non solar heat flux received by the ocean 226 & + pfrld(ji,jj) * qns(ji,jj) 226 & + pfrld(ji,jj) * qns(ji,jj) & 227 227 ! latent heat of precip (note that precip is included in qns but not in qns_ice) 228 & + ( pfrld(ji,jj)**betas - pfrld(ji,jj) ) * sprecip(ji,jj) * ( cpic * ( MIN( tatm_ice(ji,jj), rt0_snow ) - rtt ) - lfus ) & 229 & + ( 1._wp - pfrld(ji,jj) ) * ( tprecip(ji,jj) - sprecip(ji,jj) ) * rcp * ( tatm_ice(ji,jj) - rtt ) & 228 & + ( pfrld(ji,jj)**betas - pfrld(ji,jj) ) * sprecip(ji,jj) & 229 & * ( cpic * ( MIN( tatm_ice(ji,jj), rt0_snow ) - rtt ) - lfus ) & 230 & + ( 1._wp - pfrld(ji,jj) ) * ( tprecip(ji,jj) - sprecip(ji,jj) ) * rcp * ( tatm_ice(ji,jj) - rtt ) & 230 231 ! heat flux taken from the ocean where there is open water ice formation 231 & - qlead(ji,jj) * r1_rdtice 232 & - qlead(ji,jj) * r1_rdtice & 232 233 ! heat flux taken from the ocean during bottom growth/melt (fhld should be 0 while bott growth) 233 & - at_i(ji,jj) * fhtur(ji,jj) 234 & - at_i(ji,jj) * fhtur(ji,jj) & 234 235 & - at_i(ji,jj) * fhld(ji,jj) 235 236 -
branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/LIM_SRC_3/limthd_lac.F90
r4902 r4907 133 133 !Energy of melting q(S,T) [J.m-3] 134 134 zindb = 1._wp - MAX( 0._wp , SIGN( 1._wp , -v_i(ji,jj,jl) + epsi10 ) ) !0 if no ice and 1 if yes 135 e_i(ji,jj,jk,jl) = zindb * e_i(ji,jj,jk,jl) / ( area(ji,jj) * MAX( v_i(ji,jj,jl) , epsi10 ) ) * REAL( nlay_i, wp ) 135 e_i(ji,jj,jk,jl) = zindb * e_i(ji,jj,jk,jl) & 136 & / ( area(ji,jj) * MAX( v_i(ji,jj,jl) , epsi10 ) ) * REAL( nlay_i, wp ) 136 137 e_i(ji,jj,jk,jl) = e_i(ji,jj,jk,jl) * unit_fac 137 138 END DO
Note: See TracChangeset
for help on using the changeset viewer.