/[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 7 by guez, Mon Mar 31 12:24:17 2008 UTC revision 10 by guez, Fri Apr 18 14:45:53 2008 UTC
# Line 15  PROGRAM gcm Line 15  PROGRAM gcm
15    ! Pour Van-Leer plus vapeur d'eau saturée : iadv(1)=4    ! Pour Van-Leer plus vapeur d'eau saturée : iadv(1)=4
16    ! Pour Van-Leer : iadv=10    ! Pour Van-Leer : iadv=10
17    
18    USE IOIPSL, only: ioconf_calendar    USE IOIPSL, only: ioconf_calendar, histclo
19    use dimens_m, only: iim, jjm, llm, nqmx    use dimens_m, only: iim, jjm, llm, nqmx
20    use dimphy, only: klon    use dimphy, only: klon
21    use paramet_m, only: ip1jm, ip1jmp1    use paramet_m, only: ip1jm, ip1jmp1
# Line 100  PROGRAM gcm Line 100  PROGRAM gcm
100    
101    ! Lecture du fichier "start.nc" :    ! Lecture du fichier "start.nc" :
102    CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)    CALL dynetat0(vcov, ucov, teta, q, masse, ps, phis, time_0)
   ! Begin special experiment  
 !!$     print *, "This is a special experiment."  
 !!$     print *, "We are setting:"  
 !!$     print *, "ucov = vcov = 0, q = 0, ps = 101325"  
 !!$     print *, "We are averaging 'teta' horizontally."  
 !!$     ucov = 0.  
 !!$     vcov = 0.  
 !!$     q = 0.  
 !!$     ps = 101325.  
   ! Average teta over all longitudes and latitudes:  
 !!$     forall(i = 1:llm) teta(:,i) = sum(teta(:,i)) / ip1jmp1  
   ! (it would be better to weight each element with an associated  
   ! surface area)  
   ! End special experiment  
103    
104    ! Lecture des paramètres de contrôle pour la simulation :    ! Lecture des paramètres de contrôle pour la simulation :
105    ! on recalcule éventuellement le pas de temps    ! on recalcule éventuellement le pas de temps
# Line 219  PROGRAM gcm Line 205  PROGRAM gcm
205    ! Intégration temporelle du modèle :    ! Intégration temporelle du modèle :
206    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, nq, q, clesphy0, time_0)    CALL leapfrog(ucov, vcov, teta, ps, masse, phis, nq, q, clesphy0, time_0)
207    
208      call histclo
209      print *, 'Simulation finished'
210      print *, 'Everything is cool'
211    
212  END PROGRAM gcm  END PROGRAM gcm
Legend:
Removed from v.7  
changed lines
  Added in v.10
  ViewVC Help
Powered by ViewVC 1.1.21