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 8436 for branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90 – NEMO

Ignore:
Timestamp:
2017-08-14T15:22:09+02:00 (7 years ago)
Author:
dford
Message:

Implement initial version of surface chlorophyll assimilation for MEDUSA.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90

    r8400 r8436  
    4949#if defined key_lim3 
    5050   USE ice 
     51#endif 
     52#if defined key_hadocc 
     53   USE trc, ONLY: trn, & 
     54      &     pgrow_avg, & 
     55      &     ploss_avg, & 
     56      &     phyt_avg,  & 
     57      &     mld_max,   & 
     58      &     HADOCC_CHL 
     59   USE had_bgc_const, ONLY: cchl_p 
     60   USE par_hadocc 
     61#elif defined key_medusa && defined key_foam_medusa 
     62   USE trc, ONLY: trn 
     63   USE sms_medusa, ONLY: pgrow_avg, & 
     64      &                  ploss_avg, & 
     65      &                  phyt_avg,  & 
     66      &                  mld_max 
     67   USE par_medusa 
    5168#endif 
    5269   IMPLICIT NONE 
    … …  
    121138!            CALL iom_rstput( kt, nitbkg_r, inum, 'gcx'    , gcx               ) 
    122139            CALL iom_rstput( kt, nitbkg_r, inum, 'avt'    , avt               ) 
     140#if defined key_hadocc 
     141            CALL iom_rstput( kt, nitbkg_r, inum, 'pgrow_avg'     , pgrow_avg             ) 
     142            CALL iom_rstput( kt, nitbkg_r, inum, 'ploss_avg'     , ploss_avg             ) 
     143            CALL iom_rstput( kt, nitbkg_r, inum, 'phyt_avg'      , phyt_avg              ) 
     144            CALL iom_rstput( kt, nitbkg_r, inum, 'mld_max'       , mld_max               ) 
     145            CALL iom_rstput( kt, nitbkg_r, inum, 'nutrients'     , trn(:,:,:,jp_had_nut) ) 
     146            CALL iom_rstput( kt, nitbkg_r, inum, 'phytoplankton' , trn(:,:,:,jp_had_phy) ) 
     147            CALL iom_rstput( kt, nitbkg_r, inum, 'zooplankton'   , trn(:,:,:,jp_had_zoo) ) 
     148            CALL iom_rstput( kt, nitbkg_r, inum, 'detritus'      , trn(:,:,:,jp_had_pdn) ) 
     149            CALL iom_rstput( kt, nitbkg_r, inum, 'dic'           , trn(:,:,:,jp_had_dic) ) 
     150            CALL iom_rstput( kt, nitbkg_r, inum, 'alkalinity'    , trn(:,:,:,jp_had_alk) ) 
     151            CALL iom_rstput( kt, nitbkg_r, inum, 'chlorophyll'   , HADOCC_CHL(:,:,1)     ) 
     152            CALL iom_rstput( kt, nitbkg_r, inum, 'c_to_chl'      , cchl_p(:,:,1)         ) 
     153#elif defined key_medusa && defined key_foam_medusa 
     154            CALL iom_rstput( kt, nitbkg_r, inum, 'pgrow_avg'     , pgrow_avg                           ) 
     155            CALL iom_rstput( kt, nitbkg_r, inum, 'ploss_avg'     , ploss_avg                           ) 
     156            CALL iom_rstput( kt, nitbkg_r, inum, 'phyt_avg'      , phyt_avg                            ) 
     157            CALL iom_rstput( kt, nitbkg_r, inum, 'mld_max'       , mld_max                             ) 
     158            CALL iom_rstput( kt, nitbkg_r, inum, 'nutrients'     , trn(:,:,:,jpdin)                    ) 
     159            CALL iom_rstput( kt, nitbkg_r, inum, 'phytoplankton' , trn(:,:,:,jpphn) + trn(:,:,:,jpphd) ) 
     160            CALL iom_rstput( kt, nitbkg_r, inum, 'zooplankton'   , trn(:,:,:,jpzmi) + trn(:,:,:,jpzme) ) 
     161            CALL iom_rstput( kt, nitbkg_r, inum, 'detritus'      , trn(:,:,:,jpdet)                    ) 
     162            CALL iom_rstput( kt, nitbkg_r, inum, 'dic'           , trn(:,:,:,jpdic)                    ) 
     163            CALL iom_rstput( kt, nitbkg_r, inum, 'alkalinity'    , trn(:,:,:,jpalk)                    ) 
     164            CALL iom_rstput( kt, nitbkg_r, inum, 'chlorophyll'   , trn(:,:,1,jpchn) + trn(:,:,1,jpchd) ) 
     165#endif 
    123166            ! 
    124167            CALL iom_close( inum ) 
Note: See TracChangeset for help on using the changeset viewer.