- Timestamp:
- 2013-07-09T17:41:20+02:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r3406_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90
r3938 r3963 122 122 ! Normal file 123 123 !------------- 124 125 zsto = rdt_ice126 IF( ln_mskland ) THEN ; clop = "ave(only(x))" ! put 1.e+20 on land (very expensive!!)127 ELSE ; clop = "ave(x)" ! no use of the mask value (require less cpu time)128 ENDIF129 zout = nwrite * rdt_ice / nn_fsbc130 124 niter = ( nit000 - 1 ) / nn_fsbc 131 zdept(1) = 0.132 133 125 CALL ymds2ju ( nyear, nmonth, nday, rdt, zjulian ) 134 126 zjulian = zjulian - adatrj ! set calendar origin to the beginning of the experiment 135 CALL dia_nam ( clhstnam, nwrite, 'icemod_old' ) 136 CALL histbeg ( clhstnam, jpi, glamt, jpj, gphit, 1, jpi, 1, jpj, niter, zjulian, rdt_ice, & 137 & nhorid, nice, domain_id=nidom, snc4chunks=snc4set ) 138 CALL histvert( nice, "deptht", "Vertical T levels", "m", 1, zdept, ndepid, "down") 139 CALL wheneq ( jpij , tmask(:,:,1), 1, 1., ndex51, ndim) 140 141 DO jf = 1 , noumef 142 IF(lwp) WRITE(numout,*) 'jf', jf 143 IF ( nc(jf) == 1 ) THEN 144 CALL histdef( nice, nam(jf), titn(jf), uni(jf), jpi, jpj & 145 , nhorid, 1, 1, 1, -99, 32, clop, zsto, zout ) 146 IF(lwp) WRITE(numout,*) 'nice, nam(jf), titn(jf), uni(jf), nhorid, clop, zsto, zout' 147 IF(lwp) WRITE(numout,*) nice, nam(jf), titn(jf), uni(jf), nhorid, clop, zsto, zout 148 ENDIF 149 END DO 150 151 CALL histend(nice, snc4set) 152 127 !clem 128 ! zsto = rdt_ice 129 ! IF( ln_mskland ) THEN ; clop = "ave(only(x))" ! put 1.e+20 on land (very expensive!!) 130 ! ELSE ; clop = "ave(x)" ! no use of the mask value (require less cpu time) 131 ! ENDIF 132 ! zout = nwrite * rdt_ice / nn_fsbc 133 ! zdept(1) = 0. 134 ! 135 ! CALL dia_nam ( clhstnam, nwrite, 'icemod_old' ) 136 ! CALL histbeg ( clhstnam, jpi, glamt, jpj, gphit, 1, jpi, 1, jpj, niter, zjulian, rdt_ice, & 137 ! & nhorid, nice, domain_id=nidom, snc4chunks=snc4set ) 138 ! CALL histvert( nice, "deptht", "Vertical T levels", "m", 1, zdept, ndepid, "down") 139 ! CALL wheneq ( jpij , tmask(:,:,1), 1, 1., ndex51, ndim) 140 ! 141 ! DO jf = 1 , noumef 142 ! IF(lwp) WRITE(numout,*) 'jf', jf 143 ! IF ( nc(jf) == 1 ) THEN 144 ! CALL histdef( nice, nam(jf), titn(jf), uni(jf), jpi, jpj & 145 ! , nhorid, 1, 1, 1, -99, 32, clop, zsto, zout ) 146 ! IF(lwp) WRITE(numout,*) 'nice, nam(jf), titn(jf), uni(jf), nhorid, clop, zsto, zout' 147 ! IF(lwp) WRITE(numout,*) nice, nam(jf), titn(jf), uni(jf), nhorid, clop, zsto, zout 148 ! ENDIF 149 ! END DO 150 ! 151 ! CALL histend(nice, snc4set) 152 !clem 153 ! 153 154 !----------------- 154 155 ! ITD file output … … 280 281 ! ecriture d'un fichier netcdf 281 282 ! 283 282 284 niter = niter + 1 283 DO jf = 1 , noumef 284 ! 285 zfield(:,:) = zcmo(:,:,jf) * cmulti(jf) + cadd(jf) 286 ! 287 IF( jf == 7 .OR. jf == 8 .OR. jf == 15 .OR. jf == 16 ) THEN ; CALL lbc_lnk( zfield, 'T', -1. ) 288 ELSE ; CALL lbc_lnk( zfield, 'T', 1. ) 289 ENDIF 290 ! 291 IF( ln_nicep ) THEN 292 WRITE(numout,*) 293 WRITE(numout,*) 'nc(jf), nice, nam(jf), niter, ndim' 294 WRITE(numout,*) nc(jf), nice, nam(jf), niter, ndim 295 ENDIF 296 IF( nc(jf) == 1 ) CALL histwrite( nice, nam(jf), niter, zfield, ndim, ndex51 ) 297 ! 298 END DO 299 300 IF( ( nn_fsbc * niter ) >= nitend .OR. kindic < 0 ) THEN 301 IF( lwp) WRITE(numout,*) ' Closing the icemod file ' 302 CALL histclo( nice ) 303 ENDIF 304 285 !clem 286 ! DO jf = 1 , noumef 287 ! ! 288 ! zfield(:,:) = zcmo(:,:,jf) * cmulti(jf) + cadd(jf) 289 ! ! 290 ! IF( jf == 7 .OR. jf == 8 .OR. jf == 15 .OR. jf == 16 ) THEN ; CALL lbc_lnk( zfield, 'T', -1. ) 291 ! ELSE ; CALL lbc_lnk( zfield, 'T', 1. ) 292 ! ENDIF 293 ! ! 294 ! IF( ln_nicep ) THEN 295 ! WRITE(numout,*) 296 ! WRITE(numout,*) 'nc(jf), nice, nam(jf), niter, ndim' 297 ! WRITE(numout,*) nc(jf), nice, nam(jf), niter, ndim 298 ! ENDIF 299 ! IF( nc(jf) == 1 ) CALL histwrite( nice, nam(jf), niter, zfield, ndim, ndex51 ) 300 ! ! 301 ! END DO 302 ! 303 ! IF( ( nn_fsbc * niter ) >= nitend .OR. kindic < 0 ) THEN 304 ! IF( lwp) WRITE(numout,*) ' Closing the icemod file ' 305 ! CALL histclo( nice ) 306 ! ENDIF 307 !clem 308 ! 305 309 CALL iom_put ('iceconc', zcmo(:,:,1) ) ! field1: ice concentration 306 310 CALL iom_put ('icethic_cea', zcmo(:,:,2) ) ! field2: ice thickness (i.e. icethi(:,:)) … … 328 332 CALL iom_put ('isnowhc', zcmo(:,:,26) ) ! field 26: snow total heat content 329 333 CALL iom_put ('icest', zcmo(:,:,27) ) ! field 27: ice surface temperature 330 CALL iom_put ('fsbri', zcmo(:,:,28) * 86400 ) ! field 28: brine salt flux331 CALL iom_put ('fseqv', zcmo(:,:,29) * 86400 ) ! field 29: equivalent FW salt flux334 CALL iom_put ('fsbri', zcmo(:,:,28) * 86400. ) ! field 28: brine salt flux 335 CALL iom_put ('fseqv', zcmo(:,:,29) * 86400. ) ! field 29: equivalent FW salt flux 332 336 CALL iom_put ('ibrinv', zcmo(:,:,30) *100 ) ! field 30: brine volume 333 337 CALL iom_put ('icecolf', zcmo(:,:,31) ) ! field 31: frazil ice collection thickness … … 341 345 CALL iom_put ('icevolu', zcmo(:,:,39) ) ! field 39: ice volume 342 346 CALL iom_put ('snowvol', zcmo(:,:,40) ) ! field 40: snow volume 343 CALL iom_put ('fsrpo', zcmo(:,:,41) * 86400 ) ! field 41: salt flux from ridging rafting344 CALL iom_put ('fsres', zcmo(:,:,42) * 86400 ) ! field 42: salt flux from limupdate (resultant)347 CALL iom_put ('fsrpo', zcmo(:,:,41) * 86400. ) ! field 41: salt flux from ridging rafting 348 CALL iom_put ('fsres', zcmo(:,:,42) * 86400. ) ! field 42: salt flux from limupdate (resultant) 345 349 CALL iom_put ('icetrp', zcmo(:,:,43) ) ! field 43: ice volume transport 346 350
Note: See TracChangeset
for help on using the changeset viewer.