--- trunk/dyn3d/gcm.f 2014/09/16 14:00:41 108 +++ trunk/dyn3d/gcm.f 2015/02/12 16:23:33 128 @@ -13,6 +13,7 @@ use comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond use conf_gcm_m, only: day_step, iperiod, anneeref, dayref, iecri, iphysiq, & nday, raz_date, periodav, conf_gcm, iflag_phys + use conf_guide_m, only: conf_guide use dimens_m, only: iim, jjm, llm, nqmx use dimphy, only: klon USE disvert_m, ONLY : disvert @@ -35,7 +36,7 @@ use tracstoke, only: istdyn, istphy use unit_nml_m, only: unit_nml use yoethf_m, only: yoethf - use write_field, only: NbField, FieldId + use createnewfield_m, only: NbField, Ncid IMPLICIT NONE @@ -86,7 +87,7 @@ CALL iniconst ! Lecture du fichier "start.nc" : - CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0) + CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis) ! On remet le calendrier à zéro si demandé : if (raz_date) then @@ -95,7 +96,6 @@ day_ref = dayref day_ini = dayref itau_dyn = 0 - time_0 = 0. else print *, 'On garde les dates du fichier "start".' endif @@ -148,14 +148,16 @@ istdyn = day_step / 4 ! stockage toutes les 6 h = 1 jour / 4 istphy = istdyn / iphysiq + CALL conf_guide + ! Intégration temporelle du modèle : - CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0) + CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q) close(unit_nml) call histclo do i = 1, nbfield - call nf95_close(FieldId(i)) + call nf95_close(Ncid(i)) end do print *, 'Simulation finished'