Changeset 808 for trunk/LMDZ.VENUS/libf/phyvenus/physiq.F

Timestamp:

Oct 16, 2012, 12:57:35 PM (12 years ago)

Author:

slebonnois

Message:

SL: Many changes for VENUS (related to newstart) and TITAN (related to clouds). Please read DOC/chantiers/commit_importants.log (cf v808).

File:

• trunk/LMDZ.VENUS/libf/phyvenus/physiq.F (modified) (10 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/physiq.F

@@ -527 +527 @@
          ENDIF
 c
-         IF (dtime*FLOAT(radpas).GT.(RDAY*0.25).AND.cycle_diurne)
+         IF (dtime*REAL(radpas).GT.(RDAY*0.25).AND.cycle_diurne)
      $    THEN 
            WRITE(lunout,*)'Nbre d appels au rayonnement insuffisant'
@@ -718 +718 @@
 
       IF (cycle_diurne) THEN
-        zdtime=dtime*FLOAT(radpas) ! pas de temps du rayonnement (s)
+        zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s)
         CALL zenang(zlongi,gmtime,zdtime,rlatd,rlond,rmu0,fract)
       ELSE
@@ -781 +781 @@
 
 ! ADAPTATION GCM POUR CP(T)
+
       CALL clmain(dtime,itap,
      e            t_seri,u_seri,v_seri,
@@ -951 +952 @@
 c    $             ' (itaprad=',itaprad,'/radpas=',radpas,')'
 
-      dtimerad = dtime*FLOAT(radpas)  ! pas de temps du rayonnement (s)
+      dtimerad = dtime*REAL(radpas)  ! pas de temps du rayonnement (s)
 c     PRINT*,'dtimerad,dtime,radpas',dtimerad,dtime,radpas
             
@@ -958 +959 @@
 c     print*,"sollw avant radlwsw=",sollw(klon/2)
 c     print*,"avant radlwsw"
+
       CALL radlwsw
      e            (dist, rmu0, fract, 
@@ -965 +967 @@
      s             sollwdown,
      s             lwnet, swnet)
+
 c     print*,"apres radlwsw"
-
 c     print*,"radsol apres radlwsw=",radsol(klon/2)
 c     print*,"solsw apres radlwsw=",solsw(klon/2)
@@ -1027 +1029 @@
           zn2(i,k) = max(zn2(i,k),1.e-12)  ! securite
         enddo
+        zn2(i,1) = 1.e-12  ! securite
        enddo
 
@@ -1055 +1058 @@
      s                   d_t_oro, d_u_oro, d_v_oro)
 
+c       print*,"d_u_oro=",d_u_oro(klon/2,:)
 c  ajout des tendances
            t_seri(:,:) = t_seri(:,:) + d_t_oro(:,:)
@@ -1074 +1078 @@
 c ----------------------------OROLIFT
       IF (ok_orolf) THEN
+       print*,"ok_orolf NOT IMPLEMENTED !"
+       stop
 c
 c  selection des points pour lesquels le shema est actif:
@@ -1380 +1386 @@
       ENDIF
       
-
       RETURN
       END