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

Diff of /trunk/dyn3d/gcm.f

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

revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC revision 263 by guez, Wed Mar 7 14:41:46 2018 UTC
# Line 3  PROGRAM gcm Line 3  PROGRAM gcm
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 15:05:16
5    
6    ! General circulation model of LMD. Avec coordonnée verticale    ! General circulation model of LMD. Avec coordonn\'ee verticale
7    ! hybride, avec nouveaux opérateurs de dissipation "*" ("gradiv2",    ! hybride, avec nouveaux op\'erateurs de dissipation "*" ("gradiv2",
8    ! "divgrad2", "nxgraro2"). Possibilité de choisir le schéma 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 calendar, only: ioconf_calendar    use comconst, only: dtvr, iniconst
12    use comconst, only: daysec, dtvr, iniconst    use comdissnew, only: read_comdissnew
13    use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv, inigeom    use comgeom, only:  aire_2d, inigeom
14    use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond    use comgeomphy, only: airephy
15    use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, &    use conf_gcm_m, only: day_step, iperiod, iecri, iphysiq, nday, conf_gcm, &
16         nday, raz_date, periodav, conf_gcm, iflag_phys         iflag_phys
17      use conf_guide_m, only: conf_guide
18    use dimens_m, only: iim, jjm, llm, nqmx    use dimens_m, only: iim, jjm, llm, nqmx
   use dimphy, only: klon  
19    USE disvert_m, ONLY : disvert    USE disvert_m, ONLY : disvert
20    use dynetat0_m, only: dynetat0, day_ini    use dynetat0_m, only: dynetat0, day_ini
21    use dynredem0_m, only: dynredem0    use dynredem0_m, only: dynredem0
# Line 24  PROGRAM gcm Line 24  PROGRAM gcm
24    use iniadvtrac_m, only: iniadvtrac    use iniadvtrac_m, only: iniadvtrac
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
30    use jumble, only: new_unit    use jumble, only: new_unit
31    use leapfrog_m, only: leapfrog    use leapfrog_m, only: leapfrog
32      use netcdf95, only: nf95_close
33    use suphec_m, only: suphec    use suphec_m, only: suphec
   use temps, only: day_ref, annee_ref, day_end, itau_dyn  
   use tracstoke, only: istdyn, istphy  
34    use unit_nml_m, only: unit_nml    use unit_nml_m, only: unit_nml
35    use yoethf_m, only: yoethf    use yoethf_m, only: yoethf
36      use createnewfield_m, only: NbField, Ncid
37    
38    IMPLICIT NONE    IMPLICIT NONE
39    
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érature potentielle    REAL teta(iim + 1, jjm + 1, llm) ! temp\'erature potentielle
43    REAL q(iim + 1, jjm + 1, llm, nqmx) ! champs advectés    REAL q(iim + 1, jjm + 1, llm, nqmx) ! champs advect\'es
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éopotentiel au sol    REAL phis(iim + 1, jjm + 1) ! g\'eopotentiel au sol
47    
   ! Variables pour le fichier histoire :  
   REAL time_0 ! time in day, as a fraction of day, in [0, 1[  
   
   ! Calendrier :  
48    LOGICAL:: true_calendar = .false. ! default value    LOGICAL:: true_calendar = .false. ! default value
49      integer i
   logical mask_v(iim + 1, jjm)  
   ! (mask for points in the "v" grid, first index is for longitude,  
   ! second index is for latitude)  
50    
51    namelist /main_nml/true_calendar    namelist /main_nml/true_calendar
52    
# Line 63  PROGRAM gcm Line 56  PROGRAM gcm
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
59      call read_comdissnew
60    
61    print *, "Enter namelist 'main_nml'."    print *, "Enter namelist 'main_nml'."
62    read (unit=*, nml=main_nml)    read (unit=*, nml=main_nml)
# Line 75  PROGRAM gcm Line 69  PROGRAM gcm
69       call ioconf_calendar('360d')       call ioconf_calendar('360d')
70    endif    endif
71    
   ! Initialisation des traceurs  
   ! Choix du schéma pour l'advection dans le fichier "traceur.def" ou via INCA  
72    call iniadvtrac    call iniadvtrac
   
73    CALL iniconst    CALL iniconst
74      CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis)
   ! Lecture du fichier "start.nc" :  
   CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)  
   
   ! On remet le calendrier à zero si demandé :  
   if (annee_ref /= anneeref .or. day_ref /= dayref) then  
      print *, 'Attention : les dates initiales lues dans le fichier ' // &  
           '"start" ne correspondent pas à celles lues dans "gcm.def".'  
      if (raz_date) then  
         print *, 'On réinitialise à la date lue dans "gcm.def".'  
         annee_ref = anneeref  
         day_ref = dayref  
         day_ini = dayref  
         itau_dyn = 0  
         time_0 = 0.  
      else  
         print *, 'On garde les dates du fichier "start".'  
      endif  
   ELSE  
      raz_date = .false.  
   endif  
   
