/[MITgcm]/MITgcm/pkg/bling/bling_airseaflux.F
ViewVC logotype

Diff of /MITgcm/pkg/bling/bling_airseaflux.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by mmazloff, Sat Oct 15 21:30:43 2016 UTC revision 1.5 by mmazloff, Thu May 11 13:32:24 2017 UTC
# Line 66  C     FluxO2           :: air-sea flux o Line 66  C     FluxO2           :: air-sea flux o
66  C     (FluxCO2 is a global variable)  C     (FluxCO2 is a global variable)
67        _RL  SGDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  SGDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
68        _RL  SGO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  SGO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
69        _RL  FluxO2      (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  FluxO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
   
70    
71  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
72    
# Line 76  C     i,j              :: Loop counters Line 75  C     i,j              :: Loop counters
75        INTEGER i,j,klev        INTEGER i,j,klev
76  C Number of iterations for pCO2 solvers  C Number of iterations for pCO2 solvers
77        _RL co3dummy        _RL co3dummy
78        _RL Kwexch_Pre  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL Kwexch_Pre   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
79  C Solubility relation coefficients  C Solubility relation coefficients
80        _RL SchmidtNoDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL SchmidtNoDIC (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
81        _RL pCO2sat     (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL pCO2sat      (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
82        _RL Kwexch      (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL Kwexch       (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
83        _RL pisvel      (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL pisvel       (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
84  C local variables for carbon chem  C local variables for carbon chem
85        _RL surfalk     (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL surfalk      (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
86        _RL surfphos    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL surfphos     (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
87        _RL surfsi      (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL surfsi       (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
88        _RL surftemp    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL surftemp     (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
89        _RL surfsalt    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL surfsalt     (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
90        _RL surfdic     (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL surfdic      (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
91  C o2 solubility relation coefficients  C o2 solubility relation coefficients
92        _RL SchmidtNoO2 (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL SchmidtNoO2  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
93        _RL O2sat       (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL O2sat        (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
94        _RL Kwexch_o2   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL O2sat_percent(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
95          _RL Kwexch_o2    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
96        _RL aTT        _RL aTT
97        _RL aTK        _RL aTK
98        _RL aTS        _RL aTS
# Line 283  C using Garcia and Gordon (1992), L&O (m Line 283  C using Garcia and Gordon (1992), L&O (m
283  c Convert from ml/l to mol/m^3  c Convert from ml/l to mol/m^3
284                O2sat(i,j) = o2s/22391.6 _d 0 * 1. _d 3                O2sat(i,j) = o2s/22391.6 _d 0 * 1. _d 3
285    
286    cav
287                  O2sat_percent(i,j) = PTR_O2(i,j,1)/O2sat(i,j)*100
288    
289  C Determine flux, inc. correction for local atmos surface pressure  C Determine flux, inc. correction for local atmos surface pressure
290  #ifdef USE_EXF_ATMPRES  #ifdef USE_EXF_ATMPRES
291  C Atm pressure in N/m2, convert to bars  C Atm pressure in N/m2, convert to bars
# Line 312  C update surface tendencies Line 315  C update surface tendencies
315          _EXCH_XY_RL( pCO2, mythid)          _EXCH_XY_RL( pCO2, mythid)
316          _EXCH_XYZ_RL( pH, mythid)          _EXCH_XYZ_RL( pH, mythid)
317    
318    #ifdef ALLOW_DIAGNOSTICS
319          IF ( useDiagnostics ) THEN
320            CALL DIAGNOSTICS_FILL(O2sat_percent,'BLGO2SAT',0,1,2,bi,bj,
321         &       myThid)
322          ENDIF
323    #endif /* ALLOW_DIAGNOSTICS */
324    
325  #endif /* ALLOW_PTRACER */  #endif /* ALLOW_PTRACER */
326    
327          RETURN          RETURN
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22