Changeset 1456 for LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/physiq.F

Timestamp:

Nov 22, 2010, 2:16:44 PM (14 years ago)

Author:

musat

Message:

phyetat0, phyredem: correction dimension verticale pbl_tke: pbl_tke(:,1:klev+1,:)
physiq: pour pouvoir fixer la longitude solaire avec la nouvelle orbite
JYG/IM

File:

• LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/physiq.F (modified) (3 diffs)

LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/physiq.F

@@ -1250 +1250 @@
 cym Attention pbase pas initialise dans concvl !!!!
           pbase=0
+          paire_ter(:)=0.    
 cIM 180608
 c         pmflxr=0.
@@ -1846 +1847 @@
 ! doit donc etre placé avant radlwsw et pbl_surface
 
-! calcul selon la routine utilisee pour les planetes
-      if (new_orbit) then
-        call ymds2ju(year_cur, mth_eq, day_eq,0., jD_eq)
-        day_since_equinox = (jD_cur + jH_cur) - jD_eq 
-!        day_since_equinox = (jD_cur) - jD_eq 
-        call solarlong(day_since_equinox, zlongi, dist)
-      else      
-! calcul selon la routine utilisee pour l'AR4
+!!!   jyg 17 Sep 2010 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+      call ymds2ju(year_cur, mth_eq, day_eq,0., jD_eq)
+      day_since_equinox = (jD_cur + jH_cur) - jD_eq
+!
 !   choix entre calcul de la longitude solaire vraie ou valeur fixee a
 !   solarlong0
-        if (solarlong0<-999.) then
-           CALL orbite(REAL(days_elapsed+1),zlongi,dist)
-        else
+      if (solarlong0<-999.) then
+       if (new_orbit) then
+! calcul selon la routine utilisee pour les planetes
+        call solarlong(day_since_equinox, zlongi, dist)
+       else
+! calcul selon la routine utilisee pour l'AR4
+        CALL orbite(REAL(days_elapsed+1),zlongi,dist)
+       endif
+      else
            zlongi=solarlong0  ! longitude solaire vraie
            dist=1.            ! distance au soleil / moyenne
-        endif
       endif
+!!!
+
       if(prt_level.ge.1)                                                &
      &    write(lunout,*)'Longitude solaire ',zlongi,solarlong0,dist
@@ -3375 +3379 @@
      I                   cdragh,coefh,u1,v1,ftsol,pctsrf,
      I                   frac_impa, frac_nucl,
-     I                   pphis,airephy,dtime,itap,
-     I                   rlon,rlat,qx(:,:,ivap),da,phi,mp,upwd,dnwd)
+     I                   pphis,airephy,dtime,itap)