--- trunk/libf/dyn3d/gcm.f90 2012/04/20 14:58:43 61 +++ trunk/dyn3d/gcm.f90 2014/02/05 17:51:07 78 @@ -9,9 +9,8 @@ ! l'advection de "q", en modifiant "iadv" dans "traceur.def". USE calendar, only: ioconf_calendar - use clesphys2, only: read_clesphys2 use comconst, only: daysec, cpp, dtvr, g, rad, r - use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv + use comgeom, only: rlatu, aire_2d, cu_2d, cv_2d, rlonv, inigeom 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 @@ -24,7 +23,6 @@ use iniadvtrac_m, only: iniadvtrac use inidissip_m, only: inidissip use inifilr_m, only: inifilr - use inigeom_m, only: inigeom use initdynav_m, only: initdynav use inithist_m, only: inithist use init_dynzon_m, only: init_dynzon @@ -78,9 +76,6 @@ call ioconf_calendar('360d') endif - ! Lecture des fichiers "gcm.def" ou "run.def" : - call read_clesphys2 - ! Initialisation des traceurs ! Choix du schéma pour l'advection dans le fichier "traceur.def" ou via INCA call iniadvtrac @@ -118,7 +113,6 @@ ENDIF CALL iniconst - close(unit_nml) CALL inigeom ! initialisation de la géometrie CALL inifilr ! initialisation du filtre CALL inidissip @@ -169,6 +163,7 @@ ! Intégration temporelle du modèle : CALL leapfrog(ucov, vcov, teta, ps, masse, phis, q, time_0) + close(unit_nml) call histclo print *, 'Simulation finished' print *, 'Everything is cool'