/[lmdze]/trunk/dyn3d/ce0l.f
ViewVC logotype

Contents of /trunk/dyn3d/ce0l.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 113 - (show annotations)
Thu Sep 18 19:56:46 2014 UTC (9 years, 7 months ago) by guez
File size: 712 byte(s)
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 PROGRAM ce0l
2
3 ! This program sets the initial and boundary values.
4
5 use conf_gcm_m, only: conf_gcm
6 use dimens_m, only: iim, jjm
7 use etat0_mod, only: etat0
8 use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub
9 use jumble, only: new_unit
10 use limit_mod, only: limit
11 use serre, only: read_serre
12 use unit_nml_m, only: unit_nml
13
14 implicit none
15
16 REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2
17
18 !-------------------------------------
19
20 call new_unit(unit_nml)
21 open(unit_nml, file="used_namelists.txt", status="replace", action="write")
22 CALL conf_gcm
23 call read_serre
24 CALL etat0(phis)
25 CALL limit
26 CALL grilles_gcm_netcdf_sub(phis)
27 close(unit_nml)
28
29 END PROGRAM ce0l

  ViewVC Help
Powered by ViewVC 1.1.21