13 |
C fu - Zonal velocity tendency term |
C fu - Zonal velocity tendency term |
14 |
C -> read assumes N/m^2 (>0 from East to West) |
C -> read assumes N/m^2 (>0 from East to West) |
15 |
C -> transformed to m/s^2 |
C -> transformed to m/s^2 |
16 |
C via fu -> fu/(rhoNil*dR) |
C via fu -> -fu/(rhoNil*dR) |
|
C -> usage in gU: gU = gU + fu[m/s^2] |
|
17 |
C |
C |
18 |
C fv - Meridional velocity tendency term |
C fv - Meridional velocity tendency term |
19 |
C -> read assumes N/m^2 (>0 from North to South)) |
C -> read assumes N/m^2 (>0 from North to South)) |
20 |
C -> transformed to m/s^2 |
C -> transformed to m/s^2 |
21 |
C via fv -> fv/(rhoNil*dR) |
C via fv -> -fv/(rhoNil*dR) |
|
C -> usage in gU: gV = gV + fu[m/s^2] |
|
22 |
C |
C |
23 |
C EmPmR - Evaporation - Precipitation - Runoff |
C EmPmR - Evaporation - Precipitation - Runoff |
24 |
C -> read assumes m/s (>0 for ocean salting) |
C -> read assumes m/s (>0 for ocean salting) |
25 |
C -> transformed to psu/s |
C -> transformed to psu/s |
26 |
C via empmr -> -empmr*35./dR |
C via empmr -> empmr*35./dR |
|
C -> usage in gS: gS = gS + empmr[psu/s] |
|
27 |
C |
C |
28 |
C Qnet - Surface heat flux |
C Qnet - Surface heat flux |
29 |
C -> read assumes W/m^2=kg/s^3 (>0 for ocean cooling) |
C -> read assumes W/m^2=kg/s^3 (>0 for ocean cooling) |
30 |
C -> transformed to K/s |
C -> transformed to K/s |
31 |
C via Qnet -> -Qnet/(rhonil*Cp*dR) |
C via Qnet -> -Qnet/(rhonil*Cp*dR) |
|
C -> usage in gT: gT = gT + qnet[K/s] |
|
32 |
C |
C |
33 |
C Qsw - Short-wave surface heat flux |
C Qsw - Short-wave surface heat flux |
34 |
C -> read assumes W/m^2=kg/s^3 (>0 for ocean cooling) |
C -> read assumes W/m^2=kg/s^3 (>0 for ocean cooling) |
35 |
C -> transformed to K/s |
C -> transformed to K/s |
36 |
C via Qsw -> -Qsw/(rhonil*Cp*dR) |
C via Qsw -> -Qsw/(rhonil*Cp*dR) |
|
C -> usage in gT: gT = gT + Qswt[K/s] |
|
|
C only for #define SHORTWAVE_HEATING |
|
37 |
C |
C |
38 |
C SST - Sea surface temperature (degrees) for relaxation |
C SST - Sea surface temperature (degrees) for relaxation |
39 |
C SSS - Sea surface salinity (psu) for relaxation |
C SSS - Sea surface salinity (psu) for relaxation |
70 |
C -> calculate -lambda*(T(model)-T(clim)) |
C -> calculate -lambda*(T(model)-T(clim)) |
71 |
C >>> Qnet assumed to be total flux minus s/w rad. <<< |
C >>> Qnet assumed to be total flux minus s/w rad. <<< |
72 |
C -> usage in gT: gT = gT + surfaceTendencyT[K/s] |
C -> usage in gT: gT = gT + surfaceTendencyT[K/s] |
73 |
|
C |
74 |
COMMON /TENDENCY_FORCING/ |
COMMON /TENDENCY_FORCING/ |
75 |
& surfaceTendencyU, |
& surfaceTendencyU, |
76 |
& surfaceTendencyV, |
& surfaceTendencyV, |
77 |
& surfaceTendencyT, |
& surfaceTendencyT, |
78 |
& surfaceTendencyS |
& surfaceTendencyS, |
79 |
|
& tempQsw |
80 |
_RS surfaceTendencyU (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS surfaceTendencyU (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
81 |
_RS surfaceTendencyV (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS surfaceTendencyV (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
82 |
_RS surfaceTendencyT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS surfaceTendencyT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
83 |
_RS surfaceTendencyS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS surfaceTendencyS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
84 |
|
_RS tempQsw (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |