Changeset 9118
- Timestamp:
- 2017-12-18T16:17:44+01:00 (5 years ago)
- Location:
- branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3
- Files:
-
- 3 edited
-
icestp.F90 (modified) (2 diffs)
-
icevar.F90 (modified) (3 diffs)
-
icewri.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90
r9076 r9118 235 235 CALL ice_itd_init ! ice thickness distribution initialization 236 236 ! 237 IF( ln_icethd ) THEN 238 CALL ice_thd_init ! set ice thermodynics parameters (clem: important to call it first for melt ponds) 239 ENDIF 237 CALL ice_thd_init ! set ice thermodynics parameters (clem: important to call it first for melt ponds) 238 ! 240 239 ! ! Initial sea-ice state 241 240 IF( .NOT. ln_rstart ) THEN ! start from rest: sea-ice deduced from sst … … 250 249 CALL ice_forcing_init ! set ice-ocean and ice-atm. coupling parameters 251 250 ! 252 IF( ln_icedyn ) THEN 253 CALL ice_dyn_init ! set ice dynamics parameters 254 ENDIF 251 CALL ice_dyn_init ! set ice dynamics parameters 255 252 ! 256 253 CALL ice_update_init ! ice surface boundary condition -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90
r9019 r9118 124 124 tm_si(:,:) = SUM( t_si(:,:,:) * a_i(:,:,:) , dim=3 ) * z1_at_i(:,:) 125 125 om_i (:,:) = SUM( oa_i(:,:,:) , dim=3 ) * z1_at_i(:,:) 126 sm_i (:,:) = SUM( sv_i(:,:,:) , dim=3 ) * z1_vt_i(:,:) 126 127 ! 127 128 tm_i(:,:) = 0._wp 128 sm_i(:,:) = 0._wp129 129 DO jl = 1, jpl 130 130 DO jk = 1, nlay_i 131 131 tm_i(:,:) = tm_i(:,:) + r1_nlay_i * t_i (:,:,jk,jl) * v_i(:,:,jl) * z1_vt_i(:,:) 132 sm_i(:,:) = sm_i(:,:) + r1_nlay_i * sz_i(:,:,jk,jl) * v_i(:,:,jl) * z1_vt_i(:,:)133 132 END DO 134 133 END DO … … 298 297 ! !---------------------------------------! 299 298 sz_i(:,:,:,:) = rn_icesal 300 s_i (:,:,:) = rn_icesal299 s_i (:,:,:) = rn_icesal 301 300 ! 302 301 ! !---------------------------------------------! … … 313 312 ! ! Slope of the linear profile 314 313 WHERE( h_i(:,:,:) > epsi20 ) ; z_slope_s(:,:,:) = 2._wp * s_i(:,:,:) / h_i(:,:,:) 315 ELSEWHERE 314 ELSEWHERE ; z_slope_s(:,:,:) = 0._wp 316 315 END WHERE 317 316 ! -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icewri.F90
r9019 r9118 248 248 CALL histvert( kid, "ncatice", "Ice Categories","", jpl, jcat, nz_i, "up") 249 249 250 CALL histdef( kid, "sithic", "Ice thickness" , "m" , & 251 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 252 CALL histdef( kid, "siconc", "Ice concentration" , "%" , & 253 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 254 CALL histdef( kid, "sitemp", "Ice temperature" , "C" , & 255 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 256 CALL histdef( kid, "sivelu", "i-Ice speed " , "m/s" , & 257 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 258 CALL histdef( kid, "sivelv", "j-Ice speed " , "m/s" , & 259 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 260 CALL histdef( kid, "sistru", "i-Wind stress over ice " , "Pa" , & 261 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 262 CALL histdef( kid, "sistrv", "j-Wind stress over ice " , "Pa" , & 263 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 264 CALL histdef( kid, "sisflx", "Solar flux over ocean" , "w/m2" , & 265 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 266 CALL histdef( kid, "sinflx", "Non-solar flux over ocean" , "w/m2" , & 267 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 268 CALL histdef( kid, "isnowpre", "Snow precipitation" , "kg/m2/s", & 269 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 270 CALL histdef( kid, "sisali", "Ice salinity" , "PSU" , & 271 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 272 CALL histdef( kid, "sivolu", "Ice volume" , "m" , & 273 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 274 CALL histdef( kid, "sidive", "Ice divergence" , "10-8s-1", & 275 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 276 CALL histdef( kid, "si_amp", "Melt pond fraction" , "%" , & 277 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 278 CALL histdef( kid, "si_vmp", "Melt pond volume" , "m" , & 279 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 280 CALL histdef( kid, "vfxbog", "Ice bottom production" , "m/s" , & 281 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 282 CALL histdef( kid, "vfxdyn", "Ice dynamic production" , "m/s" , & 283 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 284 CALL histdef( kid, "vfxopw", "Ice open water prod" , "m/s" , & 285 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 286 CALL histdef( kid, "vfxsni", "Snow ice production " , "m/s" , & 287 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 288 CALL histdef( kid, "vfxres", "Ice prod from corrections" , "m/s" , & 289 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 290 CALL histdef( kid, "vfxbom", "Ice bottom melt" , "m/s" , & 291 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 292 CALL histdef( kid, "vfxsum", "Ice surface melt" , "m/s" , & 293 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 294 CALL histdef( kid, "vfxpnd", "Ice melt ponds flux" , "m/s" , & 295 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 296 297 CALL histdef( kid, "sithicat", "Ice thickness" , "m" , & 298 & jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 299 CALL histdef( kid, "siconcat", "Ice concentration" , "%" , & 300 & jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 301 CALL histdef( kid, "sisalcat", "Ice salinity" , "" , & 302 & jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 303 CALL histdef( kid, "snthicat", "Snw thickness" , "m" , & 304 & jpi, jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 250 CALL histdef( kid, "sithic", "Ice thickness" , "m" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 251 CALL histdef( kid, "siconc", "Ice concentration" , "%" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 252 CALL histdef( kid, "sitemp", "Ice temperature" , "C" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 253 CALL histdef( kid, "sivelu", "i-Ice speed " , "m/s" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 254 CALL histdef( kid, "sivelv", "j-Ice speed " , "m/s" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 255 CALL histdef( kid, "sistru", "i-Wind stress over ice" , "Pa" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 256 CALL histdef( kid, "sistrv", "j-Wind stress over ice" , "Pa" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 257 CALL histdef( kid, "sisflx", "Solar flx over ocean" , "W/m2" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 258 CALL histdef( kid, "sinflx", "NonSolar flx over ocean", "W/m2" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 259 CALL histdef( kid, "snwpre", "Snow precipitation" , "kg/m2/s", jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 260 CALL histdef( kid, "sisali", "Ice salinity" , "PSU" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 261 CALL histdef( kid, "sivolu", "Ice volume" , "m" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 262 CALL histdef( kid, "sidive", "Ice divergence" , "10-8s-1", jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 263 CALL histdef( kid, "si_amp", "Melt pond fraction" , "%" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 264 CALL histdef( kid, "si_vmp", "Melt pond volume" , "m" , jpi,jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 265 ! 266 CALL histdef( kid, "sithicat", "Ice thickness" , "m" , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 267 CALL histdef( kid, "siconcat", "Ice concentration" , "%" , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 268 CALL histdef( kid, "sisalcat", "Ice salinity" , "" , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 269 CALL histdef( kid, "snthicat", "Snw thickness" , "m" , jpi,jpj, kh_i, jpl, 1, jpl, nz_i, 32, "inst(x)", rdt, rdt ) 305 270 306 271 CALL histend( kid, snc4set ) ! end of the file definition 307 272 308 CALL histwrite( kid, "sithic", kt, hm_i , jpi*jpj, (/1/) ) 309 CALL histwrite( kid, "siconc", kt, at_i , jpi*jpj, (/1/) ) 310 CALL histwrite( kid, "sitemp", kt, tm_i - rt0 , jpi*jpj, (/1/) ) 311 CALL histwrite( kid, "sivelu", kt, u_ice , jpi*jpj, (/1/) ) 312 CALL histwrite( kid, "sivelv", kt, v_ice , jpi*jpj, (/1/) ) 313 CALL histwrite( kid, "sistru", kt, utau_ice , jpi*jpj, (/1/) ) 314 CALL histwrite( kid, "sistrv", kt, vtau_ice , jpi*jpj, (/1/) ) 315 CALL histwrite( kid, "sisflx", kt, qsr , jpi*jpj, (/1/) ) 316 CALL histwrite( kid, "sinflx", kt, qns , jpi*jpj, (/1/) ) 317 CALL histwrite( kid, "isnowpre", kt, sprecip , jpi*jpj, (/1/) ) 318 CALL histwrite( kid, "sisali", kt, sm_i , jpi*jpj, (/1/) ) 319 CALL histwrite( kid, "sivolu", kt, vt_i , jpi*jpj, (/1/) ) 320 CALL histwrite( kid, "sidive", kt, divu_i*1.0e8 , jpi*jpj, (/1/) ) 321 273 CALL histwrite( kid, "sithic", kt, hm_i , jpi*jpj, (/1/) ) 274 CALL histwrite( kid, "siconc", kt, at_i , jpi*jpj, (/1/) ) 275 CALL histwrite( kid, "sitemp", kt, tm_i - rt0 , jpi*jpj, (/1/) ) 276 CALL histwrite( kid, "sivelu", kt, u_ice , jpi*jpj, (/1/) ) 277 CALL histwrite( kid, "sivelv", kt, v_ice , jpi*jpj, (/1/) ) 278 CALL histwrite( kid, "sistru", kt, utau_ice , jpi*jpj, (/1/) ) 279 CALL histwrite( kid, "sistrv", kt, vtau_ice , jpi*jpj, (/1/) ) 280 CALL histwrite( kid, "sisflx", kt, qsr , jpi*jpj, (/1/) ) 281 CALL histwrite( kid, "sinflx", kt, qns , jpi*jpj, (/1/) ) 282 CALL histwrite( kid, "snwpre", kt, sprecip , jpi*jpj, (/1/) ) 283 CALL histwrite( kid, "sisali", kt, sm_i , jpi*jpj, (/1/) ) 284 CALL histwrite( kid, "sivolu", kt, vt_i , jpi*jpj, (/1/) ) 285 CALL histwrite( kid, "sidive", kt, divu_i*1.0e8 , jpi*jpj, (/1/) ) 322 286 CALL histwrite( kid, "si_amp", kt, at_ip , jpi*jpj, (/1/) ) 323 287 CALL histwrite( kid, "si_vmp", kt, vt_ip , jpi*jpj, (/1/) ) 324 325 CALL histwrite( kid, "vfxbog", kt, wfx_bog , jpi*jpj, (/1/) ) 326 CALL histwrite( kid, "vfxdyn", kt, wfx_dyn , jpi*jpj, (/1/) ) 327 CALL histwrite( kid, "vfxopw", kt, wfx_opw , jpi*jpj, (/1/) ) 328 CALL histwrite( kid, "vfxsni", kt, wfx_sni , jpi*jpj, (/1/) ) 329 CALL histwrite( kid, "vfxres", kt, wfx_res , jpi*jpj, (/1/) ) 330 CALL histwrite( kid, "vfxbom", kt, wfx_bom , jpi*jpj, (/1/) ) 331 CALL histwrite( kid, "vfxsum", kt, wfx_sum , jpi*jpj, (/1/) ) 332 CALL histwrite( kid, "vfxpnd", kt, wfx_pnd , jpi*jpj, (/1/) ) 333 288 ! 334 289 CALL histwrite( kid, "sithicat", kt, h_i , jpi*jpj*jpl, (/1/) ) 335 290 CALL histwrite( kid, "siconcat", kt, a_i , jpi*jpj*jpl, (/1/) )
Note: See TracChangeset
for help on using the changeset viewer.
