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

Diff of /trunk/dyn3d/gcm.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/dyn3d/gcm.f revision 265 by guez, Tue Mar 20 09:35:59 2018 UTC trunk/dyn3d/gcm.f90 revision 344 by guez, Tue Nov 12 15:18:14 2019 UTC
# Line 1  Line 1 
1  PROGRAM gcm  PROGRAM gcm
2    
3    ! Authors: P. Le Van, L. Fairhead, F. Hourdin    ! Authors: P. Le Van, L. Fairhead, F. Hourdin
4    ! From "gcm.F", version 1.4, 2006/04/04 15:05:16    ! From "gcm.F", version 1.4, 2006/04/04
5    
6    ! General circulation model of LMD. Avec coordonn\'ee verticale    ! General circulation model of LMD. Avec coordonn\'ee verticale
7    ! hybride, avec nouveaux op\'erateurs de dissipation "*" ("gradiv2",    ! hybride, avec nouveaux op\'erateurs de dissipation "*" ("gradiv2",
8    ! "divgrad2", "nxgraro2"). Possibilit\'e de choisir le sch\'ema pour    ! "divgrad2", "nxgraro2"). Possibilit\'e de choisir le sch\'ema pour
9    ! l'advection de "q", en modifiant "iadv" dans "traceur.def".    ! 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    use comconst, only: dtvr, iniconst
15    use comdissnew, only: read_comdissnew    use comdissnew, only: read_comdissnew
16    use comgeom, only:  aire_2d, inigeom    use comgeom, only:  aire_2d, inigeom
17    use comgeomphy, only: airephy    use comgeomphy, only: airephy
18    use conf_gcm_m, only: day_step, iperiod, iecri, iphysiq, nday, conf_gcm, &    use conf_gcm_m, only: day_step, iperiod, iphysiq, nday, conf_gcm, iflag_phys
        iflag_phys  
19    use conf_guide_m, only: conf_guide    use conf_guide_m, only: conf_guide
20    use dimensions, only: iim, jjm, llm, nqmx    use dimensions, only: iim, jjm, llm, nqmx
21    USE disvert_m, ONLY : disvert    USE disvert_m, ONLY : disvert
22    use dynetat0_m, only: dynetat0, day_ini    use dynetat0_m, only: dynetat0, day_ini
23      use dynetat0_chosen_m, only: dynetat0_chosen
24    use dynredem0_m, only: dynredem0    use dynredem0_m, only: dynredem0
25    use grid_change, only: dyn_phy, init_dyn_phy    use grid_change, only: dyn_phy, init_dyn_phy
26    use histclo_m, only: histclo    use histclo_m, only: histclo
27    use iniadvtrac_m, only: iniadvtrac    use infotrac_init_m, only: infotrac_init
28    use inidissip_m, only: inidissip    use inidissip_m, only: inidissip
29    use inifilr_m, only: inifilr    use inifilr_m, only: inifilr
30    use inithist_m, only: inithist    use inithist_m, only: inithist
31    use init_dynzon_m, only: init_dynzon    use init_dynzon_m, only: init_dynzon
32    USE ioconf_calendar_m, only: ioconf_calendar    USE ioconf_calendar_m, only: ioconf_calendar
   use jumble, only: new_unit  
33    use leapfrog_m, only: leapfrog    use leapfrog_m, only: leapfrog
   use netcdf95, only: nf95_close  
34    use suphec_m, only: suphec    use suphec_m, only: suphec
35    use unit_nml_m, only: unit_nml    use unit_nml_m, only: unit_nml, set_unit_nml
   use yoethf_m, only: yoethf  
36    use createnewfield_m, only: NbField, Ncid    use createnewfield_m, only: NbField, Ncid
37    
38    IMPLICIT NONE    IMPLICIT NONE
# Line 40  PROGRAM gcm Line 40  PROGRAM gcm
40    ! Variables dynamiques :    ! Variables dynamiques :
41    REAL ucov(iim + 1, jjm + 1, llm), vcov(iim + 1, jjm, llm)  ! vent covariant    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    REAL teta(iim + 1, jjm + 1, llm) ! temp\'erature potentielle
43    REAL q(iim + 1, jjm + 1, llm, nqmx) ! champs advect\'es    REAL q(iim + 1, jjm + 1, llm, nqmx) ! mass fraction of advected species
44    REAL ps(iim + 1, jjm + 1) ! pression au sol (Pa)    REAL ps(iim + 1, jjm + 1) ! pression au sol (Pa)
45    REAL masse(iim + 1, jjm + 1, llm) ! masse d'air    REAL masse(iim + 1, jjm + 1, llm) ! masse d'air
46    REAL phis(iim + 1, jjm + 1) ! g\'eopotentiel au sol    REAL phis(iim + 1, jjm + 1) ! g\'eopotentiel au sol
# Line 52  PROGRAM gcm Line 52  PROGRAM gcm
52    
53    !------------------------------------------------------------    !------------------------------------------------------------
54    
55    call new_unit(unit_nml)    call set_unit_nml
56    open(unit_nml, file="used_namelists.txt", status="replace", action="write")    open(unit_nml, file="used_namelists.txt", status="replace", action="write")
57    
58    CALL conf_gcm    CALL conf_gcm
# Line 69  PROGRAM gcm Line 69  PROGRAM gcm
69       call ioconf_calendar('360d')       call ioconf_calendar('360d')
70    endif    endif
71    
72    call iniadvtrac    call infotrac_init
73    CALL iniconst    CALL iniconst
74      CALL dynetat0_chosen
75    CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)    CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)
76    CALL disvert    CALL disvert
77    CALL inigeom ! initialisation de la g\'eometrie    CALL inigeom ! initialisation de la g\'eometrie
# Line 82  PROGRAM gcm Line 83  PROGRAM gcm
83    IF (iflag_phys) THEN    IF (iflag_phys) THEN
84       airephy = pack(aire_2d, dyn_phy)       airephy = pack(aire_2d, dyn_phy)
85       CALL suphec       CALL suphec
      call yoethf  
86    ENDIF    ENDIF
87    
88    ! Initialisation des entr\'ees-sorties :    ! Initialisation des entr\'ees-sorties :
89    CALL dynredem0(day_ini + nday, phis)    CALL dynredem0(day_ini + nday, phis)
90    CALL inithist(t_ops = iecri * dtvr, t_wrt = iecri * dtvr)    CALL inithist(t_ops = dtvr, t_wrt = dtvr)
91    call init_dynzon(dt_app = dtvr * iperiod)    call init_dynzon(dt_app = dtvr * iperiod)
92    
93    CALL conf_guide    CALL conf_guide
Legend:
Removed from v.265  
changed lines
  Added in v.344
  ViewVC Help
Powered by ViewVC 1.1.21