- Timestamp:
- 2017-10-18T19:14:32+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/LIM_SRC_3/icewri.F90
r8586 r8637 128 128 CALL iom_put( "snowvol" , vt_s * zswi ) ! snow volume 129 129 130 IF ( ln_pnd ) THEN 131 CALL iom_put( "iceamp" , at_ip * zswi ) ! melt pond total fraction 132 CALL iom_put( "icevmp" , vt_ip * zswi ) ! melt pond total volume per unit area 133 ENDIF 130 CALL iom_put( "iceamp" , at_ip * zswi ) ! melt pond total fraction 131 CALL iom_put( "icevmp" , vt_ip * zswi ) ! melt pond total volume per unit area 134 132 135 133 !---------------------------------- … … 145 143 IF ( iom_use('brinevol_cat') ) CALL iom_put( "brinevol_cat", bv_i * 100. * zswi2 ) ! brine volume 146 144 147 IF ( ln_pnd ) THEN 148 IF ( iom_use('iceamp_cat') ) CALL iom_put( "iceamp_cat" , a_ip * zswi2 ) ! melt pond frac for categories 149 IF ( iom_use('icevmp_cat') ) CALL iom_put( "icevmp_cat" , v_ip * zswi2 ) ! melt pond frac for categories 150 IF ( iom_use('icehmp_cat') ) CALL iom_put( "icehmp_cat" , h_ip * zswi2 ) ! melt pond frac for categories 151 IF ( iom_use('iceafp_cat') ) CALL iom_put( "iceafp_cat" , a_ip_frac * zswi2 ) ! melt pond frac for categories 152 ENDIF 145 IF ( iom_use('iceamp_cat') ) CALL iom_put( "iceamp_cat" , a_ip * zswi2 ) ! melt pond frac for categories 146 IF ( iom_use('icevmp_cat') ) CALL iom_put( "icevmp_cat" , v_ip * zswi2 ) ! melt pond frac for categories 147 IF ( iom_use('icehmp_cat') ) CALL iom_put( "icehmp_cat" , h_ip * zswi2 ) ! melt pond frac for categories 148 IF ( iom_use('iceafp_cat') ) CALL iom_put( "iceafp_cat" , a_ip_frac * zswi2 ) ! melt pond frac for categories 153 149 154 150 !-------------------------------- … … 316 312 CALL histdef( kid, "sidive", "Ice divergence" , "10-8s-1", & 317 313 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 318 319 ! MV MP 2016 320 IF ( ln_pnd ) THEN 321 CALL histdef( kid, "si_amp", "Melt pond fraction" , "%" , & 322 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 323 CALL histdef( kid, "si_vmp", "Melt pond volume" , "m" , & 324 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 325 ENDIF 326 ! END MV MP 2016 327 314 CALL histdef( kid, "si_amp", "Melt pond fraction" , "%" , & 315 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 316 CALL histdef( kid, "si_vmp", "Melt pond volume" , "m" , & 317 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 328 318 CALL histdef( kid, "vfxbog", "Ice bottom production" , "m/s" , & 329 319 & jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) … … 370 360 CALL histwrite( kid, "sidive", kt, divu_i*1.0e8 , jpi*jpj, (/1/) ) 371 361 372 ! MV MP 2016 373 IF ( ln_pnd ) THEN 374 CALL histwrite( kid, "si_amp", kt, at_ip , jpi*jpj, (/1/) ) 375 CALL histwrite( kid, "si_vmp", kt, vt_ip , jpi*jpj, (/1/) ) 376 ENDIF 377 ! END MV MP 2016 362 CALL histwrite( kid, "si_amp", kt, at_ip , jpi*jpj, (/1/) ) 363 CALL histwrite( kid, "si_vmp", kt, vt_ip , jpi*jpj, (/1/) ) 378 364 379 365 CALL histwrite( kid, "vfxbog", kt, wfx_bog , jpi*jpj, (/1/) ) … … 384 370 CALL histwrite( kid, "vfxbom", kt, wfx_bom , jpi*jpj, (/1/) ) 385 371 CALL histwrite( kid, "vfxsum", kt, wfx_sum , jpi*jpj, (/1/) ) 386 IF ( ln_pnd ) & 387 CALL histwrite( kid, "vfxpnd", kt, wfx_pnd , jpi*jpj, (/1/) ) 372 CALL histwrite( kid, "vfxpnd", kt, wfx_pnd , jpi*jpj, (/1/) ) 388 373 389 374 CALL histwrite( kid, "sithicat", kt, h_i , jpi*jpj*jpl, (/1/) )
Note: See TracChangeset
for help on using the changeset viewer.