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

Diff of /trunk/dyn3d/gcm.f

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.f revision 254 by guez, Mon Feb 5 10:39:38 2018 UTC
# Line 9  PROGRAM gcm Line 9  PROGRAM gcm
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: daysec, 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, iecri, iphysiq, nday, periodav, &
# Line 31  PROGRAM gcm Line 32  PROGRAM gcm
32    use leapfrog_m, only: leapfrog    use leapfrog_m, only: leapfrog
33    use netcdf95, only: nf95_close    use netcdf95, only: nf95_close
34    use suphec_m, only: suphec    use suphec_m, only: suphec
   use tracstoke, only: istdyn, istphy  
35    use unit_nml_m, only: unit_nml    use unit_nml_m, only: unit_nml
36    use yoethf_m, only: yoethf    use yoethf_m, only: yoethf
37    use createnewfield_m, only: NbField, Ncid    use createnewfield_m, only: NbField, Ncid
# Line 57  PROGRAM gcm Line 57  PROGRAM gcm
57    open(unit_nml, file="used_namelists.txt", status="replace", action="write")    open(unit_nml, file="used_namelists.txt", status="replace", action="write")
58    
59    CALL conf_gcm    CALL conf_gcm
60      call read_comdissnew
61    
62    print *, "Enter namelist 'main_nml'."    print *, "Enter namelist 'main_nml'."
63    read (unit=*, nml=main_nml)    read (unit=*, nml=main_nml)
# Line 79  PROGRAM gcm Line 80  PROGRAM gcm
80    call init_dyn_phy    call init_dyn_phy
81    
82    ! Initialisation de la physique :    ! Initialisation de la physique :
83    IF (iflag_phys == 1) THEN    IF (iflag_phys) THEN
84       airephy = pack(aire_2d, dyn_phy)       airephy = pack(aire_2d, dyn_phy)
85       CALL suphec       CALL suphec
86       call yoethf       call yoethf
# Line 91  PROGRAM gcm Line 92  PROGRAM gcm
92    CALL initdynav(dtvr, nqmx, t_ops = iperiod * dtvr, t_wrt = periodav * daysec)    CALL initdynav(dtvr, nqmx, t_ops = iperiod * dtvr, t_wrt = periodav * daysec)
93    call init_dynzon(dt_app = dtvr * iperiod)    call init_dynzon(dt_app = dtvr * iperiod)
94    
   ! 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      
   
95    CALL conf_guide    CALL conf_guide
96    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)
97    
Legend:
Removed from v.178  
changed lines
  Added in v.254
  ViewVC Help
Powered by ViewVC 1.1.21