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 10047 for branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90 – NEMO

Ignore:
Timestamp:
2018-08-08T17:12:35+02:00 (6 years ago)
Author:
jpalmier
Message:

merge with GO6_package_branch 9385-10020 ; plus debug OMIP_DIC

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90

    r10045 r10047  
    5959                                   zalk, zdic, zoxy, zsal, ztmp,          & 
    6060#  if defined key_omip_dic 
    61                                    pi_fgco2,                               & 
    62                                    f_pi_co2flux, f_pi_co2flux2d,           & 
    63                                    f_pi_co2starair_2d, f_pi_co3,           & 
    64                                    f_pi_fco2a_2d, f_pi_fco2w_2d,           & 
    65                                    f_pi_h2co3, f_pi_hco3,                  & 
    66                                    f_pi_ocndpco2_2d,                       & 
    67                                    f_pi_omarg, f_pi_omcal,                 & 
    68                                    f_pi_pco2a2d, f_pi_pco2atm, f_pi_pco2w, & 
    69                                    f_pi_pco2w2d, f_pi_ph,                  & 
    70                                    f_pi_TDIC, f_pi_TDIC2d,                 & 
    71                                    f_pi_xco2a, f_pi_xco2a_2d,              & 
    72                                    zomd,                                   & 
    73 #  endif                                    
    74 # endif 
     61                                   pi_fgco2,                                 & 
     62                                   f_pi_co2flux, f_pi_co2flux_2d,            & 
     63                                   f_pi_co2starair_2d, f_pi_co3,             & 
     64                                   f_pi_fco2a_2d, f_pi_fco2w_2d,             & 
     65                                   f_pi_h2co3, f_pi_hco3,                    & 
     66                                   f_pi_ocndpco2_2d,                         & 
     67                                   f_pi_omarg, f_pi_omcal,                   & 
     68                                   f_pi_pco2a_2d, f_pi_pco2atm, f_pi_pco2w,  & 
     69                                   f_pi_pco2w_2d, f_pi_ph,                   & 
     70                                   f_pi_TDIC, f_pi_TDIC_2d,                  & 
     71                                   f_pi_xco2a, f_pi_xco2a_2d,                & 
     72               f_pi_ph_2d, f_pi_omcal_2d, f_pi_omarg_2d, & 
     73               f_pi_h2co3_2d, f_pi_hco3_2d, f_pi_co3_2d, & 
     74                                   zomd,                                     & 
     75#  endif   !! key_omip_dic                                 
     76# endif  !! key_roam 
    7577# if defined key_mocsy 
    7678                                   zpho,                                  & 
     79# endif  
     80                                   zchd, zchn, zdin, zsil 
     81      USE dom_oce,           ONLY: e3t_0, gphit, tmask, mig, mjg 
     82# if defined key_vvl 
     83      USE dom_oce,           ONLY: e3t_n 
    7784# endif 
    78                                    zchd, zchn, zdin, zsil 
    79       USE dom_oce,           ONLY: e3t_0, e3t_n, gphit, tmask, mig, mjg 
    80 # if defined key_iomput 
    8185      USE iom,               ONLY: lk_iomput 
    82 # endif 
    8386      USE in_out_manager,    ONLY: lwp, numout 
    8487      USE par_kind,          ONLY: wp 
    … …  
    101104#  else 
    102105      USE trcco2_medusa,     ONLY: trc_co2_medusa 
    103 #  endif 
     106#  endif  
    104107      USE trcdms_medusa,     ONLY: trc_dms_medusa 
    105108      USE trcoxy_medusa,     ONLY: trc_oxy_medusa 
    … …  
    136139#  if defined key_omip_dic 
    137140      !! AXY (06/08/18): output variables that are not different for OMIP DIC 
     141      REAL(wp), DIMENSION(jpi,jpj) :: f_pi_fco2w 
     142      REAL(wp), DIMENSION(jpi,jpj) :: f_pi_fco2atm 
     143      REAL(wp), DIMENSION(jpi,jpj) :: f_pi_co2starair, f_pi_dpco2 
    138144      REAL(wp) :: f_rhosw_dum, f_schmidtco2_dum, f_kwco2_dum,f_K0_dum 
    139145#  endif 
    … …  
    208214      call trc_rst_dia_stat(zpho(:,:), 'zpho') 
    209215#  endif 
    210  
    211216#  if defined key_axy_carbchem 
    212217#   if defined key_mocsy 
    … …  
    330335               !! mmol / m3 -> umol / kg 
    331336               f_pi_TDIC(ji,jj) = (zomd(ji,jj) / f_rhosw(ji,jj)) * 1000. 
    332 #    endif 
     337#    endif  !! omip 
    333338            ENDIF 
    334339         ENDDO 
    … …  
    355360#    endif 
    356361                
    357 #   else    
     362#   else   !! mocsy 
    358363      !!----------------------------------------------------------- 
    359364      !! Blackford et al. (2007) carbonate chemistry 
    … …  
    392397      ENDDO 
    393398 
    394 #   endif 
    395 #  else 
     399#   endif !! mocsy  
     400#  else  !! carbchem 
    396401      !!----------------------------------------------------------- 
    397402      !! No carbonate chemistry 
    … …  
    434439         ENDDO 
    435440      ENDDO 
    436 #  endif 
     441#  endif  !! carbchem  
    437442 
    438443#  if defined key_axy_killco2flux 
    … …  
    634639                  !! AXY (06/08/18): diagnostics for OMIP PI DIC 
    635640                  IF( med_diag%PI_ATM_PCO2%dgsave ) THEN 
    636                      f_pi_pco2a2d(ji,jj) = f_pi_pco2atm(ji,jj) 
    637                   ENDIF 
    638                   IF( med_diag%OCN_PCO2%dgsave ) THEN 
    639                      f_pi_pco2w2d(ji,jj) = f_pi_pco2w(ji,jj) 
     641                     f_pi_pco2a_2d(ji,jj) = f_pi_pco2atm(ji,jj) 
     642                  ENDIF 
     643                  IF( med_diag%PI_OCN_PCO2%dgsave ) THEN 
     644                     f_pi_pco2w_2d(ji,jj) = f_pi_pco2w(ji,jj) 
    640645                  ENDIF 
    641646                  IF( med_diag%PI_CO2FLUX%dgsave ) THEN 
    642647                     !! mmol/m3/d -> mmol/m2/d 
    643                      f_pi_co2flux2d(ji,jj) = f_pi_co2flux(ji,jj) * fse3t(ji,jj,1) 
     648                     f_pi_co2flux_2d(ji,jj) = f_pi_co2flux(ji,jj) * fse3t(ji,jj,1) 
    644649                  ENDIF 
    645650                  IF( med_diag%PI_TCO2%dgsave ) THEN 
    646                      f_PI_TDIC2d(ji,jj) = f_PI_TDIC(ji,jj) 
     651                     f_PI_TDIC_2d(ji,jj) = f_PI_TDIC(ji,jj) 
    647652                  ENDIF 
    648653                  IF( med_diag%PI_ATM_XCO2%dgsave ) THEN 
    … …  
    667672                          CO2flux_conv 
    668673                  ENDIF 
     674                  IF( med_diag%PI_OCN_PH%dgsave ) THEN 
     675                     f_pi_ph_2d(ji,jj) = f_pi_ph(ji,jj) 
     676                  ENDIF                   
     677                  IF( med_diag%PI_OCNH2CO3%dgsave ) THEN 
     678                     f_pi_h2co3_2d(ji,jj) = f_pi_h2co3(ji,jj) 
     679                  ENDIF                   
     680                  IF( med_diag%PI_OCN_HCO3%dgsave ) THEN 
     681                     f_pi_hco3_2d(ji,jj) = f_pi_hco3(ji,jj) 
     682                  ENDIF                   
     683                  IF( med_diag%PI_OCN_CO3%dgsave ) THEN 
     684                     f_pi_co3_2d(ji,jj) = f_pi_co3(ji,jj) 
     685                  ENDIF                   
     686                  IF( med_diag%PI_OM_CAL%dgsave ) THEN 
     687                     f_pi_omcal_2d(ji,jj) = f_pi_omcal(ji,jj) 
     688                  ENDIF                   
     689                  IF( med_diag%PI_OM_ARG%dgsave ) THEN 
     690                     f_pi_omarg_2d(ji,jj) = f_pi_omarg(ji,jj) 
     691                  ENDIF                   
    669692#   endif                   
    670693               ENDIF 
    … …  
    672695         ENDDO 
    673696      ENDDO 
    674  
    675 #  endif 
    676697 
    677698      !!----------------------------------------------------------------- 
    … …  
    766787         ENDDO 
    767788      ENDDO 
    768 # endif 
     789# endif  !! Roam 
    769790       
    770791   END SUBROUTINE air_sea 
    771792 
    772 #else 
     793#else  !! medusa 
    773794   !!====================================================================== 
    774795   !!  Dummy module :                                   No MEDUSA bio-model 
    … …  
    778799      WRITE(*,*) 'air_sea: You should not have seen this print! error?' 
    779800   END SUBROUTINE air_sea 
    780 #endif  
     801#endif  !! medusa 
    781802 
    782803   !!====================================================================== 
Note: See TracChangeset for help on using the changeset viewer.