/[lmdze]/trunk/dyn3d/gcm.f90
ViewVC logotype

Contents of /trunk/dyn3d/gcm.f90

Parent Directory Parent Directory | Revision Log Revision Log


Revision 344 - (show annotations)
Tue Nov 12 15:18:14 2019 UTC (4 years, 7 months ago) by guez
File size: 3117 byte(s)
Replace pi / 180 by `deg_to_rad`

In procedure etat0, rename variable tsoil to ftsoil, which is the
corresponding name in the gcm program.

In `laplacien_gam`, replace call to scopy by array assignment.

Replace pi / 180 by `deg_to_rad` in `start_init_phys`.

Encapsulate diagcld1 and orolift in modules.

Avoid duplicated computation in `interfsurf_hq`.

Promote internal function fz of procedure soil to function of module
`soil_m`.  Use `new_unit` in procedure soil.

1 PROGRAM gcm
2
3 ! Authors: P. Le Van, L. Fairhead, F. Hourdin
4 ! From "gcm.F", version 1.4, 2006/04/04
5
6 ! General circulation model of LMD. Avec coordonn\'ee verticale
7 ! hybride, avec nouveaux op\'erateurs de dissipation "*" ("gradiv2",
8 ! "divgrad2", "nxgraro2"). Possibilit\'e de choisir le sch\'ema pour
9 ! l'advection de "q", en modifiant "iadv" dans "traceur.def".
10
11 ! Libraries:
12 use netcdf95, only: nf95_close
13
14 use comconst, only: dtvr, iniconst
15 use comdissnew, only: read_comdissnew
16 use comgeom, only: aire_2d, inigeom
17 use comgeomphy, only: airephy
18 use conf_gcm_m, only: day_step, iperiod, iphysiq, nday, conf_gcm, iflag_phys
19 use conf_guide_m, only: conf_guide
20 use dimensions, only: iim, jjm, llm, nqmx
21 USE disvert_m, ONLY : disvert
22 use dynetat0_m, only: dynetat0, day_ini
23 use dynetat0_chosen_m, only: dynetat0_chosen
24 use dynredem0_m, only: dynredem0
25 use grid_change, only: dyn_phy, init_dyn_phy
26 use histclo_m, only: histclo
27 use infotrac_init_m, only: infotrac_init
28 use inidissip_m, only: inidissip
29 use inifilr_m, only: inifilr
30 use inithist_m, only: inithist
31 use init_dynzon_m, only: init_dynzon
32 USE ioconf_calendar_m, only: ioconf_calendar
33 use leapfrog_m, only: leapfrog
34 use suphec_m, only: suphec
35 use unit_nml_m, only: unit_nml, set_unit_nml
36 use createnewfield_m, only: NbField, Ncid
37
38 IMPLICIT NONE
39
40 ! Variables dynamiques :
41 REAL ucov(iim + 1, jjm + 1, llm), vcov(iim + 1, jjm, llm) ! vent covariant
42 REAL teta(iim + 1, jjm + 1, llm) ! temp\'erature potentielle
43 REAL q(iim + 1, jjm + 1, llm, nqmx) ! mass fraction of advected species
44 REAL ps(iim + 1, jjm + 1) ! pression au sol (Pa)
45 REAL masse(iim + 1, jjm + 1, llm) ! masse d'air
46 REAL phis(iim + 1, jjm + 1) ! g\'eopotentiel au sol
47
48 LOGICAL:: true_calendar = .false. ! default value
49 integer i
50
51 namelist /main_nml/true_calendar
52
53 !------------------------------------------------------------
54
55 call set_unit_nml
56 open(unit_nml, file="used_namelists.txt", status="replace", action="write")
57
58 CALL conf_gcm
59 call read_comdissnew
60
61 print *, "Enter namelist 'main_nml'."
62 read (unit=*, nml=main_nml)
63 write(unit_nml, nml=main_nml)
64
65 ! Choix du calendrier :
66 if (true_calendar) then
67 call ioconf_calendar('gregorian')
68 else
69 call ioconf_calendar('360d')
70 endif
71
72 call infotrac_init
73 CALL iniconst
74 CALL dynetat0_chosen
75 CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)
76 CALL disvert
77 CALL inigeom ! initialisation de la g\'eometrie
78 CALL inifilr ! initialisation du filtre
79 CALL inidissip
80 call init_dyn_phy
81
82 ! Initialisation de la physique :
83 IF (iflag_phys) THEN
84 airephy = pack(aire_2d, dyn_phy)
85 CALL suphec
86 ENDIF
87
88 ! Initialisation des entr\'ees-sorties :
89 CALL dynredem0(day_ini + nday, phis)
90 CALL inithist(t_ops = dtvr, t_wrt = dtvr)
91 call init_dynzon(dt_app = dtvr * iperiod)
92
93 CALL conf_guide
94 CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)
95
96 close(unit_nml)
97 call histclo
98
99 do i = 1, nbfield
100 call nf95_close(Ncid(i))
101 end do
102
103 print *, 'Simulation finished'
104 print *, 'Everything is cool'
105
106 END PROGRAM gcm

  ViewVC Help
Powered by ViewVC 1.1.21