Changeset 1458 for LMDZ5/trunk

Timestamp:

Nov 22, 2010, 4:32:14 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

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phyetat0.F (modified) (2 diffs)
• phyredem.F (modified) (1 diff)
• physiq.F (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/phyetat0.F

@@ -817 +817 @@
           ENDIF
           WRITE(str2,'(i2.2)') nsrf
-          CALL get_field("TKE"//str2,pbl_tke(:,1:klev,nsrf),found)
+          CALL get_field("TKE"//str2,pbl_tke(:,1:klev+1,nsrf),found)
           IF (.NOT. found) THEN
             PRINT*, "phyetat0: <TKE"//str2//"> est absent"
@@ -824 +824 @@
           xmin = 1.0E+20
           xmax = -1.0E+20
-          DO k = 1, klev
+          DO k = 1, klev+1
             DO i = 1, klon
               xmin = MIN(pbl_tke(i,k,nsrf),xmin)

LMDZ5/trunk/libf/phylmd/phyredem.F

@@ -289 +289 @@
             WRITE(str2,'(i2.2)') nsrf
             CALL put_field("TKE"//str2,"Energ. Cineti. Turb."//str2,
-     .                     pbl_tke(:,1:klev,nsrf))
+     .                     pbl_tke(:,1:klev+1,nsrf))
           ELSE
             PRINT*, "Trop de sous-mailles"

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -1846 +1846 @@
 ! doit donc etre placé avant radlwsw et pbl_surface
 
+!!!   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
+       if (new_orbit) then
 ! 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      
+       else
 ! calcul selon la routine utilisee pour l'AR4
-!   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
+        CALL orbite(REAL(days_elapsed+1),zlongi,dist)
+       endif
+      else
            zlongi=solarlong0  ! longitude solaire vraie
-           dist=1.            ! distance au soleil / moyenne
-        endif
+           dist=1.            ! distance au soleil / moyenne 
       endif
       if(prt_level.ge.1)                                                &