Parent Directory | Revision Log
Moved the call to read_serre out of conf_gcm so that it can be called only in the program ce0l, not in gcm. In gcm, variables of module serre are read from start file. Added reading of dzoomx, dzoomy, taux, tauy from start file, in dynetat0. Those variables were written by dynredem0 but not read. Removed possibility fxyhypb = false, because the geometric part of the program is such a mess. Could then remove variables transx, transy, alphax, alphay, pxo, pyo of module serre. Bug fix in tau2alpha: missing save attributes. The first call to tau2alpha needs to compute dxdyu and dxdyv regardless of value of argument type, because they will be needed for subsequent calls to tau2alpha with various values of argument type.
1 | guez | 68 | PROGRAM ce0l |
2 | guez | 3 | |
3 | ! This program sets the initial and boundary values. | ||
4 | |||
5 | use conf_gcm_m, only: conf_gcm | ||
6 | guez | 107 | use dimens_m, only: iim, jjm |
7 | guez | 3 | use etat0_mod, only: etat0 |
8 | guez | 107 | use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub |
9 | guez | 57 | use jumble, only: new_unit |
10 | guez | 3 | use limit_mod, only: limit |
11 | guez | 113 | use serre, only: read_serre |
12 | guez | 57 | use unit_nml_m, only: unit_nml |
13 | guez | 3 | |
14 | implicit none | ||
15 | |||
16 | guez | 107 | REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2 |
17 | |||
18 | guez | 3 | !------------------------------------- |
19 | |||
20 | guez | 57 | call new_unit(unit_nml) |
21 | guez | 79 | open(unit_nml, file="used_namelists.txt", status="replace", action="write") |
22 | guez | 3 | CALL conf_gcm |
23 | guez | 113 | call read_serre |
24 | guez | 107 | CALL etat0(phis) |
25 | guez | 3 | CALL limit |
26 | guez | 107 | CALL grilles_gcm_netcdf_sub(phis) |
27 | guez | 57 | close(unit_nml) |
28 | guez | 3 | |
29 | guez | 68 | END PROGRAM ce0l |
ViewVC Help | |
Powered by ViewVC 1.1.21 |