Changeset 1005


Ignore:
Timestamp:
2008-05-30T11:48:49+02:00 (13 years ago)
Author:
cetlod
Message:

Update LOBSTER modules to take into account the re-organization of TOP initialization phase, see ticket 169

Location:
trunk/NEMO/TOP_SRC/LOBSTER
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/LOBSTER/trcini_lobster.F90

    r934 r1005  
    1919   USE sms             ! Source Minus Sink variables 
    2020   USE oce_trc         ! ocean variables 
     21   USE trc 
    2122 
    2223   IMPLICIT NONE 
     
    142143      END DO 
    143144      ! 
     145 
     146      ! Initialization of tracer concentration in case of  no restart  
     147      !------------------------------------------------------------- 
     148      IF( .NOT. lrsttr ) THEN     
     149 
     150# if defined key_eel_r6 || defined key_eel_r2 
     151         ! LOBSTER initialisation for EEL 
     152         ! ---------------------- 
     153         ! here: analytical initialisation used in Levy et al. (2001) 
     154          
     155         DO jk = 1, 7 
     156            trn(:,:,jk,jpdet) = 0.016 * tmask(:,:,jk) 
     157            trn(:,:,jk,jpzoo) = 0.018 * tmask(:,:,jk) 
     158            trn(:,:,jk,jpphy) = 0.036 * tmask(:,:,jk)  
     159            trn(:,:,jk,jpno3) = 1.e-5 * tmask(:,:,jk) 
     160            trn(:,:,jk,jpnh4) = 5.e-4 * tmask(:,:,jk) 
     161            trn(:,:,jk,jpdom) = 0.017 * tmask(:,:,jk) 
     162         END DO 
     163          
     164         trn(:,:, 8,jpdet) = 0.020   * tmask(:,:, 8) 
     165         trn(:,:, 8,jpzoo) = 0.027   * tmask(:,:, 8) 
     166         trn(:,:, 8,jpphy) = 0.041   * tmask(:,:, 8) 
     167         trn(:,:, 8,jpno3) = 0.00022 * tmask(:,:, 8) 
     168         trn(:,:, 8,jpnh4) = 0.0033  * tmask(:,:, 8) 
     169         trn(:,:, 8,jpdom) = 0.021   * tmask(:,:, 8) 
     170          
     171         trn(:,:, 9,jpdet) = 0.0556  * tmask(:,:, 9) 
     172         trn(:,:, 9,jpzoo) = 0.123   * tmask(:,:, 9) 
     173         trn(:,:, 9,jpphy) = 0.122   * tmask(:,:, 9) 
     174         trn(:,:, 9,jpno3) = 0.028   * tmask(:,:, 9) 
     175         trn(:,:, 9,jpnh4) = 0.024   * tmask(:,:, 9) 
     176         trn(:,:, 9,jpdom) = 0.06    * tmask(:,:, 9) 
     177          
     178         trn(:,:,10,jpdet) = 0.025   * tmask(:,:,10) 
     179         trn(:,:,10,jpzoo) = 0.016   * tmask(:,:,10) 
     180         trn(:,:,10,jpphy) = 0.029   * tmask(:,:,10) 
     181         trn(:,:,10,jpno3) = 2.462   * tmask(:,:,10) 
     182         trn(:,:,10,jpnh4) = 0.04    * tmask(:,:,10) 
     183         trn(:,:,10,jpdom) = 0.022   * tmask(:,:,10) 
     184          
     185         trn(:,:,11,jpdet) = 0.0057  * tmask(:,:,11) 
     186         trn(:,:,11,jpzoo) = 0.0005  * tmask(:,:,11) 
     187         trn(:,:,11,jpphy) = 0.0006  * tmask(:,:,11) 
     188         trn(:,:,11,jpno3) = 3.336   * tmask(:,:,11) 
     189         trn(:,:,11,jpnh4) = 0.005   * tmask(:,:,11) 
     190         trn(:,:,11,jpdom) = 0.004   * tmask(:,:,11) 
     191          
     192         trn(:,:,12,jpdet) = 0.002   * tmask(:,:,12) 
     193         trn(:,:,12,jpzoo) = 1.e-6   * tmask(:,:,12) 
     194         trn(:,:,12,jpphy) = 5.e-6   * tmask(:,:,12) 
     195         trn(:,:,12,jpno3) = 4.24    * tmask(:,:,12) 
     196         trn(:,:,12,jpnh4) = 0.001   * tmask(:,:,12) 
     197         trn(:,:,12,jpdom) = 3.e-5   * tmask(:,:,12) 
     198          
     199         DO jk=13,jpk 
     200            trn(:,:,jk,jpdet) = 0.e0 
     201            trn(:,:,jk,jpzoo) = 0.e0 
     202            trn(:,:,jk,jpphy) = 0.e0 
     203            trn(:,:,jk,jpnh4) = 0.e0 
     204            trn(:,:,jk,jpdom) = 0.e0 
     205         END DO 
     206          
     207         trn(:,:,13,jpno3) = 5.31  * tmask(:,:,13) 
     208         trn(:,:,14,jpno3) = 6.73  * tmask(:,:,14) 
     209         trn(:,:,15,jpno3) = 8.32  * tmask(:,:,15) 
     210         trn(:,:,16,jpno3) = 10.13 * tmask(:,:,16) 
     211         trn(:,:,17,jpno3) = 11.95 * tmask(:,:,17) 
     212         trn(:,:,18,jpno3) = 13.57 * tmask(:,:,18) 
     213         trn(:,:,19,jpno3) = 15.08 * tmask(:,:,19) 
     214         trn(:,:,20,jpno3) = 16.41 * tmask(:,:,20) 
     215         trn(:,:,21,jpno3) = 17.47 * tmask(:,:,21) 
     216         trn(:,:,22,jpno3) = 18.29 * tmask(:,:,22) 
     217         trn(:,:,23,jpno3) = 18.88 * tmask(:,:,23) 
     218         trn(:,:,24,jpno3) = 19.30 * tmask(:,:,24) 
     219         trn(:,:,25,jpno3) = 19.68 * tmask(:,:,25) 
     220         trn(:,:,26,jpno3) = 19.91 * tmask(:,:,26) 
     221         trn(:,:,27,jpno3) = 19.99 * tmask(:,:,27) 
     222         trn(:,:,28,jpno3) = 20.01 * tmask(:,:,28) 
     223         trn(:,:,29,jpno3) = 20.01 * tmask(:,:,29) 
     224         trn(:,:,30,jpno3) = 20.01 * tmask(:,:,30) 
     225 
     226# elif defined key_gyre 
     227         ! LOBSTER initialisation for GYRE 
     228         ! ---------------------- 
     229         ! here:  init NO3=f(density) by asklod AS Kremeur 2005-07 
     230         trn(:,:,:,jpdet) = 0.1 * tmask(:,:,:) 
     231         trn(:,:,:,jpzoo) = 0.1 * tmask(:,:,:) 
     232         trn(:,:,:,jpnh4) = 0.1 * tmask(:,:,:) 
     233         trn(:,:,:,jpphy) = 0.1 * tmask(:,:,:) 
     234         trn(:,:,:,jpdom) = 1.0 * tmask(:,:,:) 
     235         DO jk = 1, jpk 
     236            DO jj = 1, jpj 
     237               DO ji = 1, jpi 
     238                  IF( rhd(ji,jj,jk) <= 24.5e-3 ) THEN 
     239                     trn(ji,jj,jk,jpno3) = 2. * tmask(ji,jj,jk) 
     240                  ELSE 
     241                     trn(ji,jj,jk,jpno3) = ( 15.55 * ( rhd(ji,jj,jk) * 1000. ) - 380.11 ) * tmask(ji,jj,jk) 
     242                  ENDIF 
     243               END DO 
     244            END DO 
     245         END DO 
     246#endif 
     247 
     248      ENDIF 
     249 
     250      !  initialize the POC in sediments 
     251      sedpocb(:,:) = 0.e0 
     252      sedpocn(:,:) = 0.e0 
     253      sedpoca(:,:) = 0.e0 
     254 
     255 
     256      IF(lwp) WRITE(numout,*) 'Initialization of LOBSTER tracers done' 
     257      IF(lwp) WRITE(numout,*) ' ' 
     258 
     259 
    144260   END SUBROUTINE trc_ini_lobster 
    145261 
  • trunk/NEMO/TOP_SRC/LOBSTER/trclsm_lobster.F90

    r934 r1005  
    6262      !                               ! Open the namelist file 
    6363      !                               ! ---------------------- 
    64       clname ='namelist.trc.sms' 
     64      clname ='namelist.sms.lobster' 
    6565      CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL',   & 
    6666         &           1, numout, .FALSE., 1 ) 
  • trunk/NEMO/TOP_SRC/LOBSTER/trcsed.F90

    r934 r1005  
    7979      ! tracer flux at w-point: we use -vsed (downward flux)  with simplification : no e1*e2 
    8080 
    81       DO jk = 2, jpkm1 
     81      DO jk = 2, jpk 
    8282         zwork(:,:,jk) = -vsed * trn(:,:,jk-1,jpdet) 
    8383      END DO 
Note: See TracChangeset for help on using the changeset viewer.