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

Contents of /trunk/dyn3d/serre.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: 1096 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 module serre
2
3 implicit none
4
5 REAL:: clon = 0. ! longitude of the center of the zoom, in degrees
6 real:: clat = 0. ! latitude of the center of the zoom, in degrees
7
8 real:: grossismx = 1. ! facteur de grossissement du zoom, selon la longitude
9 real:: grossismy = 1. ! facteur de grossissement du zoom, selon la latitude
10
11 real:: dzoomx = 0.
12 ! extension en longitude de la zone du zoom (fraction de la zone totale)
13
14 real:: dzoomy = 0.
15 ! extension en latitude de la zone du zoom (fraction de la zone totale)
16
17 real:: taux = 3. ! raideur du zoom en X
18 real:: tauy = 3. ! raideur du zoom en Y
19
20 contains
21
22 subroutine read_serre
23
24 use unit_nml_m, only: unit_nml
25 use nr_util, only: assert
26
27 namelist /serre_nml/ clon, clat, grossismx, grossismy, dzoomx, dzoomy, &
28 taux, tauy
29
30 !-------------------------------------------------
31
32 print *, "Enter namelist 'serre_nml'."
33 read(unit=*, nml=serre_nml)
34 write(unit_nml, nml=serre_nml)
35
36 call assert(grossismx >= 1. .and. grossismy >= 1., "read_serre grossism")
37
38 end subroutine read_serre
39
40 end module serre
  ViewVC Help
Powered by ViewVC 1.1.21