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 3318 for branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/LOBSTER – NEMO

Ignore:
Timestamp:
2012-02-25T16:50:01+01:00 (12 years ago)
Author:
gm
Message:

Ediag branche: #927 split TRA/DYN trd computation

Location:
branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/LOBSTER
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcbio.F90

    r3294 r3318  
    22   !!====================================================================== 
    33   !!                         ***  MODULE trcbio  *** 
    4    !! TOP :   LOBSTER 
     4   !! TOP :   LOBSTER bio model 
    55   !!====================================================================== 
    6    !! History :    -   !  1999-07  (M. Levy) Original code 
    7    !!              -   !  2000-12  (E. Kestenare) assign a parameter to name individual tracers 
    8    !!              -   !  2001-03  (M. Levy)  LNO3 + dia2d  
    9    !!             2.0  !  2007-12  (C. Deltel, G. Madec)  F90 
     6   !! History :  OPA  !  1999-07  (M. Levy) Original code 
     7   !!             -   !  2000-12  (E. Kestenare) assign a parameter to name individual tracers 
     8   !!             -   !  2001-03  (M. Levy)  LNO3 + dia2d  
     9   !!   NEMO     2.0  !  2007-12  (C. Deltel, G. Madec)  F90 
    1010   !!---------------------------------------------------------------------- 
    1111#if defined key_lobster 
     
    1313   !!   'key_lobster'                                     LOBSTER bio-model 
    1414   !!---------------------------------------------------------------------- 
    15    !!   trc_bio         
    16    !!---------------------------------------------------------------------- 
    17    USE oce_trc         ! 
    18    USE trc             !  
    19    USE sms_lobster     !  
    20    USE lbclnk          !  
    21    USE prtctl_trc      ! Print control for debbuging 
    22    USE trdmod_oce 
    23    USE trdmod_trc 
     15   !!   trc_bio        
     16   !!---------------------------------------------------------------------- 
     17   USE oce_trc        ! 
     18   USE trc            !  
     19   USE sms_lobster    !  
     20   USE lbclnk         !  
     21   USE prtctl_trc     ! Print control for debbuging 
     22   USE trd_oce        ! trends: ocean variables 
     23   USE trdmod_trc     ! trends: passive tracer variables 
    2424   USE iom 
    2525    
     
    2727   PRIVATE 
    2828 
    29    PUBLIC   trc_bio    ! called in ??? 
     29   PUBLIC   trc_bio   ! called in ??? 
    3030 
    3131   !!* Substitution 
     
    3636   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3737   !!---------------------------------------------------------------------- 
    38  
    3938CONTAINS 
    4039 
     
    6059      !!              for passive tracers are saved for futher diagnostics. 
    6160      !!--------------------------------------------------------------------- 
    62       !! 
    6361      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    64       !! 
     62      ! 
    6563      INTEGER  ::   ji, jj, jk, jl 
    6664      REAL(wp) ::   zdet, zzoo, zphy, zno3, znh4, zdom      ! now concentrations 
  • branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcexp.F90

    r3294 r3318  
    44   !! TOP :   PISCES Compute loss of organic matter in the sediments 
    55   !!====================================================================== 
    6    !! History :    -   !  1999    (O. Aumont, C. Le Quere)  original code 
     6   !! History :   OPA  !  1999    (O. Aumont, C. Le Quere)  original code 
    77   !!              -   !  2001-05 (O. Aumont, E. Kestenare) add sediment computations 
    8    !!             1.0  !  2005-06 (A.-S. Kremeur) new temporal integration for sedpoc 
     8   !!   NEMO      1.0  !  2005-06 (A.-S. Kremeur) new temporal integration for sedpoc 
    99   !!             2.0  !  2007-12  (C. Deltel, G. Madec)  F90 
    1010   !!---------------------------------------------------------------------- 
     
    2121   USE trcnam_trp 
    2222   USE prtctl_trc      ! Print control for debbuging 
    23    USE trdmod_oce 
    24    USE trdmod_trc 
     23   USE trd_oce         ! trends: ocean variables 
     24   USE trdmod_trc      ! trends: passive tracer variables 
    2525   USE iom 
    2626 
     
    3737   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3838   !!---------------------------------------------------------------------- 
    39  
    4039CONTAINS 
    4140 
     
    5352      !!              COLUMN BELOW THE SURFACE LAYER. 
    5453      !!--------------------------------------------------------------------- 
    55       !! 
    5654      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    57       !! 
     55      ! 
    5856      INTEGER  ::   ji, jj, jk, jl, ikt, ierr 
    5957      REAL(wp) ::   zgeolpoc, zfact, zwork, ze3t, zsedpocd 
  • branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcsed.F90

    r3294 r3318  
    1818   USE sms_lobster 
    1919   USE lbclnk 
    20    USE trdmod_oce 
    21    USE trdmod_trc 
     20   USE trd_oce         ! trends: ocean variables 
     21   USE trdmod_trc      ! trends: passive tracer variables 
    2222   USE iom 
    2323   USE prtctl_trc      ! Print control for debbuging 
     
    3535   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3636   !!---------------------------------------------------------------------- 
    37  
    3837CONTAINS 
    3938 
     
    5655      !!              trend of passive tracers is saved for futher diagnostics. 
    5756      !!--------------------------------------------------------------------- 
    58       !! 
    5957      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    60       !! 
     58      ! 
    6159      INTEGER  ::   ji, jj, jk, jl, ierr 
    6260      CHARACTER (len=25) :: charout 
    63       REAL(wp), POINTER, DIMENSION(:,:  ) :: zw2d 
    64       REAL(wp), POINTER, DIMENSION(:,:,:) :: zwork, ztra, ztrbio 
     61      REAL(wp), POINTER, DIMENSION(:,:  ) ::   zw2d 
     62      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zwork, ztra, ztrbio 
    6563      !!--------------------------------------------------------------------- 
    6664      ! 
     
    7674      CALL wrk_alloc( jpi, jpj, jpk, zwork, ztra ) 
    7775 
    78       IF( ln_diatrc )  THEN 
    79          CALL wrk_alloc( jpi, jpj, zw2d ) 
    80       ENDIF 
     76      IF( ln_diatrc )   CALL wrk_alloc( jpi, jpj, zw2d ) 
    8177 
    8278      IF( l_trdtrc ) THEN 
     
    109105      IF( ln_diatrc ) THEN  
    110106         DO jk = 1, jpkm1 
    111             DO jj = 1, jpj 
    112                DO ji = 1, jpi 
    113                   zw2d(ji,jj) = zw2d(ji,jj) + ztra(ji,jj,jk) * fse3t(ji,jj,jk) * 86400. 
    114                END DO 
    115             END DO 
     107            zw2d(:,:) = zw2d(:,:) + ztra(:,:,jk) * fse3t(:,:,jk) * 86400._wp 
    116108         END DO 
    117          IF( lk_iomput )  THEN 
    118            CALL iom_put( "TDETSED", zw2d ) 
    119          ELSE 
    120            trc2d(:,:,jp_lob0_2d + 7) = zw2d(:,:) 
     109         IF( lk_iomput ) THEN   ;   CALL iom_put( "TDETSED", zw2d ) 
     110         ELSE                   ;   trc2d(:,:,jp_lob0_2d + 7) = zw2d(:,:) 
    121111         ENDIF 
    122112         CALL wrk_dealloc( jpi, jpj, zw2d ) 
    123113      ENDIF 
    124114      ! 
    125       IF( ln_diabio )  trbio(:,:,:,jp_lob0_trd + 7) = ztra(:,:,:) 
     115      IF( ln_diabio )   trbio(:,:,:,jp_lob0_trd + 7) = ztra(:,:,:) 
    126116      CALL wrk_dealloc( jpi, jpj, jpk, zwork, ztra ) 
    127117      ! 
     
    133123      ENDIF 
    134124 
    135       IF(ln_ctl)   THEN  ! print mean trends (used for debugging) 
     125      IF(ln_ctl) THEN      ! print mean trends (used for debugging) 
    136126         WRITE(charout, FMT="('sed')") 
    137127         CALL prt_ctl_trc_info(charout) 
     
    139129      ENDIF 
    140130      ! 
    141       IF( nn_timing == 1 )  CALL timing_stop('trc_sed') 
     131      IF( nn_timing == 1 )   CALL timing_stop('trc_sed') 
    142132      ! 
    143133   END SUBROUTINE trc_sed 
  • branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcsms_lobster.F90

    r3294 r3318  
    44   !! TOP :   Time loop of LOBSTER model 
    55   !!====================================================================== 
    6    !! History :   1.0  !            M. Levy 
     6   !! History :   1.0  !  2004-02  (M. Levy)  original code 
    77   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  revised architecture 
    88   !!---------------------------------------------------------------------- 
     
    1313   !!   trcsms_lobster        :  Time loop of passive tracers sms 
    1414   !!---------------------------------------------------------------------- 
    15    USE oce_trc          ! 
     15   USE oce_trc       ! 
    1616   USE trc 
    1717   USE trcbio 
     
    1919   USE trcsed 
    2020   USE trcexp 
    21    USE trdmod_oce 
    22    USE trdmod_trc_oce 
     21   USE trd_oce       ! trends: ocean variables 
     22   USE trdmod_trc    ! trends: passive tracer variables 
    2323   USE trdmod_trc 
    2424   USE trdmld_trc 
     
    6060      IF( l_trdtrc ) THEN 
    6161         DO jn = jp_lob0, jp_lob1 
    62            CALL trd_mod_trc( tra(:,:,:,jn), jn, jptra_trd_sms, kt )   ! save trends 
     62           CALL trd_mod_trc( tra(:,:,:,jn), jn, jptra_sms, kt )   ! save trends 
    6363         END DO 
    6464      END IF 
Note: See TracChangeset for help on using the changeset viewer.