/[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 108 by guez, Tue Sep 16 14:00:41 2014 UTC revision 115 by guez, Fri Sep 19 17:36:20 2014 UTC
# Line 13  PROGRAM gcm Line 13  PROGRAM gcm
13    use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond    use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond
14    use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, &    use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, &
15         nday, raz_date, periodav, conf_gcm, iflag_phys         nday, raz_date, periodav, conf_gcm, iflag_phys
16      use conf_guide_m, only: conf_guide
17    use dimens_m, only: iim, jjm, llm, nqmx    use dimens_m, only: iim, jjm, llm, nqmx
18    use dimphy, only: klon    use dimphy, only: klon
19    USE disvert_m, ONLY : disvert    USE disvert_m, ONLY : disvert
# Line 35  PROGRAM gcm Line 36  PROGRAM gcm
36    use tracstoke, only: istdyn, istphy    use tracstoke, only: istdyn, istphy
37    use unit_nml_m, only: unit_nml    use unit_nml_m, only: unit_nml
38    use yoethf_m, only: yoethf    use yoethf_m, only: yoethf
39    use write_field, only: NbField, FieldId    use createnewfield_m, only: NbField, Ncid
40    
41    IMPLICIT NONE    IMPLICIT NONE
42    
# Line 148  PROGRAM gcm Line 149  PROGRAM gcm
149    istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4    istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4
150    istphy = istdyn / iphysiq        istphy = istdyn / iphysiq    
151    
152      CALL conf_guide
153    
154    ! Intégration temporelle du modèle :    ! Intégration temporelle du modèle :
155    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0)
156    
# Line 155  PROGRAM gcm Line 158  PROGRAM gcm
158    call histclo    call histclo
159    
160    do i = 1, nbfield    do i = 1, nbfield
161       call nf95_close(FieldId(i))       call nf95_close(Ncid(i))
162    end do    end do
163    
164    print *, 'Simulation finished'    print *, 'Simulation finished'
Legend:
Removed from v.108  
changed lines
  Added in v.115
  ViewVC Help
Powered by ViewVC 1.1.21