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 8607 for branches/NERC – NEMO

Changeset 8607 for branches/NERC


Ignore:
Timestamp:
2017-10-10T12:43:01+02:00 (5 years ago)
Author:
jpalmier
Message:

JPALM -- enable to scale Chl sent to the UM and choose between surf or MLD chl

Location:
branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_diag_slice.F90

    r8442 r8607  
    4848                                   jdms, ocal_ccd, xpar, xze,              & 
    4949                                   zb_co2_flx, zb_dms_srf,                 & 
    50                                    zn_co2_flx, zn_dms_srf, zn_chl_srf 
     50                                   zn_co2_flx, zn_dms_srf, zn_chl_srf,     & 
     51                                   scl_chl, chl_out 
    5152      USE trc,               ONLY: med_diag 
    5253 
     
    7071         !! no need to output, just pass to cpl var 
    7172         IF (lk_oasis) THEN 
    72             zn_chl_srf(:,:) = (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * 1.0E-6  !! surf Chl in Kg-chl/m3 as needed for cpl 
     73            IF (chl_out.eq.1) THEN 
     74               !! export and scale surface chl 
     75               zn_chl_srf(:,:) = (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * scl_chl * 1.0E-6   
     76                                 !! surf Chl in Kg-chl/m3 as needed for cpl 
     77            ELSEIF (chl_out.eq.2) THEN 
     78               !! export and scale mld chl 
     79               zn_chl_srf(:,:) = fchl_ml * scl_chl * 1.0E-6   
     80                                 !! mld Chl in Kg-chl/m3 as needed for cpl 
     81            ENDIF 
    7382            chloro_out_cpl(:,:) = zn_chl_srf(:,:)        !! Coupling Chl 
    7483         END IF 
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90

    r8521 r8607  
    5151                                   zb_sed_n, zb_sed_si,                 & 
    5252                                   zn_sed_c, zn_sed_ca, zn_sed_fe,      & 
    53                                    zn_sed_n, zn_sed_si 
     53                                   zn_sed_n, zn_sed_si, zn_chl_srf 
    5454      USE trc,               ONLY: med_diag, nittrc000  
    5555      USE trcnam_trp,        ONLY: ln_trcadv_cen2, ln_trcadv_tvd 
     
    241241#   endif 
    242242         IF ( med_diag%INVTN%dgsave ) THEN 
    243             CALL iom_put( "INVTN"  , ftot_n ) 
     243            !! HACK - Coupled Chl sent to UM     
     244            CALL iom_put( "INVTN"  , zn_chl_srf ) 
     245            !! CALL iom_put( "INVTN"  , ftot_n ) 
    244246         ENDIF 
    245247         IF ( med_diag%INVTSI%dgsave ) THEN 
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90

    r8132 r8607  
    177177   INTEGER  ::  jdms_model   !: choice of DMS model passed to atmosphere 
    178178!!                              1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL 
    179 !! 
     179!! FOR UKESM    
     180   REAL(wp) ::  scl_chl      !: scaling factor for tuned Chl passed to the UM  
     181   INTEGER  ::  chl_out      !: select Chl field exported and scaled for UM: 
     182                             !: 1- Surface Chl ; 2- MLD Chl 
    180183!! 
    181184   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   remdmp   !: depth dependent damping coefficient of passive tracers  
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90

    r8521 r8607  
    8585      &  jriver_n,jriver_si,jriver_c,jriver_alk,jriver_dep,   & 
    8686      &  xsdiss,                                              & 
    87       &  sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model 
     87      &  sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model,   & 
     88      &  scl_chl, chl_out 
    8889#if defined key_roam 
    8990      NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit,        & 
     
    246247      jdms_input  = 0 
    247248      jdms_model  = 0 
     249      scl_chl     = 0.5 
     250      chl_out     = 1 
    248251             
    249252      !REWIND(numnatm) 
     
    400403!!                      1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL, 5 = ANDM 
    401404!! 
     405!! UKESM1 - exported Chl to UM 
     406!!       scl_chl     : scaling factor to tune the chl field sent to the UM 
     407!!       chl_out     : select the chl field to send at the UM: 
     408!!                     1- Surf Chl ; 2- MLD Chl  
     409 
    402410      IF(lwp) THEN 
    403411!! 
     
    909917!! 
    910918!! UKESM1 - new diagnostics  !! Jpalm; AXY (08/07/15) 
    911          WRITE(numout,*) '=== UKESM1-related parameters' 
     919         WRITE(numout,*) '=== UKESM1-related parameters ===' 
     920         WRITE(numout,*) ' ---- --- ---' 
     921         WRITE(numout,*) '=== UKESM1 --  coupled DMS to the UM' 
    912922         WRITE(numout,*)     & 
    913923         &   ' include DMS diagnostic?,                                   jdms        = ', jdms 
     
    935945            &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Anderson et al. (2001; modified)' 
    936946         endif 
     947         WRITE(numout,*) '=== UKESM1 --  coupled Chl to the UM' 
     948         WRITE(numout,*)        & 
     949            &   ' Scaling factor to export tuned Chl to the UM               scl_chl  = ', scl_chl 
     950         IF (chl_out .eq. 1) THEN 
     951            WRITE(numout,*)        & 
     952            &   ' Chl field to be scaled and sent to the UM:                 chl_out  = Surface Chl field ' 
     953         ELSEIF (chl_out .eq. 2) THEN 
     954            WRITE(numout,*)        & 
     955            &   ' Chl field to be scaled and sent to the UM:                 chl_out  = MLD Chl field ' 
     956         ENDIF 
    937957!! 
    938958      ENDIF 
  • branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

    r8427 r8607  
    282282      ELSE 
    283283         IF(lwp) WRITE(numout,*) 'Chl surf concentration - setting to zero ...' 
    284          zn_chl_srf(:,:)  = (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * 1.E-6 
     284         zn_chl_srf(:,:)  = (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * scl_chl * 1.E-6 
    285285      ENDIF 
    286286      IF (lk_oasis) THEN 
Note: See TracChangeset for help on using the changeset viewer.