/[MITgcm]/MITgcm/model/inc/FFIELDS.h
ViewVC logotype

Contents of /MITgcm/model/inc/FFIELDS.h

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

Revision 1.9 - (show annotations) (download)
Wed Nov 29 22:29:23 2000 UTC (23 years, 5 months ago) by adcroft
Branch: MAIN
CVS Tags: branch-atmos-merge-freeze, branch-atmos-merge-start, branch-atmos-merge-shapiro, checkpoint33, checkpoint34, branch-atmos-merge-zonalfilt, branch-atmos-merge-phase5, branch-atmos-merge-phase4, branch-atmos-merge-phase7, branch-atmos-merge-phase6, branch-atmos-merge-phase1, branch-atmos-merge-phase3, branch-atmos-merge-phase2
Branch point for: branch-atmos-merge
Changes since 1.8: +17 -27 lines
File MIME type: text/plain 
Fixed confusion about units of forcing arrays in FFIELDS.h
namely Fu,Fv,Qnet,Qsw,EmPmR:
  - Removed verification/natl_box/code/external_fields_scale.F
        (did not differ from that in model/src)
  - Changed units of fu,fv,Qnet,Qsw,EmPmR back to proper units
     (see FFIELDS.h for description)
  - Scale fu,fv,Qnet,Qsw,EmPmR when used in external_forcing_surf.F,
    kpp_calc.F and kpp_transport_t.F
  - Removed model/src/external_fields_scale.F and calls to it
  - verification/natl_box uses flux data with "atmospheric" sign so
    a special version of external_fields_load.F is used to
    change the data as it's read in. This way, the arrays
    have the right units and signs at all times tha a user could
    possibly use them.
1 C $Header: /u/gcmpack/models/MITgcmUV/model/inc/FFIELDS.h,v 1.8 2000/11/13 16:32:57 heimbach Exp $
2 C
3 C /==========================================================\
4 C | FFIELDS.h |
5 C | o Model forcing fields |
6 C |==========================================================|
7 C | The arrays here will need changing and customising for a |
8 C | particular experiment. |
9 C \==========================================================/
10 C
11 C-- For a classical "gyre" type experiment just one term is needed.
12 C
13 C fu - Zonal surface wind stress
14 C Units are N/m^2 (>0 from East to West)
15 C
16 C fv - Meridional surface wind stress
17 C Units are N/m^2 (>0 from North to South))
18 C
19 C EmPmR - Evaporation - Precipitation - Runoff
20 C Units are m/s (>0 for ocean salting)
21 C
22 C Qnet - Upward surface heat flux
23 C Units are W/m^2=kg/s^3 (>0 for ocean cooling)
24 C
25 C Qsw - Upward short-wave surface heat flux
26 C Units are W/m^2=kg/s^3 (>0 for ocean cooling)
27 C
28 C SST - Sea surface temperature (degrees) for relaxation
29 C SSS - Sea surface salinity (psu) for relaxation
30
31 COMMON /FFIELDS/
32 & fu,
33 & fv,
34 & Qnet,
35 & EmPmR,
36 & SST,
37 & SSS,
38 & Qsw
39 _RS fu (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
40 _RS fv (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
41 _RS Qnet (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
42 _RS EmPmR (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
43 _RS SST (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
44 _RS SSS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
45 _RS Qsw (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
46
47 C surfaceTendencyU (units are m/s^2)
48 C -> usage in gU: gU = gU + surfaceTendencyU[m/s^2]
49 C
50 C surfaceTendencyV (units are m/s^2)
51 C -> usage in gV: gV = gV + surfaceTendencyV[m/s^2]
52 C
53 C surfaceTendencyS (units are psu/s)
54 C - EmPmR plus salinity relaxation term
55 C -> calculate -lambda*(S(model)-S(clim))
56 C -> usage in gS: gS = gS + surfaceTendencyS[psu/s]
57 C
58 C surfaceTendencyT (units are degrees/s)
59 C - Qnet plus temp. relaxation
60 C -> calculate -lambda*(T(model)-T(clim))
61 C >>> Qnet assumed to be total flux minus s/w rad. <<<
62 C -> usage in gT: gT = gT + surfaceTendencyT[K/s]
63 C
64 COMMON /TENDENCY_FORCING/
65 & surfaceTendencyU,
66 & surfaceTendencyV,
67 & surfaceTendencyT,
68 & surfaceTendencyS,
69 & tempQsw
70 _RS surfaceTendencyU (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
71 _RS surfaceTendencyV (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
72 _RS surfaceTendencyT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
73 _RS surfaceTendencyS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
74 _RS tempQsw (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
  ViewVC Help
Powered by ViewVC 1.1.22