/[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/Sources/dyn3d/gcm.f revision 178 by guez, Fri Mar 11 18:47:26 2016 UTC trunk/dyn3d/gcm.f90 revision 328 by guez, Thu Jun 13 14:40:06 2019 UTC
# Line 8  PROGRAM gcm Line 8  PROGRAM gcm
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    use comconst, only: daysec, dtvr, iniconst    use comconst, only: dtvr, iniconst
12      use comdissnew, only: read_comdissnew
13    use comgeom, only:  aire_2d, inigeom    use comgeom, only:  aire_2d, inigeom
14    use comgeomphy, only: airephy    use comgeomphy, only: airephy
15    use conf_gcm_m, only: day_step, iperiod, iecri, iphysiq, nday, periodav, &    use conf_gcm_m, only: day_step, iperiod, iphysiq, nday, conf_gcm, iflag_phys
        conf_gcm, iflag_phys  
16    use conf_guide_m, only: conf_guide    use conf_guide_m, only: conf_guide
17    use dimens_m, only: iim, jjm, llm, nqmx    use dimensions, only: iim, jjm, llm, nqmx
18    USE disvert_m, ONLY : disvert    USE disvert_m, ONLY : disvert
19    use dynetat0_m, only: dynetat0, day_ini    use dynetat0_m, only: dynetat0, day_ini
20      use dynetat0_chosen_m, only: dynetat0_chosen
21    use dynredem0_m, only: dynredem0    use dynredem0_m, only: dynredem0
22    use grid_change, only: dyn_phy, init_dyn_phy    use grid_change, only: dyn_phy, init_dyn_phy
23    use histclo_m, only: histclo    use histclo_m, only: histclo
24    use iniadvtrac_m, only: iniadvtrac    use infotrac_init_m, only: infotrac_init
25    use inidissip_m, only: inidissip    use inidissip_m, only: inidissip
26    use inifilr_m, only: inifilr    use inifilr_m, only: inifilr
   use initdynav_m, only: initdynav  
27    use inithist_m, only: inithist    use inithist_m, only: inithist
28    use init_dynzon_m, only: init_dynzon    use init_dynzon_m, only: init_dynzon
29    USE ioconf_calendar_m, only: ioconf_calendar    USE ioconf_calendar_m, only: ioconf_calendar
   use jumble, only: new_unit  
30    use leapfrog_m, only: leapfrog    use leapfrog_m, only: leapfrog
31    use netcdf95, only: nf95_close    use netcdf95, only: nf95_close
32    use suphec_m, only: suphec    use suphec_m, only: suphec
33    use tracstoke, only: istdyn, istphy    use unit_nml_m, only: unit_nml, set_unit_nml
   use unit_nml_m, only: unit_nml  
   use yoethf_m, only: yoethf  
34    use createnewfield_m, only: NbField, Ncid    use createnewfield_m, only: NbField, Ncid
35    
36    IMPLICIT NONE    IMPLICIT NONE
# Line 53  PROGRAM gcm Line 50  PROGRAM gcm
50    
51    !------------------------------------------------------------    !------------------------------------------------------------
52    
53    call new_unit(unit_nml)    call set_unit_nml
54    open(unit_nml, file="used_namelists.txt", status="replace", action="write")    open(unit_nml, file="used_namelists.txt", status="replace", action="write")
55    
56    CALL conf_gcm    CALL conf_gcm
57      call read_comdissnew
58    
59    print *, "Enter namelist 'main_nml'."    print *, "Enter namelist 'main_nml'."
60    read (unit=*, nml=main_nml)    read (unit=*, nml=main_nml)
# Line 69  PROGRAM gcm Line 67  PROGRAM gcm
67       call ioconf_calendar('360d')       call ioconf_calendar('360d')
68    endif    endif
69    
70    call iniadvtrac    call infotrac_init
71    CALL iniconst    CALL iniconst
72      CALL dynetat0_chosen
73    CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)    CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)
74    CALL disvert    CALL disvert
75    CALL inigeom ! initialisation de la g\'eometrie    CALL inigeom ! initialisation de la g\'eometrie
# Line 79  PROGRAM gcm Line 78  PROGRAM gcm
78    call init_dyn_phy    call init_dyn_phy
79    
80    ! Initialisation de la physique :    ! Initialisation de la physique :
81    IF (iflag_phys == 1) THEN    IF (iflag_phys) THEN
82       airephy = pack(aire_2d, dyn_phy)       airephy = pack(aire_2d, dyn_phy)
83       CALL suphec       CALL suphec
      call yoethf  
84    ENDIF    ENDIF
85    
86    ! Initialisation des entr\'ees-sorties :    ! Initialisation des entr\'ees-sorties :
87    CALL dynredem0(day_ini + nday, phis)    CALL dynredem0(day_ini + nday, phis)
88    CALL inithist(dtvr, nqmx, t_ops = iecri * daysec, t_wrt = iecri * daysec)    CALL inithist(t_ops = dtvr, t_wrt = dtvr)
   CALL initdynav(dtvr, nqmx, t_ops = iperiod * dtvr, t_wrt = periodav * daysec)  
89    call init_dynzon(dt_app = dtvr * iperiod)    call init_dynzon(dt_app = dtvr * iperiod)
90    
   ! Choix des fr\'equences de stockage pour le hors-ligne :  
   istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4  
   istphy = istdyn / iphysiq      
   
91    CALL conf_guide    CALL conf_guide
92    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)
93    
Legend:
Removed from v.178  
changed lines
  Added in v.328
  ViewVC Help
Powered by ViewVC 1.1.21