New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 1011 for trunk/NEMO/TOP_SRC/trcrst.F90 – NEMO

Ignore:
Timestamp:
2008-05-30T11:57:12+02:00 (16 years ago)
Author:
cetlod
Message:

re-organization of TOP initialization and outputs phases see ticket 171

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/trcrst.F90

    r945 r1011  
    3232   INTEGER, PUBLIC ::   numrtr, numrtw   !: logical unit for trc restart (read and write) 
    3333 
     34#if defined key_pisces 
    3435   REAL(wp) ::  & 
    3536     alkmean = 2426. ,  & ! mean value of alkalinity ( Glodap ; for Goyet 2391. ) 
     
    3738     no3mean = 30.90 ,  & ! mean value of nitrate 
    3839     siomean = 91.51      ! mean value of silicate 
     40#endif 
    3941 
    4042   !! * Substitutions 
     
    158160      CALL iom_get( numrtr, jpdom_local, 'SEDB'//ctrcnm(jpdet), sedpocb(:,:) )  
    159161      CALL iom_get( numrtr, jpdom_local, 'SEDN'//ctrcnm(jpdet), sedpocn(:,:) )  
    160 # elif defined key_pisces 
     162#endif 
     163#if defined key_pisces 
    161164      CALL iom_get( numrtr, jpdom_local, 'Silicalim', xksi(:,:) )  
    162165      CALL iom_get( numrtr, jpdom_local, 'Silicamax', xksimax(:,:) ) 
    163 # elif defined key_cfc 
    164       DO jn = 1, jptra 
     166#endif 
     167#if defined key_cfc 
     168      DO jn = jp_cfc0, jp_cfc1 
    165169         CALL iom_get( numrtr, jpdom_local, 'qint'//ctrcnm(jn), qint(:,:,jn) )  
    166170         CALL iom_get( numrtr, jpdom_local, 'qtr'//ctrcnm(jn) , qtr( :,:,jn) )  
     
    325329         CALL iom_rstput( kt, nitrst, numrtw, 'SEDB'//ctrcnm(jpdet), sedpocb(:,:) ) 
    326330         CALL iom_rstput( kt, nitrst, numrtw, 'SEDN'//ctrcnm(jpdet), sedpocn(:,:) ) 
    327 #elif defined key_pisces 
     331#endif 
     332#if defined key_pisces 
    328333         CALL iom_rstput( kt, nitrst, numrtw, 'Silicalim', xksi(:,:) ) 
    329334         CALL iom_rstput( kt, nitrst, numrtw, 'Silicamax', xksimax(:,:) ) 
    330  
    331 #elif defined key_cfc 
    332          DO jn = 1, jptra 
     335#endif 
     336#if defined key_cfc 
     337         DO jn = jp_cfc0, jp_cfc1 
    333338            CALL iom_rstput( kt, nitrst, numrtw, 'qint'//ctrcnm(jn), qint(:,:,jn) ) 
    334339            CALL iom_rstput( kt, nitrst, numrtw, 'qtr'//ctrcnm(jn) , qtr (:,:,jn) ) 
Note: See TracChangeset for help on using the changeset viewer.