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 10302 for branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90 – NEMO

Ignore:
Timestamp:
2018-11-13T18:21:16+01:00 (5 years ago)
Author:
dford
Message:

Merge in revisions 8447:10159 of dev_r5518_GO6_package.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90

    r8132 r10302  
    177177   INTEGER  ::  jdms_model   !: choice of DMS model passed to atmosphere 
    178178!!                              1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL 
    179 !! 
     179!! JPALM --19-12-2017 -- UM people need to tune the Anderson DMS scheme     
     180   REAL(wp) ::  dmsmin       !: Anderson DMS scheme - DMS minimum value 
     181   REAL(wp) ::  dmscut       !: Anderson DMS scheme - DMS cutoff value  
     182   REAL(wp) ::  dmsslp       !: Anderson DMS scheme - DMS slope  
     183!! FOR UKESM    
     184   REAL(wp) ::  scl_chl      !: scaling factor for tuned Chl passed to the UM  
     185   INTEGER  ::  chl_out      !: select Chl field exported and scaled for UM: 
     186                             !: 1- Surface Chl ; 2- MLD Chl 
    180187!! 
    181188   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   remdmp   !: depth dependent damping coefficient of passive tracers  
    … …  
    205212   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   :: f2_ccd_arg  !: 2D aragonite CCD depth 
    206213!! 
     214!! 2D fields of pCO2 and fCO2 for observation operator 
     215   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   :: f2_pco2w    !: 2D pCO2 
     216   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   :: f2_fco2w    !: 2D fCO2 
     217!! 
    207218!! 2D fields of organic and inorganic material sedimented on the seafloor 
    208219   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   :: zb_sed_n    !: 2D organic nitrogen   (before) 
    … …  
    276287#if defined key_roam 
    277288!!---------------------------------------------------------------------- 
    278 !! Atmospheric pCO2 data (1859 to 2100 inclusive) 
    279 !!---------------------------------------------------------------------- 
    280 !! 
    281    REAL(wp), DIMENSION(242)         ::   hist_pco2 !: pCO2 
    282  
    283 # if defined key_rcp26 
    284       !! UKMO, run AJKKH + KAAEC, RCP 2.6, pCO2 time evolution 
    285       DATA hist_pco2 / 286.0230, 286.1730, 286.3230, 286.4480, 286.5730, & 
    286       & 286.7230, 286.8480, 286.9480, 287.0480, 287.1730, & 
    287       & 287.3230, 287.4730, 287.6480, 287.8480, 288.0730, & 
    288       & 288.3480, 288.6480, 288.9730, 289.3470, 289.7470, & 
    289       & 290.1730, 290.6470, 291.1470, 291.6220, 292.0720, & 
    290       & 292.5220, 292.9220, 293.2470, 293.5220, 293.7470, & 
    291       & 293.9470, 294.1220, 294.2720, 294.4220, 294.5470, & 
    292       & 294.6470, 294.7470, 294.8470, 294.9710, 295.1710, & 
    293       & 295.4460, 295.7470, 296.0720, 296.4210, 296.7710, & 
    294       & 297.1460, 297.5710, 298.0210, 298.4460, 298.8460, & 
    295       & 299.2460, 299.6450, 300.0210, 300.3710, 300.7200, & 
    296       & 301.0450, 301.3460, 301.6710, 302.0200, 302.3450, & 
    297       & 302.6450, 302.9700, 303.3450, 303.7200, 304.0700, & 
    298       & 304.4700, 304.9200, 305.3440, 305.7700, 306.2450, & 
    299       & 306.7190, 307.1700, 307.6440, 308.1190, 308.5440, & 
    300       & 308.9440, 309.3440, 309.6940, 309.9440, 310.1190, & 
    301       & 310.2440, 310.3190, 310.3190, 310.2440, 310.1440, & 
    302       & 310.0690, 310.0440, 310.0690, 310.1440, 310.2690, & 
    303       & 310.4440, 310.6940, 311.0430, 311.4440, 311.8690, & 
    304       & 312.3680, 312.9430, 313.5430, 314.1680, 314.7900, & 
    305       & 315.4430, 316.2150, 317.0170, 317.7370, 318.3400, & 
    306       & 318.8680, 319.5900, 320.5890, 321.5470, 322.5770, & 
    307       & 323.8440, 324.9260, 325.7960, 327.0810, 328.6180, & 
    308       & 329.6830, 330.5250, 331.6880, 333.2120, 334.7870, & 
    309       & 336.4640, 338.2990, 339.6660, 340.7310, 342.1360, & 
    310       & 343.7200, 345.2200, 346.7350, 348.5820, 350.6740, & 
    311       & 352.4230, 353.7910, 354.9530, 355.8210, 356.7130, & 
    312       & 358.0630, 359.7720, 361.3970, 363.0900, 365.2560, & 
    313       & 367.2810, 368.7980, 370.4000, 372.4550, 374.6920, & 
    314       & 376.7440, 378.7440, 380.7580, 382.7080, 384.7300, & 
    315       & 386.9310, 389.2150, 391.4910, 393.7710, 396.0460, & 
    316       & 398.3240, 400.6080, 402.8950, 405.1780, 407.4550, & 
    317       & 409.7260, 411.9930, 414.2500, 416.4410, 418.5280, & 
    318       & 420.5250, 422.4390, 424.2720, 426.0200, 427.6750, & 
    319       & 429.2360, 430.7050, 432.0850, 433.3580, 434.5140, & 
    320       & 435.5740, 436.5490, 437.4420, 438.2550, 438.9810, & 
    321       & 439.6110, 440.1430, 440.5770, 440.9450, 441.2660, & 
    322       & 441.5410, 441.7840, 442.0050, 442.2040, 442.3780, & 
    323       & 442.5210, 442.6200, 442.6720, 442.6810, 442.6540, & 
    324       & 442.5830, 442.4670, 442.3270, 442.1680, 441.9960, & 
    325       & 441.8060, 441.5930, 441.3440, 441.0540, 440.7230, & 
    326       & 440.3510, 439.9300, 439.4650, 438.9730, 438.4630, & 
    327       & 437.9400, 437.4020, 436.8400, 436.2640, 435.6850, & 
    328       & 435.1030, 434.5160, 433.9170, 433.3060, 432.7010, & 
    329       & 432.1110, 431.5380, 430.9810, 430.4320, 429.8860, & 
    330       & 429.3370, 428.7810, 428.2220, 427.6490, 427.0660, & 
    331       & 426.4890, 425.9270, 425.3840, 424.8610, 424.3540, & 
    332       & 423.8540, 423.3540, 422.8530, 422.3510, 421.8410, & 
    333       & 421.3250, 420.8190 / 
    334 # else 
    335       !! UKMO, run AJKKH + KAAEF, RCP 8.5, pCO2 time evolution 
    336       DATA hist_pco2 / 286.0230, 286.1730, 286.3230, 286.4480, 286.5730, & 
    337       & 286.7230, 286.8480, 286.9480, 287.0480, 287.1730, & 
    338       & 287.3230, 287.4730, 287.6480, 287.8480, 288.0730, & 
    339       & 288.3480, 288.6480, 288.9730, 289.3470, 289.7470, & 
    340       & 290.1730, 290.6470, 291.1470, 291.6220, 292.0720, & 
    341       & 292.5220, 292.9220, 293.2470, 293.5220, 293.7470, & 
    342       & 293.9470, 294.1220, 294.2720, 294.4220, 294.5470, & 
    343       & 294.6470, 294.7470, 294.8470, 294.9710, 295.1710, & 
    344       & 295.4460, 295.7470, 296.0720, 296.4210, 296.7710, & 
    345       & 297.1460, 297.5710, 298.0210, 298.4460, 298.8460, & 
    346       & 299.2460, 299.6450, 300.0210, 300.3710, 300.7200, & 
    347       & 301.0450, 301.3460, 301.6710, 302.0200, 302.3450, & 
    348       & 302.6450, 302.9700, 303.3450, 303.7200, 304.0700, & 
    349       & 304.4700, 304.9200, 305.3440, 305.7700, 306.2450, & 
    350       & 306.7190, 307.1700, 307.6440, 308.1190, 308.5440, & 
    351       & 308.9440, 309.3440, 309.6940, 309.9440, 310.1190, & 
    352       & 310.2440, 310.3190, 310.3190, 310.2440, 310.1440, & 
    353       & 310.0690, 310.0440, 310.0690, 310.1440, 310.2690, & 
    354       & 310.4440, 310.6940, 311.0430, 311.4440, 311.8690, & 
    355       & 312.3680, 312.9430, 313.5430, 314.1680, 314.7900, & 
    356       & 315.4430, 316.2150, 317.0170, 317.7370, 318.3400, & 
    357       & 318.8680, 319.5900, 320.5890, 321.5470, 322.5770, & 
    358       & 323.8440, 324.9260, 325.7960, 327.0810, 328.6180, & 
    359       & 329.6830, 330.5250, 331.6880, 333.2120, 334.7870, & 
    360       & 336.4640, 338.2990, 339.6660, 340.7310, 342.1360, & 
    361       & 343.7200, 345.2200, 346.7350, 348.5820, 350.6740, & 
    362       & 352.4230, 353.7910, 354.9530, 355.8210, 356.7130, & 
    363       & 358.0630, 359.7720, 361.3970, 363.0900, 365.2560, & 
    364       & 367.2810, 368.7980, 370.4000, 372.4550, 374.6920, & 
    365       & 376.7440, 378.7440, 380.7580, 382.7080, 384.7300, & 
    366       & 386.9420, 389.2540, 391.5670, 393.9370, 396.3920, & 
    367       & 398.9320, 401.5550, 404.2550, 407.0220, 409.8530, & 
    368       & 412.7470, 415.7050, 418.7210, 421.7880, 424.9180, & 
    369       & 428.1200, 431.3970, 434.7470, 438.1650, 441.6410, & 
    370       & 445.1700, 448.7530, 452.3920, 456.0950, 459.8810, & 
    371       & 463.7680, 467.7660, 471.8750, 476.0960, 480.4210, & 
    372       & 484.8390, 489.3470, 493.9430, 498.6400, 503.4380, & 
    373       & 508.3410, 513.3630, 518.5160, 523.8050, 529.2290, & 
    374       & 534.7780, 540.4450, 546.2230, 552.1120, 558.1110, & 
    375       & 564.2110, 570.4130, 576.7390, 583.1990, 589.7980, & 
    376       & 596.5390, 603.4110, 610.4060, 617.4940, 624.6500, & 
    377       & 631.8800, 639.1750, 646.5360, 653.9800, 661.5230, & 
    378       & 669.1840, 676.9570, 684.8290, 692.7790, 700.7690, & 
    379       & 708.8050, 716.8870, 725.0020, 733.1770, 741.3900, & 
    380       & 749.6700, 758.0480, 766.5050, 775.0350, 783.6110, & 
    381       & 792.2200, 800.8740, 809.5680, 818.2760, 827.0090, & 
    382       & 835.8020, 844.6550, 853.5730, 862.5690, 871.6190, & 
    383       & 880.7020, 889.8240, 898.9590, 908.1270, 917.3080, & 
    384       & 926.4960, 935.7040 / 
    385 # endif 
     289!! JPALM -- change hist_pco2 init 
     290!!---------------------------------------------------------------------- 
     291!! 
     292   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:)   ::   hist_pco2 !: pCO2 
     293   INTEGER  :: co2_rec 
     294   REAL(wp) :: co2_yinit, co2_yend    !: First and Last year read in the xCO2.atm file 
     295   REAL(wp) :: xobs_xco2a             !: Observed atmospheric xCO2, from namelist 
     296#endif 
     297 
     298!!---------------------------------------------------------------------- 
     299!! JPALM -- PI CO2 key 
     300!!---------------------------------------------------------------------- 
     301!! 
     302#if defined key_axy_pi_co2 
     303   LOGICAL , PUBLIC ::   lk_pi_co2 = .TRUE.  !: PI xCO2 used 
     304#else 
     305   LOGICAL , PUBLIC ::   lk_pi_co2 = .FALSE. !: PI xCO2 unused 
    386306#endif 
    387307 
    … …  
    434354   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   :: cmask       !: ??? 
    435355 
     356!!---------------------------------------------------------------------- 
     357!! Parameters required for ocean colour assimilation 
     358!!---------------------------------------------------------------------- 
     359!! 
     360   LOGICAL :: ln_foam_medusa                                 !: Flag to calculate and save diagnostics 
     361   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: pgrow_avg  !: Mixed layer average phytoplankton growth 
     362   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: ploss_avg  !: Mixed layer average phytoplankton loss 
     363   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: phyt_avg   !: Mixed layer average phytoplankton 
     364   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: mld_max    !: Maximum mixed layer depth 
     365!! 
     366 
    436367   !!---------------------------------------------------------------------- 
    437368   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    … …  
    446377      !!---------------------------------------------------------------------- 
    447378      USE lib_mpp , ONLY: ctl_warn 
    448       INTEGER ::   ierr(8)        ! Local variables 
     379      INTEGER ::   ierr(9)        ! Local variables 
    449380      !!---------------------------------------------------------------------- 
    450381      ierr(:) = 0 
    … …  
    456387      !* 2D and 3D fields of carbonate system parameters 
    457388      ALLOCATE( f2_ccd_cal(jpi,jpj)  , f2_ccd_arg(jpi,jpj)  ,       & 
     389                f2_pco2w(jpi,jpj)    , f2_fco2w(jpi,jpj)    ,       & 
    458390         &      f3_pH(jpi,jpj,jpk)   , f3_h2co3(jpi,jpj,jpk),       & 
    459391         &      f3_hco3(jpi,jpj,jpk) , f3_co3(jpi,jpj,jpk)  ,       & 
    … …  
    504436         &      ffln(jpi,jpj,jpk)    , fflf(jpi,jpj,jpk)    ,       & 
    505437         &      ffls(jpi,jpj,jpk)    , cmask(jpi,jpj)       ,    STAT=ierr(8) )  
     438      !* Fields for ocean colour data assimilation 
     439      ALLOCATE( pgrow_avg(jpi,jpj)   , ploss_avg(jpi,jpj)   ,       & 
     440         &      phyt_avg(jpi,jpj)    , mld_max(jpi,jpj)     ,    STAT=ierr(9) ) 
    506441#endif 
    507442      ! 
Note: See TracChangeset for help on using the changeset viewer.