Changeset 13107


Ignore:
Timestamp:
2020-06-12T20:34:30+02:00 (6 months ago)
Author:
lovato
Message:

add trcopt.F90 containing generalized module for vertical light propagation (ticket #2489)

Location:
NEMO/branches/2020/dev_r12985_TOP-04_IMMERSE_BGC_interface/src/TOP
Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12985_TOP-04_IMMERSE_BGC_interface/src/TOP/PISCES/sms_pisces.F90

    r12377 r13107  
    7070 
    7171   !!*  Biological fluxes for light : variables shared by pisces & lobster 
    72    INTEGER , ALLOCATABLE, SAVE, DIMENSION(:,:)   ::  neln  !: number of T-levels + 1 in the euphotic layer 
    73    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::  heup  !: euphotic layer depth 
    74    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::  etot  !: par (photosynthetic available radiation) 
    75    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::  etot_ndcy      !: PAR over 24h in case of diurnal cycle 
    7672   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::  enano, ediat   !: PAR for phyto, nano and diat  
    7773   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::  enanom, ediatm !: PAR for phyto, nano and diat  
     
    7975   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::  epicom         !: PAR for pico 
    8076   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::  emoy           !: averaged PAR in the mixed layer 
    81    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::  heup_01 !: Absolute euphotic layer depth 
    8277   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::  xksi  !:  LOBSTER : zooplakton closure 
    8378 
     
    139134      ierr(:) = 0 
    140135      !*  Biological fluxes for light : shared variables for pisces & lobster 
    141       ALLOCATE( etot(jpi,jpj,jpk), neln(jpi,jpj), heup(jpi,jpj),    & 
    142         &       heup_01(jpi,jpj) , xksi(jpi,jpj)               ,  STAT=ierr(1) ) 
     136      ALLOCATE( xksi(jpi,jpj)  ,                                  STAT=ierr(1) )  
    143137      ! 
    144    
    145138      IF( ln_p4z .OR. ln_p5z ) THEN 
    146139         !*  Biological fluxes for light  
    147140         ALLOCATE(  enano(jpi,jpj,jpk)    , ediat(jpi,jpj,jpk) ,   & 
    148141           &        enanom(jpi,jpj,jpk)   , ediatm(jpi,jpj,jpk),   & 
    149            &        etot_ndcy(jpi,jpj,jpk), emoy(jpi,jpj,jpk)  ,  STAT=ierr(2) )  
     142           &        emoy(jpi,jpj,jpk)     ,                       STAT=ierr(2) )  
    150143 
    151144         !*  Biological fluxes for primary production 
  • NEMO/branches/2020/dev_r12985_TOP-04_IMMERSE_BGC_interface/src/TOP/trc.F90

    r12489 r13107  
    8282   CHARACTER(len=2), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   cn_trc_o         !: choice of ocean tracer cc 
    8383 
     84   !! Information for the optics module 
     85   !! --------------------------------- 
     86   INTEGER , ALLOCATABLE, SAVE, DIMENSION(:,:)   ::  neln       !: number of T-levels + 1 in the euphotic layer 
     87   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::  heup       !: euphotic layer depth 
     88   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::  heup_01    !: Absolute euphotic layer depth 
     89   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::  etot       !: par (photosynthetic available radiation) 
     90   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::  etot_ndcy  !: PAR over 24h in case of diurnal cycle 
     91 
    8492 
    8593   !! information for outputs 
     
    154162         &      gtrui(jpi,jpj,jptra)  , gtrvi(jpi,jpj,jptra)                          ,       & 
    155163         &      trc_ice_ratio(jptra)  , trc_ice_prescr(jptra) , cn_trc_o(jptra)       ,       & 
     164         &      neln(jpi,jpj)         , heup(jpi,jpj)         , heup_01(jpi,jpj)      ,       & 
     165         &      etot(jpi,jpj,jpk)     , etot_ndcy(jpi,jpj,jpk),                               & 
    156166         &      sbc_trc_b(jpi,jpj,jptra), sbc_trc(jpi,jpj,jptra)                      ,       &   
    157167         &      cvol(jpi,jpj,jpk)     , trai(jptra)           , qsr_mean(jpi,jpj)     ,       & 
Note: See TracChangeset for help on using the changeset viewer.