/[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 108 by guez, Tue Sep 16 14:00:41 2014 UTC
# Line 8  PROGRAM gcm Line 8  PROGRAM gcm
8    ! "divgrad2", "nxgraro2"). Possibilité de choisir le schéma pour    ! "divgrad2", "nxgraro2"). Possibilité de choisir le schéma pour
9    ! l'advection de "q", en modifiant "iadv" dans "traceur.def".    ! l'advection de "q", en modifiant "iadv" dans "traceur.def".
10    
   USE calendar, only: ioconf_calendar  
11    use comconst, only: daysec, dtvr, iniconst    use comconst, only: daysec, dtvr, iniconst
12    use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv, inigeom    use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv, inigeom
13    use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond    use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond
# Line 27  PROGRAM gcm Line 26  PROGRAM gcm
26    use initdynav_m, only: initdynav    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
34    use temps, only: day_ref, annee_ref, day_end, itau_dyn    use temps, only: day_ref, annee_ref, day_end, itau_dyn
35    use tracstoke, only: istdyn, istphy    use tracstoke, only: istdyn, istphy
36    use unit_nml_m, only: unit_nml    use unit_nml_m, only: unit_nml
37    use yoethf_m, only: yoethf    use yoethf_m, only: yoethf
38      use write_field, only: NbField, FieldId
39    
40    IMPLICIT NONE    IMPLICIT NONE
41    
# Line 55  PROGRAM gcm Line 57  PROGRAM gcm
57    ! (mask for points in the "v" grid, first index is for longitude,    ! (mask for points in the "v" grid, first index is for longitude,
58    ! second index is for latitude)    ! second index is for latitude)
59    
60      integer i
61    
62    namelist /main_nml/true_calendar    namelist /main_nml/true_calendar
63    
64    !------------------------------------------------------------    !------------------------------------------------------------
# Line 84  PROGRAM gcm Line 88  PROGRAM gcm
88    ! Lecture du fichier "start.nc" :    ! Lecture du fichier "start.nc" :
89    CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)    CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)
90    
91    ! On remet le calendrier à zero si demandé :    ! On remet le calendrier à zéro si demandé :
92    if (annee_ref /= anneeref .or. day_ref /= dayref) then    if (raz_date) then
93       print *, 'Attention : les dates initiales lues dans le fichier ' // &       print *, 'On réinitialise à la date lue dans la namelist.'
94            '"start" ne correspondent pas à celles lues dans "gcm.def".'       annee_ref = anneeref
95       if (raz_date) then       day_ref = dayref
96          print *, 'On réinitialise à la date lue dans "gcm.def".'       day_ini = dayref
97          annee_ref = anneeref       itau_dyn = 0
98          day_ref = dayref       time_0 = 0.
99          day_ini = dayref    else
100          itau_dyn = 0       print *, 'On garde les dates du fichier "start".'
         time_0 = 0.  
      else  
         print *, 'On garde les dates du fichier "start".'  
      endif  
   ELSE  
      raz_date = .false.  
101    endif    endif
102    
103    CALL disvert    CALL disvert
# Line 155  PROGRAM gcm Line 153  PROGRAM gcm
153    
154    close(unit_nml)    close(unit_nml)
155    call histclo    call histclo
156    
157      do i = 1, nbfield
158         call nf95_close(FieldId(i))
159      end do
160    
161    print *, 'Simulation finished'    print *, 'Simulation finished'
162    print *, 'Everything is cool'    print *, 'Everything is cool'
163    
Legend:
Removed from v.82  
changed lines
  Added in v.108
  ViewVC Help
Powered by ViewVC 1.1.21