Ignore:
Timestamp:
Jul 19, 2010, 5:11:24 PM (14 years ago)
Author:
jghattas
Message:

Following bug corrections are already done at branch LMDZ4_AR5 rev
1417,1416 :

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 compiled with vsafe at mercure).
  • the results are the same as the sequential mode in previous revision.
  • test done only with guide_u=y,guide_v=y
  • copied optimized version of pres2lev.F90 from dyn3dpar to dyn3d


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

Corrected bug in calculation of the diagnostic variable ec550aer
"Extinction at 550nm" (aeropt_5wv.F90) (Maria Raffaella Vuolo, LSCE)

Added stop if use_filtre_fft=y with dyn3d (conf_gcm.F) : this option is
not implemented in dyn3d.

Location:
LMDZ4/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/aeropt_5wv.F90

    r1337 r1418  
    850850!  ENDDO
    851851
    852    DO i=1, klon
    853       od550aer(i)=SUM(tausum(i,2,:))
    854       od865aer(i)=SUM(tausum(i,5,:))
    855       DO k=1, KLEV
    856          ec550aer(i,k)=SUM(tau(i,k,2,:))/dh(i,k)
    857       END DO   
    858    END DO
     852  DO i=1, klon
     853     od550aer(i)=0.
     854     DO m=1,naero_spc
     855        od550aer(i)=od550aer(i)+tausum(i,2,m)
     856     END DO
     857  END DO
     858  DO i=1, klon
     859     od865aer(i)=0.
     860     DO m=1,naero_spc
     861        od865aer(i)=od865aer(i)+tausum(i,5,m)
     862     END DO
     863  END DO
     864  DO i=1, klon
     865     DO k=1, KLEV
     866        ec550aer(i,k)=0.
     867        DO m=1,naero_spc
     868           ec550aer(i,k)=ec550aer(i,k)+tau(i,k,2,m)/dh(i,k)
     869        END DO
     870     END DO
     871  END DO
     872 
    859873   od550lt1aer(:)=tausum(:,2,id_ASSO4M)+tausum(:,2,id_ASBCM)+tausum(:,2,id_AIBCM)+ &
    860874        tausum(:,2,id_ASPOMM)+tausum(:,2,id_AIPOMM)+tausum(:,2,id_ASSSM)+ &
  • LMDZ4/trunk/libf/phylmd/calcul_STDlev.h

    r1398 r1418  
    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/trunk/libf/phylmd/undefSTD.F

    r1403 r1418  
    22! $Id$
    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.