23 |
c | |
c | |
24 |
c |-- exf_getclim (get climatological fields used e.g. for relax.) |
c |-- exf_getclim (get climatological fields used e.g. for relax.) |
25 |
c | |
c | |
26 |
c |-- exf_getsurfacefluxes |
c |-- exf_getffields <- this one does almost everything |
|
c | | |
|
|
c | |-- exf_getffields <- this one does almost everything |
|
27 |
c | | 1. reads in fields, either flux or atmos. state, |
c | | 1. reads in fields, either flux or atmos. state, |
28 |
c | | depending on CPP options |
c | | depending on CPP options |
29 |
c | | 2. If forcing and control is flux, we're already done here |
c | | 2. If forcing and control is flux, we're already done here |
36 |
c | | * ctrl_getvwind |
c | | * ctrl_getvwind |
37 |
c | | 4. flux fields are interpolated on current time |
c | | 4. flux fields are interpolated on current time |
38 |
c | | (b) next the flux fields are computed via bulk formulae |
c | | (b) next the flux fields are computed via bulk formulae |
|
c | | 5. if forcing and control is flux, then the control vector |
|
|
c | | anomalies are read here |
|
|
c | | * ctrl_getheatflux |
|
|
c | | * ctrl_getsaltflux |
|
|
c | | * ctrl_getzonstress |
|
|
c | | * call ctrl_getmerstress |
|
39 |
c | | |
c | | |
40 |
c | |-- exf_obcs |
c | |-- exf_obcs |
41 |
c | | If open boundaries are prescribed externally, |
c | If open boundaries are prescribed externally, |
42 |
c | | OB values for current time step are read here |
c | OB values for current time step are read here |
43 |
c | | 1. for each boundary (north, south, west, east) |
c | 1. for each boundary (north, south, west, east) |
44 |
c | | sliced fields are read for T,S,U,V |
c | sliced fields are read for T,S,U,V |
45 |
c | | 2. if OB's are part of control vector, |
c | 2. if OB's are part of control vector, |
46 |
c | | control variable anomalies are added to OB values |
c | control variable anomalies are added to OB values |
47 |
c | | * ctrl_getobcsn |
c | * ctrl_getobcsn |
48 |
c | | * ctrl_getobcss |
c | * ctrl_getobcss |
49 |
c | | * ctrl_getobcsw |
c | * ctrl_getobcsw |
50 |
c | | * ctrl_getobcse |
c | * ctrl_getobcse |
51 |
c | | |
c | |
52 |
c | |-- exf_mapfields |
c |-- exf_getsurfacefluxes |
53 |
c | | Forcing fields from exf package are mapped onto |
c | If forcing and control is flux, then the |
54 |
c | | mitgcm forcing arrays. |
c | control vector anomalies are read here |
55 |
c | | Mapping enables a runtime rescaling of fields |
c | * ctrl_getheatflux |
56 |
|
c | * ctrl_getsaltflux |
57 |
|
c | * ctrl_getzonstress |
58 |
|
c | * call ctrl_getmerstress |
59 |
|
c | |
60 |
|
c |-- exf_mapfields |
61 |
|
c | Forcing fields from exf package are mapped onto |
62 |
|
c | mitgcm forcing arrays. |
63 |
|
c | Mapping enables a runtime rescaling of fields |
64 |
c |
c |
65 |
c \ev |
c \ev |
66 |
CEOI |
CEOI |
84 |
c this way. A positive zonal wind stress accelerates the |
c this way. A positive zonal wind stress accelerates the |
85 |
c model ocean towards the east. |
c model ocean towards the east. |
86 |
c started: eckert@mit.edu, heimbach@mit.edu, ralf@ocean.mit.edu |
c started: eckert@mit.edu, heimbach@mit.edu, ralf@ocean.mit.edu |
87 |
|
c changed: menemenlis@jpl.nas.gov |
88 |
c *================================================================= |
c *================================================================= |
89 |
c | SUBROUTINE exf_getforcing |
c | SUBROUTINE exf_getforcing |
90 |
c *================================================================= |
c *================================================================= |