--- trunk/libf/dyn3d/gcm.f90 2012/01/30 12:54:02 57 +++ trunk/libf/dyn3d/gcm.f90 2013/02/18 16:33:12 69 @@ -9,7 +9,6 @@ ! 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 comgeomphy, only: airephy, cuphy, cvphy, rlatd, rlond @@ -20,7 +19,7 @@ use dynetat0_m, only: dynetat0, day_ini use dynredem0_m, only: dynredem0 use grid_change, only: dyn_phy, init_dyn_phy - use histcom, only: histclo + use histclo_m, only: histclo use iniadvtrac_m, only: iniadvtrac use inidissip_m, only: inidissip use inifilr_m, only: inifilr @@ -78,9 +77,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 +114,6 @@ ENDIF CALL iniconst - close(unit_nml) CALL inigeom ! initialisation de la géometrie CALL inifilr ! initialisation du filtre CALL inidissip @@ -169,6 +164,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'