Changeset 684 for LMDZ4/trunk/libf/phylmd/calcul_simulISCCP.h

Timestamp:

Apr 4, 2006, 5:00:40 PM (18 years ago)

Author:

lmdzadmin

Message:

Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def +
Lecture freq. ecriture en nombre de jours
IM

File:

• LMDZ4/trunk/libf/phylmd/calcul_simulISCCP.h (modified) (8 diffs)

LMDZ4/trunk/libf/phylmd/calcul_simulISCCP.h

@@ -2 +2 @@
 c $Header$
 c
-cIM 090704 BEG
-c     nbapp_isccp=48
-c     nbapp_isccp=8
-c     nbapp_isccp=6
-c     nbapp_isccp=4 !CPU < 30min pour  9pdt/jour
-      nbapp_isccp=3 !CPU ??           10pdt/jour
-c     nbapp_isccp=2
-c     nbapp_isccp=1
-      isccppas=NINT(86400./dtime/nbapp_isccp)
-cIM 010904 BEG
-cIM   IF (MOD(itap,isccppas).EQ.0) THEN
-c      PRINT*,'itap,isccppas,xjour',itap,isccppas,xjour
+c on appelle le simulateur ISCCP toutes les 3h
+c et on fait des sorties 1 fois par jour 
+c
+c ATTENTION : le temps de calcul peut augmenter considerablement !
+c =============================================================== c
+      DO n=1, napisccp
+c
+      nbapp_isccp=30 !appel toutes les 15h
+      isccppas=NINT((itap*dtime)/3600.) !Nb. d'heures de la physique
+      freqin_pdt(n)=ifreq_isccp(n)
 c
 cIM initialisation nbsunlit pour calculs simulateur ISCCP pdt la journee
@@ -20 +18 @@
        sunlit(i)=1 
        IF(rmu0(i).EQ.0.) sunlit(i)=0
-       nbsunlit(1,i)=FLOAT(sunlit(i))
+       nbsunlit(1,i,n)=FLOAT(sunlit(i))
       ENDDO
 c
@@ -27 +25 @@
       convfra(:,:)=rnebcon(:,:)
       convliq(:,:)=rnebcon(:,:)*clwcon(:,:)
-cIM Amip2 beg
 c
       CALL newmicro (paprs, pplay,ok_newmicro,
@@ -38 +35 @@
      s            cldtaupi, re, fl)
 c
-cIM Amip2 end
-
-c
 cIM calcul tau, emissivite nuages startiformes
-c
-cIM Amip2 beg
 c
       CALL newmicro (paprs, pplay,ok_newmicro,
@@ -53 +45 @@
      e            bl95_b0, bl95_b1,
      s            cldtaupi, re, fl)
-c
-cIM Amip2 end
 c
       cldtot(:,:)=min(max(cldfra(:,:),rnebcon(:,:)),1.)
@@ -86 +76 @@
       close(99)
 c
-cIM: calcul coordonnees regions pour statistiques distribution 
-cIM: nuages en ftion du regime dynamique pour regions oceaniques
-c
-      IF (ok_regdyn) THEN !histREGDYN
-c
-#include "ini_coord_REGDYN.h"
-c
-      ENDIF !ok_regdyn
-c 
       ENDIF !debut
 c
 cIM: initialisation de seed
+c
         DO i=1, klon
-          seed(i)=i+100
+c
+         aa=ABS(paprs(i,2)-NINT(paprs(i,2)))
+         seed_re(i,n)=1000.*aa+1.
+         seed(i,n)=NINT(seed_re(i,n))
+c
+         IF(seed(i,n).LT.50) THEN
+c          print*,'seed<50 avant i seed itap paprs',i,
+c    .     seed(i,n),itap,paprs(i,2)
+           seed(i,n)=50+seed(i,n)+i+itap
+           seed_old(i,n)=seed(i,n)
+c
+           IF(itap.GT.1) then
+            IF(seed(i,n).EQ.seed_old(i,n)) THEN
+             seed(i,n)=seed(i,n)+10
+             seed_old(i,n)=seed(i,n)
+            ENDIF
+           ENDIF
+c
+c          print*,'seed<50 apres i seed itap paprs',i,
+c    .     seed(i,n),itap,paprs(i,2)
+c
+          ELSE IF(seed(i,n).EQ.0) THEN
+           print*,'seed=0 i paprs aa seed_re',
+     .     i,paprs(i,2),aa,seed_re(i,n)
+           STOP
+          ELSE IF(seed(i,n).LT.0) THEN
+           print*,'seed < 0, i seed itap paprs',i,
+     .     seed(i,n),itap,paprs(i,2)
+           STOP
+          ENDIF
+c
         ENDDO
 c     
@@ -126 +138 @@
      &     debug,
      &     debugcol,
-cIM 300704    &     itap, debut,
-cIM 300604 klon !BAD
      &     klon,
      &     sunlit,
      &     klev,
-     &     ncol,
-     &     seed,
+     &     ncol(n),
+     &     seed(:,n),
      &     pfull,
      &     phalf,
@@ -143 +153 @@
      &     emsfc_lw,
      &     at, dem_sH2B, dem_cH2B,
-     &     fq_isccp,
-     &     totalcldarea,
-     &     meanptop,
-     &     meantaucld,
-     &     boxtau,
-     &     boxptop)
+     &     fq_isccp(:,:,:,n),
+     &     totalcldarea(:,n),
+     &     meanptop(:,n),
+     &     meantaucld(:,n),
+     &     boxtau(:,:,n),
+     &     boxptop(:,:,n))
 c
-c calcul regime dynamique sur les regions fixees
-c
-       IF (ok_regdyn) THEN !histREGDYN
-c
-#include "calcul_REGDYN.h"
-c
-       ENDIF !(ok_regdyn) THEN !histREGDYN
-cIM    ENDIF !(MOD(itaprad,radpas).EQ.0) THEN
-cIM 010904 END
+      ENDDO !n=1, napisccp
+