75    CALL disvert    CALL disvert
76    CALL inigeom ! initialisation de la géometrie    CALL inigeom ! initialisation de la g\'eometrie
77    CALL inifilr ! initialisation du filtre    CALL inifilr ! initialisation du filtre
78    CALL inidissip    CALL inidissip
79    call init_dyn_phy    call init_dyn_phy
80    
81    ! Initialisation de la physique :    ! Initialisation de la physique :
82    IF (iflag_phys == 1) THEN    IF (iflag_phys) THEN
      rlatd(1)=rlatu(1)  
      rlatd(2:klon-1) = pack(spread(rlatu(2:jjm), 1, iim), .true.)  
      rlatd(klon)= rlatu(jjm + 1)  
   
      rlond(1)=0.  
      rlond(2:klon-1) = pack(spread(rlonv(:iim), 2, jjm - 1), .true.)  
      rlond(klon)= 0.  
   
      cuphy = pack(cu_2d, dyn_phy)  
   
      ! Construct a mask for points in the "v" grid:  
      mask_v = .true.  
      mask_v(2:, 1) = .false.  
      mask_v(iim + 1, 2:) = .false.  
   
      cvphy(:klon - 1) = pack(cv_2d, mask_v)  
      cvphy(klon) = cv_2d(1, jjm)  
      ! (that value of "cv_2d" is used twice in "cvphy")  
   
83       airephy = pack(aire_2d, dyn_phy)       airephy = pack(aire_2d, dyn_phy)
84       CALL suphec       CALL suphec
85       call yoethf       call yoethf
86    ENDIF    ENDIF
87    
88    ! Initialisation des entrées-sorties :    ! Initialisation des entr\'ees-sorties :
89    day_end = day_ini + nday    CALL dynredem0(day_ini + nday, phis)
90    print *, "day_ini = ", day_ini    CALL inithist(t_ops = iecri * dtvr, t_wrt = iecri * dtvr)
   print *, "day_end = ", day_end  
   
   CALL dynredem0("restart.nc", day_end, phis)  
   CALL inithist(day_ref, annee_ref, dtvr, nqmx, t_ops = iecri * daysec, &  
        t_wrt = iecri * daysec)  
   CALL initdynav(day_ref, annee_ref, dtvr, nqmx, t_ops = iperiod * dtvr, &  
        t_wrt = periodav * daysec)  
91    call init_dynzon(dt_app = dtvr * iperiod)    call init_dynzon(dt_app = dtvr * iperiod)
92    
93    ! Choix des fréquences de stockage pour le hors-ligne :    CALL conf_guide
94    istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q)
   istphy = istdyn / iphysiq      
   
   ! Intégration temporelle du modèle :  
   CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)  
95    
96    close(unit_nml)    close(unit_nml)
97    call histclo    call histclo
98    
99      do i = 1, nbfield
100         call nf95_close(Ncid(i))
101      end do
102    
103    print *, 'Simulation finished'    print *, 'Simulation finished'
104    print *, 'Everything is cool'    print *, 'Everything is cool'
105    
Legend:
Removed from v.82  
changed lines
  Added in v.263
  ViewVC Help
Powered by ViewVC 1.1.21