Changeset 3124
- Timestamp:
- 2011-11-16T10:25:18+01:00 (12 years ago)
- Location:
- branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 18 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/C14b/trcnam_c14b.F90
r2977 r3124 38 38 !! 39 39 !! ** Method : Read the namc14 namelist and check the parameter 40 !! values called at the first timestep (nit 000)40 !! values called at the first timestep (nittrc000) 41 41 !! 42 42 !! ** input : Namelist namelist_c14b -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcexp.F90
r3116 r3124 126 126 ! Time filter and swap of arrays 127 127 ! ------------------------------ 128 IF( neuler == 0 .AND. kt == nit 000 ) THEN ! Euler time-stepping at first time-step128 IF( neuler == 0 .AND. kt == nittrc000 ) THEN ! Euler time-stepping at first time-step 129 129 ! ! (only swap) 130 130 sedpocn(:,:) = sedpoca(:,:) -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zflx.F90
r3028 r3124 252 252 !! 253 253 !! ** Method : Read the nampisext namelist and check the parameters 254 !! called at the first timestep (nit 000)254 !! called at the first timestep (nittrc000) 255 255 !! ** input : Namelist nampisext 256 256 !!---------------------------------------------------------------------- … … 326 326 327 327 ! ! -------------------- ! 328 IF( kt == nit 000 ) THEN ! First call kt=nit000 !328 IF( kt == nittrc000 ) THEN ! First call kt=nittrc000 ! 329 329 ! ! -------------------- ! 330 330 ! !* set file information (default values) -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zlim.F90
r2977 r3124 182 182 !! 183 183 !! ** Method : Read the nampislim namelist and check the parameters 184 !! called at the first timestep (nit 000)184 !! called at the first timestep (nittrc000) 185 185 !! 186 186 !! ** input : Namelist nampislim -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zlys.F90
r2977 r3124 181 181 !! 182 182 !! ** Method : Read the nampiscal namelist and check the parameters 183 !! called at the first timestep (nit 000)183 !! called at the first timestep (nittrc000) 184 184 !! 185 185 !! ** input : Namelist nampiscal -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zmeso.F90
r2977 r3124 237 237 !! 238 238 !! ** Method : Read the nampismes namelist and check the parameters 239 !! called at the first timestep (nit 000)239 !! called at the first timestep (nittrc000) 240 240 !! 241 241 !! ** input : Namelist nampismes -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zmicro.F90
r2977 r3124 201 201 !! 202 202 !! ** Method : Read the nampiszoo namelist and check the parameters 203 !! called at the first timestep (nit 000)203 !! called at the first timestep (nittrc000) 204 204 !! 205 205 !! ** input : Namelist nampiszoo -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zprod.F90
r2979 r3124 466 466 !! 467 467 !! ** Method : Read the nampisprod namelist and check the parameters 468 !! called at the first timestep (nit 000)468 !! called at the first timestep (nittrc000) 469 469 !! 470 470 !! ** input : Namelist nampisprod -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zsed.F90
r2977 r3124 337 337 !!--------------------------------------------------------------------- 338 338 339 ! Compute dust at nit 000 or only if there is more than 1 time record in dust file339 ! Compute dust at nittrc000 or only if there is more than 1 time record in dust file 340 340 IF( ln_dust ) THEN 341 IF( kt == nit 000 .OR. ( kt /= nit000 .AND. ntimes_dust > 1 ) ) THEN341 IF( kt == nittrc000 .OR. ( kt /= nittrc000 .AND. ntimes_dust > 1 ) ) THEN 342 342 CALL fld_read( kt, 1, sf_dust ) 343 343 dust(:,:) = sf_dust(1)%fnow(:,:,1) … … 346 346 347 347 ! N/P and Si releases due to coastal rivers 348 ! Compute river at nit 000 or only if there is more than 1 time record in river file348 ! Compute river at nittrc000 or only if there is more than 1 time record in river file 349 349 ! ----------------------------------------- 350 350 IF( ln_river ) THEN 351 IF( kt == nit 000 .OR. ( kt /= nit000 .AND. ntimes_riv > 1 ) ) THEN351 IF( kt == nittrc000 .OR. ( kt /= nittrc000 .AND. ntimes_riv > 1 ) ) THEN 352 352 CALL fld_read( kt, 1, sf_riverdic ) 353 353 CALL fld_read( kt, 1, sf_riverdoc ) … … 362 362 ENDIF 363 363 364 ! Compute N deposition at nit 000 or only if there is more than 1 time record in N deposition file364 ! Compute N deposition at nittrc000 or only if there is more than 1 time record in N deposition file 365 365 IF( ln_ndepo ) THEN 366 IF( kt == nit 000 .OR. ( kt /= nit000 .AND. ntimes_ndep > 1 ) ) THEN366 IF( kt == nittrc000 .OR. ( kt /= nittrc000 .AND. ntimes_ndep > 1 ) ) THEN 367 367 CALL fld_read( kt, 1, sf_ndepo ) 368 368 DO jj = 1, jpj … … 384 384 !! 385 385 !! ** method : read the files and compute the budget 386 !! called at the first timestep (nit 000)386 !! called at the first timestep (nittrc000) 387 387 !! 388 388 !! ** input : external netcdf files -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/PISCES/trcrst_pisces.F90
r2977 r3124 195 195 ENDIF 196 196 ! 197 CALL trc_dta( nit 000, ztrcdta ) ! read tracer data at nit000197 CALL trc_dta( nittrc000, ztrcdta ) ! read tracer data at nittrc000 198 198 ! 199 199 DO jn = 1, jptra -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv.F90
r3116 r3124 85 85 IF( neuler == 0 .AND. kt == nittrc000 ) THEN ! at nittrc000 86 86 r2dt(:) = rdttrc(:) ! = rdttrc (restarting with Euler time stepping) 87 ELSEIF( kt <= nittrc000 + nn_dttrc) THEN ! at nittrc000 or nittrc000+187 ELSEIF( kt <= nittrc000 + 1 ) THEN ! at nittrc000 or nittrc000+1 88 88 r2dt(:) = 2. * rdttrc(:) ! = 2 rdttrc (leapfrog) 89 89 ENDIF -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90
r3116 r3124 119 119 120 120 ! set time step size (Euler/Leapfrog) 121 IF( neuler == 0 .AND. kt == nittrc000 ) THEN ; r2dt(:) = rdttrc(:) ! at nittrc000 (Euler)122 ELSEIF( kt <= nittrc000 + nn_dttrc )THEN ; r2dt(:) = 2.* rdttrc(:) ! at nit000 or nit000+1 (Leapfrog)121 IF( neuler == 0 .AND. kt == nittrc000 ) THEN ; r2dt(:) = rdttrc(:) ! at nittrc000 (Euler) 122 ELSEIF( kt <= nittrc000 + 1 ) THEN ; r2dt(:) = 2.* rdttrc(:) ! at nit000 or nit000+1 (Leapfrog) 123 123 ENDIF 124 124 … … 139 139 ELSE 140 140 ! Leap-Frog + Asselin filter time stepping 141 IF( lk_vvl ) THEN ; CALL tra_nxt_vvl( kt, nittrc000, 'TRC', trb, trn, tra, jptra ) ! variable volume level (vvl)142 ELSE ; CALL tra_nxt_fix( kt, nittrc000, 'TRC', trb, trn, tra, jptra ) ! fixed volume level141 IF( lk_vvl ) THEN ; CALL tra_nxt_vvl( kt, nittrc000, 'TRC', trb, trn, tra, jptra ) ! variable volume level (vvl) 142 ELSE ; CALL tra_nxt_fix( kt, nittrc000, 'TRC', trb, trn, tra, jptra ) ! fixed volume level 143 143 ENDIF 144 144 ENDIF -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/TRP/trcrad.F90
r3116 r3124 188 188 IF( l_trdtrc ) THEN 189 189 ! 190 zs2rdt = 1. / ( 2. * rdt * FLOAT( nn_dttrc) )190 zs2rdt = 1. / ( 2. * rdt * FLOAT( nn_dttrc ) ) 191 191 ztrtrdb(:,:,:) = ( ptrb(:,:,:,jn) - ztrtrdb(:,:,:) ) * zs2rdt 192 192 ztrtrdn(:,:,:) = ( ptrn(:,:,:,jn) - ztrtrdn(:,:,:) ) * zs2rdt -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/TRP/trczdf.F90
r3116 r3124 74 74 IF( neuler == 0 .AND. kt == nittrc000 ) THEN ! at nittrc000 75 75 r2dt(:) = rdttrc(:) ! = rdttrc (restarting with Euler time stepping) 76 ELSEIF( kt <= nittrc000 + nn_dttrc ) THEN ! at nittrc000 or nittrc000+nn_dttrc76 ELSEIF( kt <= nittrc000 + 1 ) THEN ! at nittrc000 or nittrc000+1 77 77 r2dt(:) = 2. * rdttrc(:) ! = 2 rdttrc (leapfrog) 78 78 ENDIF … … 88 88 SELECT CASE ( nzdf ) ! compute lateral mixing trend and add it to the general trend 89 89 CASE ( -1 ) ! esopa: test all possibility with control print 90 CALL tra_zdf_exp( kt, nittrc000, 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra )90 CALL tra_zdf_exp( kt, nittrc000, 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra ) 91 91 WRITE(charout, FMT="('zdf1 ')") ; CALL prt_ctl_trc_info(charout) 92 92 CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 93 CALL tra_zdf_imp( kt, nittrc000, 'TRC', r2dt, trb, tra, jptra )93 CALL tra_zdf_imp( kt, nittrc000, 'TRC', r2dt, trb, tra, jptra ) 94 94 WRITE(charout, FMT="('zdf2 ')") ; CALL prt_ctl_trc_info(charout) 95 95 CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 96 CASE ( 0 ) ; CALL tra_zdf_exp( kt, nittrc000, 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra ) ! explicit scheme97 CASE ( 1 ) ; CALL tra_zdf_imp( kt, nittrc000, 'TRC', r2dt, trb, tra, jptra ) ! implicit scheme96 CASE ( 0 ) ; CALL tra_zdf_exp( kt, nittrc000, 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra ) ! explicit scheme 97 CASE ( 1 ) ; CALL tra_zdf_imp( kt, nittrc000, 'TRC', r2dt, trb, tra, jptra ) ! implicit scheme 98 98 99 99 END SELECT -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/trcdia.F90
r3116 r3124 71 71 !!--------------------------------------------------------------------- 72 72 ! 73 IF( kt == nit 000 ) THEN73 IF( kt == nittrc000 ) THEN 74 74 ALLOCATE( ndext50(jpij*jpk), ndext51(jpij), STAT=ierr ) 75 75 IF( ierr > 0 ) THEN … … 145 145 itmod = kt - nittrc000 + 1 146 146 it = kt 147 iiter = ( nit 000 - 1 ) / nn_dttrc147 iiter = ( nittrc000 - 1 ) / nn_dttrc 148 148 149 149 ! Define NETCDF files and fields at beginning of first time step … … 292 292 itmod = kt - nittrc000 + 1 293 293 it = kt 294 iiter = ( nit 000 - 1 ) / nn_dttrc294 iiter = ( nittrc000 - 1 ) / nn_dttrc 295 295 296 296 ! 1. Define NETCDF files and fields at beginning of first time step … … 433 433 itmod = kt - nittrc000 + 1 434 434 it = kt 435 iiter = ( nit 000 - 1 ) / nn_dttrc435 iiter = ( nittrc000 - 1 ) / nn_dttrc 436 436 437 437 ! Define NETCDF files and fields at beginning of first time step -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/trcdta.F90
r2977 r3124 177 177 IF( ln_sco ) THEN !== s- or mixed s-zps-coordinate ==! 178 178 ! 179 IF( kt == nit 000 .AND. lwp )THEN179 IF( kt == nittrc000 .AND. lwp )THEN 180 180 WRITE(numout,*) 181 181 WRITE(numout,*) 'trc_dta: interpolates passive tracer data onto the s- or mixed s-z-coordinate mesh' … … 232 232 ENDDO 233 233 ! 234 IF( lwp .AND. kt == nit 000 ) THEN234 IF( lwp .AND. kt == nittrc000 ) THEN 235 235 DO jn = 1, ntra 236 236 clndta = TRIM( sf_trcdta(jn)%clvar ) -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/trcini.F90
r3116 r3124 121 121 ENDIF 122 122 ! 123 CALL trc_dta( nit 000, ztrcdta ) ! read tracer data at nit000123 CALL trc_dta( nittrc000, ztrcdta ) ! read tracer data at nittrc000 124 124 ! 125 125 DO jn = 1, jptra -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r3116 r3124 66 66 67 67 IF( MOD( kt - 1, nstock ) == 0 ) THEN 68 ! we use kt - 1 and not kt - nit 000 to keep the same periodicity from the beginning of the experiment68 ! we use kt - 1 and not kt - nittrc000 to keep the same periodicity from the beginning of the experiment 69 69 nitrst = kt + nstock - 1 ! define the next value of nitrst for restart writing 70 70 IF( nitrst > nitend ) nitrst = nitend ! make sure we write a restart at the end of the run … … 200 200 !! In both those options, the exact duration of the experiment 201 201 !! since the beginning (cumulated duration of all previous restart runs) 202 !! is not stored in the restart and is assumed to be (nit 000-1)*rdt.202 !! is not stored in the restart and is assumed to be (nittrc000-1)*rdt. 203 203 !! This is valid is the time step has remained constant. 204 204 !!
Note: See TracChangeset
for help on using the changeset viewer.