Ignore:
Timestamp:
Nov 28, 2006, 3:40:13 PM (18 years ago)
Author:
lmdzadmin
Message:

Correction bug: le seed pour le simulateur ISCCP doit etre aleatoire.
FC/IM

Location:
LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/calcul_simulISCCP.h

    r644 r738  
    9797      ENDIF !debut
    9898c
    99 cIM: initialisation de seed
     99cIM: initialisation de seed
     100c       DO i=1, klon
     101c         seed(i)=i+100
     102c       ENDDO
     103c
    100104        DO i=1, klon
    101           seed(i)=i+100
     105c
     106         aa=ABS(paprs(i,2)-NINT(paprs(i,2)))
     107         seed_re(i)=1000.*aa+1.
     108         seed(i)=NINT(seed_re(i))
     109c
     110         IF(seed(i).LT.50) THEN
     111c          print*,'seed<50 avant i seed itap paprs',i,
     112c    .     seed(i),itap,paprs(i,2)
     113           seed(i)=50+seed(i)+i+itap
     114           seed_old(i)=seed(i)
     115c
     116           IF(itap.GT.1) then
     117            IF(seed(i).EQ.seed_old(i)) THEN
     118             seed(i)=seed(i)+10
     119             seed_old(i)=seed(i)
     120            ENDIF   
     121           ENDIF
     122c
     123c          print*,'seed<50 apres i seed itap paprs',i,
     124c    .     seed(i),itap,paprs(i,2)
     125c
     126          ELSE IF(seed(i).EQ.0) THEN
     127           print*,'seed=0 i paprs aa seed_re',
     128     .     i,paprs(i,2),aa,seed_re(i)
     129           STOP   
     130          ELSE IF(seed(i).LT.0) THEN
     131           print*,'seed < 0, i seed itap paprs',i,
     132     .     seed(i),itap,paprs(i,2)
     133           STOP   
     134          ENDIF
     135c
    102136        ENDDO
    103137c     
  • LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/physiq.F

    r723 r738  
    381381      REAL nbsunlit(nregISCtot,klon)  !nbsunlit : moyenne de sunlit
    382382      INTEGER ncol, seed(klon)
     383cIM ajout seed_re (REAL)et seed_old (integer)
     384      REAL seed_re(klon), aa
     385      INTEGER seed_old(klon)
     386      SAVE seed_old
     387
    383388
    384389c ncol = nb. de sous-colonnes pour chaque maille du GCM
Note: See TracChangeset for help on using the changeset viewer.