Changeset 1320 for LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d

Timestamp:

Feb 25, 2010, 3:50:51 PM (15 years ago)

Author:

idelkadi

Message:

Decoupage du pas de temps physique en N sous-pas (splitting)

Location:

LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d

Files:

• calfis.F (modified) (3 diffs)
• conf_gcm.F (modified) (1 diff)
• control_mod.F90 (modified) (1 diff)
• gcm.F (modified) (1 diff)

LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/calfis.F

@@ -150 +150 @@
       REAL zdpsrf(ngridmx)
 c
+      REAL zdufic(ngridmx,llm),zdvfic(ngridmx,llm)
+      REAL zdtfic(ngridmx,llm),zdqfic(ngridmx,llm,nqtot)
+      REAL jH_cur_split,zdt_split
+      LOGICAL debut_split,lafin_split
+      INTEGER isplit
+
       REAL zsin(iim),zcos(iim),z1(iim)
       REAL zsinbis(iim),zcosbis(iim),z1bis(iim)
@@ -443 +449 @@
 
       if (planet_type=="earth") then
+
+      print*,'PHYSIQUE AVEC NSPLIT_PHYS=',nsplit_phys
+      zdt_split=dtphys/nsplit_phys
+      zdufic(:,:)=0.
+      zdvfic(:,:)=0.
+      zdtfic(:,:)=0.
+      zdqfic(:,:,:)=0.
+
+      do isplit=1,nsplit_phys
+
+         jH_cur_split=jH_cur+(isplit-1) * dtvr / (daysec *nsplit_phys)
+         debut_split=debut.and.isplit==1
+         lafin_split=lafin.and.isplit==nsplit_phys
+
 #ifdef CPP_EARTH
-      CALL physiq (ngridmx,
+         CALL physiq (ngridmx,
      .             llm,
-     .             debut,
-     .             lafin,
+     .             debut_split,
+     .             lafin_split,
      .             jD_cur,
-     .             jH_cur,
-     .             dtphys,
+     .             jH_cur_split,
+     .             zdt_split,
      .             zplev,
      .             zplay,
@@ -470 +490 @@
      .             pducov,
      .             PVteta)
+
+         zufi(:,:)=zufi(:,:)+zdufi(:,:)*zdt_split
+         zvfi(:,:)=zvfi(:,:)+zdvfi(:,:)*zdt_split
+         ztfi(:,:)=ztfi(:,:)+zdtfi(:,:)*zdt_split
+         zqfi(:,:,:)=zqfi(:,:,:)+zdqfi(:,:,:)*zdt_split
+
+         zdufic(:,:)=zdufic(:,:)+zdufi(:,:)
+         zdvfic(:,:)=zdvfic(:,:)+zdvfi(:,:)
+         zdtfic(:,:)=zdtfic(:,:)+zdtfi(:,:)
+         zdqfic(:,:,:)=zdqfic(:,:,:)+zdqfi(:,:,:)
+
+      enddo
+      zdufi(:,:)=zdufic(:,:)/nsplit_phys
+      zdvfi(:,:)=zdvfic(:,:)/nsplit_phys
+      zdtfi(:,:)=zdtfic(:,:)/nsplit_phys
+      zdqfi(:,:,:)=zdqfic(:,:,:)/nsplit_phys
+
 #endif
       endif !of if (planet_type=="earth")

LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/conf_gcm.F

@@ -163 +163 @@
        day_step = 240 
        CALL getin('day_step',day_step)
+
+!Config  Key  = nsplit_phys
+!Config  Desc = nombre de pas par jour
+!Config  Def  = 1 
+!Config  Help = nombre de pas par jour (multiple de iperiod) (
+!Config          ici pour  dt = 1 min ) 
+       nsplit_phys = 1 
+       CALL getin('nsplit_phys',nsplit_phys)
 
 !Config  Key  = iperiod

LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/control_mod.F90

@@ -11 +11 @@
 
   REAL    :: periodav
-  INTEGER :: nday,day_step,iperiod,iapp_tracvl
+  INTEGER :: nday,day_step,iperiod,iapp_tracvl,nsplit_phys
   INTEGER :: iconser,iecri,idissip,iphysiq,iecrimoy
   INTEGER :: dayref,anneeref, raz_date, ip_ebil_dyn

LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3d/gcm.F

@@ -406 +406 @@
          if (planet_type.eq."earth") then
 #ifdef CPP_EARTH
-         CALL iniphysiq(ngridmx,llm,daysec,day_ini,dtphys ,
+         CALL iniphysiq(ngridmx,llm,daysec,day_ini,dtphys/nsplit_phys ,
      ,                latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp     )
 #endif