Changeset 1529

Timestamp:

May 26, 2011, 5:17:33 PM (13 years ago)

Author:

Laurent Fairhead

Message:

Necessary modifications to use the model in an aquaplanet or terraplanet configuration

---

Modifications nécessaires pour l'utilisation du moèle en configuration aquaplanète
ou terraplanète

Dans la dynamique:
Changement pour le paramètre iflag_phys:

auparavant : 0 pas de physique, 1 phyique, 2 rappel
ici on ajoute iflag_phys>=100 pour les aqua et terra

Du coup, dans leapfrog.F et gcm.F on appelle la physique si
iflag_phys=1.or.iflag_phys>=100
Dans iniacademic, on initialise si iflag_phys>=2 au lieu de =2
Dans gcm.F, on appelle en plus iniaqua (sous une clef CPP_EARTH)
Dans iniacademic, on met ps=108080 pour les aqua et terra pour répondre
à une specification CFMIP.

Dans la physique:
On ajoute phyaqua.F qui contient :
iniaqua : initialise les startphy.nc et limit.nc pour la phyique
zenang_an : calcule un ensolleillement moyen sur l'année en fonction de

la latiude (A. Campoy).

profil_sst : qui calcule différents profils latitudinaux de SSTs.
writelim : écrit le fichier limit.nc

Dans physiq.F
On ajoute la possibilité d'appeller un calcul d'ensoleillement moyen
sur l'année quand solarlong0=1000.

Location:

LMDZ5/trunk/libf

Files:

• dyn3d/gcm.F (modified) (2 diffs)
• dyn3d/iniacademic.F90 (modified) (2 diffs)
• dyn3d/leapfrog.F (modified) (4 diffs)
• phylmd/phyaqua.F (added)
• phylmd/physiq.F (modified) (1 diff)

LMDZ5/trunk/libf/dyn3d/gcm.F

@@ -408 +408 @@
 c   -------------------------------
 
-      IF (call_iniphys.and.(iflag_phys.eq.1)) THEN
+      IF (call_iniphys.and.(iflag_phys==1.or.iflag_phys>=100)) THEN
          latfi(1)=rlatu(1)
          lonfi(1)=0.
@@ -459 +459 @@
 #endif
 
+#ifdef CPP_EARTH
+! Create start file (startphy.nc) and boundary conditions (limit.nc) 
+! for the Earth verstion
+       if (iflag_phys>=100) then
+          call iniaqua(ngridmx,latfi,lonfi,iflag_phys)
+       endif
+#endif
+
 !      if (planet_type.eq."earth") then
 ! Write an Earth-format restart file
+
         CALL dynredem0("restart.nc", day_end, phis)
 !      endif

LMDZ5/trunk/libf/dyn3d/iniacademic.F90

@@ -115 +115 @@
   endif
 
-  academic_case: if (iflag_phys == 2) then
+  academic_case: if (iflag_phys >= 2) then
      ! initializations
 
@@ -208 +208 @@
      IF (.NOT. read_start) THEN
         ! surface pressure
-        ps(:)=preff
+        if (iflag_phys>2) then
+           ps(:)=preff
+        else
+           ps(:)=101080.
+        endif
         ! ground geopotential
         phis(:)=0.

LMDZ5/trunk/libf/dyn3d/leapfrog.F

@@ -149 +149 @@
       logical ok_sync
       parameter (ok_sync = .true.) 
+      logical physic
 
       data callinigrads/.true./
@@ -193 +194 @@
       itaufin   = nday*day_step
       itaufinp1 = itaufin +1
-      
-
       itau = 0
+      physic=.true.
+      if (iflag_phys==0.or.iflag_phys==2) physic=.false.
+
 c      iday = day_ini+itau/day_step
 c      time = REAL(itau-(iday-day_ini)*day_step)/day_step+time_0
@@ -281 +283 @@
      s        apdiss = .TRUE.
          IF( MOD(itau,iphysiq ).EQ.0.AND..NOT.forward 
-     s          .and. iflag_phys.EQ.1                 ) apphys = .TRUE.
+     s          .and. physic                        ) apphys = .TRUE.
       ELSE
       ! Leapfrog/Matsuno time stepping 
@@ -287 +289 @@
          IF( MOD(itau+1,dissip_period).EQ.0 .AND. .NOT. forward )
      s        apdiss = .TRUE.
-         IF( MOD(itau+1,iphysiq).EQ.0.AND.iflag_phys.EQ.1) apphys=.TRUE.
+         IF( MOD(itau+1,iphysiq).EQ.0.AND.physic       ) apphys=.TRUE.
       END IF
 

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -1881 +1881 @@
      &    write(lunout,*)'Longitude solaire ',zlongi,solarlong0,dist
 
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! Calcul de l'ensoleillement :
+! ============================
+! Pour une solarlong0=1000., on calcule un ensoleillement moyen sur
+! l'annee a partir d'une formule analytique.
+! Cet ensoleillement est symmétrique autour de l'équateur et
+! non nul aux poles.
+      IF (abs(solarlong0-1000.)<1.e-4) then
+         call zenang_an(cycle_diurne,jH_cur,rlat,rlon,rmu0,fract)
+      ELSE
 !  Avec ou sans cycle diurne
-      IF (cycle_diurne) THEN
-        zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s)
-        CALL zenang(zlongi,jH_cur,zdtime,rlat,rlon,rmu0,fract)
-      ELSE
-        CALL angle(zlongi, rlat, fract, rmu0)
+         IF (cycle_diurne) THEN
+           zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s)
+           CALL zenang(zlongi,jH_cur,zdtime,rlat,rlon,rmu0,fract)
+         ELSE
+           CALL angle(zlongi, rlat, fract, rmu0)
+         ENDIF
       ENDIF