Changeset 1416


Ignore:
Timestamp:
Jul 16, 2010, 11:29:10 AM (14 years ago)
Author:
jghattas
Message:

Bug corrections for nudged run (pres2lev.F90, guide_p_mod.F90) :

  • now the results are the same for sequentiel and parallel mode(if adjust=n and use_filtre_fft=n).
  • the results are the same as the sequential mode in previous revision.
  • test done only with guide_u=y,guide_v=y

Added condition read_climoz for the variable O3daySTD(calcul_STDlev.h,
undefSTD.F)

ACo, YM, JG

Location:
LMDZ4/branches/LMDZ4_AR5/libf
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4_AR5/libf/dyn3dpar/guide_p_mod.F90

    r1304 r1416  
    390390              stop
    391391          ELSE
    392               IF (guide_v) vnat1(jjb_v:jje_v,:,:)=vnat2(jjb_v:jje_v,:,:)
    393               IF (guide_u) unat1(jjb_u:jje_u,:,:)=unat2(jjb_u:jje_u,:,:)
    394               IF (guide_T) tnat1(jjb_u:jje_u,:,:)=tnat2(jjb_u:jje_u,:,:)
    395               IF (guide_Q) qnat1(jjb_u:jje_u,:,:)=qnat2(jjb_u:jje_u,:,:)
    396               IF (guide_P.OR.guide_modele) psnat1(jjb_u:jje_u,:)=psnat2(jjb_u:jje_u,:)
     392              IF (guide_v) vnat1(:,jjb_v:jje_v,:)=vnat2(:,jjb_v:jje_v,:)
     393              IF (guide_u) unat1(:,jjb_u:jje_u,:)=unat2(:,jjb_u:jje_u,:)
     394              IF (guide_T) tnat1(:,jjb_u:jje_u,:)=tnat2(:,jjb_u:jje_u,:)
     395              IF (guide_Q) qnat1(:,jjb_u:jje_u,:)=qnat2(:,jjb_u:jje_u,:)
     396              IF (guide_P.OR.guide_modele) psnat1(:,jjb_u:jje_u)=psnat2(:,jjb_u:jje_u)
    397397              step_rea=step_rea+1
    398398              itau_test=itau
  • LMDZ4/branches/LMDZ4_AR5/libf/dyn3dpar/pres2lev.F90

    r1345 r1416  
    6060      DO k=1,ni*nj
    6161        IF (pn(k,ln) >= zpo(k,1) ) THEN
    62           varn(k,ln) = varo(k,1)
     62          varn(k,ln) = zvaro(k,1)
    6363        ELSE IF (pn(k,ln) <= zpo(k,lmo)) THEN
    6464          varn(k,ln) = zvaro(k,lmo)
  • LMDZ4/branches/LMDZ4_AR5/libf/phylmd/calcul_STDlev.h

    r1398 r1416  
    128128cIM on somme les valeurs toutes les freq_calNMC secondes
    129129c
    130        CALL undefSTD(itap,freq_calNMC)
     130       CALL undefSTD(itap,freq_calNMC, read_climoz)
    131131c
    132132cIM on moyenne a la fin du mois ou du jour (toutes les freq_outNMC secondes)
  • LMDZ4/branches/LMDZ4_AR5/libf/phylmd/undefSTD.F

    r1398 r1416  
    22! $Header$
    33!
    4       SUBROUTINE undefSTD(itap,freq_calNMC)
     4      SUBROUTINE undefSTD(itap,freq_calNMC, read_climoz)
    55      USE netcdf
    66      USE dimphy
     
    4343      INTEGER i, k, n
    4444      REAL freq_calNMC(nout)
     45      INTEGER read_climoz
    4546c
    4647c variables Output
     
    8586           T2sumSTD(i,k,n)=T2sumSTD(i,k,n)+T2STD(i,k)
    8687           O3sumSTD(i,k,n)=O3sumSTD(i,k,n)+O3STD(i,k)
    87            O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)+O3daySTD(i,k)
     88           IF (read_climoz==2)
     89     &          O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)+O3daySTD(i,k)
     90
    8891          ENDIF
    8992         ENDDO !i
Note: See TracChangeset for help on using the changeset viewer.