Changeset 883 for LMDZ4/trunk/libf/phylmd/physiq.F

Timestamp:

Jan 31, 2008, 6:20:42 PM (16 years ago)

Author:

Laurent Fairhead

Message:

modifications pour faire de l'aquaplanète FH
LF

File:

• LMDZ4/trunk/libf/phylmd/physiq.F (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -687 +687 @@
 c
       REAL qsol(klon)
+      REAL,save :: qsol0
+      REAL,save ::  solarlong0
 c
       REAL,allocatable,save :: falbe(:,:)
@@ -1478 +1480 @@
 c======================================================================
 
-      if (klon.eq.1) then
-          print*,'WARNING !!!! omega=0'
-          omega=0.
-          igout=1
+      if (1.eq.1) then
+          igout=klon/2
          write(lunout,*) 'DEBUT DE PHYSIQ !!!!!!!!!!!!!!!!!!!!'
          write(lunout,*)
@@ -1738 +1738 @@
 c
          call conf_phys(ocean, ok_veget, ok_journe, ok_mensuel,
-     .                  ok_instan, ok_hf, seuil_inversion,
+     .                  ok_instan, ok_hf,
+     .                  solarlong0,qsol0,seuil_inversion,
      .                  fact_cldcon, facttemps,ok_newmicro,iflag_radia,
      .                  iflag_cldcon,iflag_ratqs,ratqsbas,ratqshaut,
@@ -2256 +2257 @@
 
 c
-C calculs necessaires au calcul de l'albedo dans l'interface
-c
-      CALL orbite(FLOAT(julien),zlongi,dist)
+c=========================================================================
+! Calculs de l'orbite.
+! Necessaires pour le rayonnement et la surface (calcul de l'albedo).
+! doit donc etre placé avant radlwsw et pbl_surface
+
+!   choix entre calcul de la longitude solaire vraie ou valeur fixee a
+!   solarlong0
+
+      if (solarlong0<-999.) then
+         CALL orbite(FLOAT(julien),zlongi,dist)
+      else
+         zlongi=solarlong0  ! longitude solaire vraie
+         dist=1.            ! distance au soleil / moyenne
+      endif
+
+      print*,'Longitude solaire ',zlongi,solarlong0
+
+!  Avec ou sans cycle diurne
       IF (cycle_diurne) THEN
         zdtime=dtime*FLOAT(radpas) ! pas de temps du rayonnement (s)
@@ -2292 +2308 @@
 c   dsens, devap, zxsnow, zxfluxt, zxfluxq, q2m, fluxq
 c
+
+      if (qsol0>0.) qsol(:)=qsol0
+
       CALL pbl_surface( 
      e     dtime,     date0,     itap,    julien,