Changeset 2110 for LMDZ5/trunk


Ignore:
Timestamp:
Aug 27, 2014, 5:54:44 PM (10 years ago)
Author:
lguez
Message:

Abort if surface pressure becomes negative. The call to abort_gcm
already was in the sequential version but not in dyn3dpar nor in
dyn3dmem. In dyn3dmem, there is a variable checksum_all, which is
always true. The call to MPI_ALLREDUCE which should update
checksum_all is commented out. I do not know why (performance?).

Non-ASCII characters in comments are not always rendered properly and
they risk being lost. See revision [1740].

Bug fix in procedure convect2: the dimension len must be declared
before the array idcum which has this dimension. Bug fix in procedure
icefrac: the dimensions nl and len must be declared before the arrays
which have these dimensions.

Location:
LMDZ5/trunk/libf
Files:
7 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/dyn3dmem/integrd_loc.F

    r2094 r2110  
    147147!     &                   MPI_LOGICAL,MPI_LOR,COMM_LMDZ,ierr)
    148148      IF( .NOT. checksum ) THEN
    149           write(lunout,*) "integrd: negative surface pressure ",
    150      &                                                ps(stop_it)
     149         write(lunout,*) "integrd: ps = ", ps(stop_it)
    151150         write(lunout,*) " at node ij =", stop_it
    152151         ! since ij=j+(i-1)*jjp1 , we have
     
    155154!         write(lunout,*) " lon = ",rlonv(i)*180./pi, " deg",
    156155!     &                   " lat = ",rlatu(j)*180./pi, " deg"
     156         call abort_gcm("integrd_loc", "negative surface pressure", 1)
    157157      ENDIF
    158158
  • LMDZ5/trunk/libf/dyn3dpar/integrd_p.F

    r1907 r2110  
    137137       
    138138        IF( .NOT. checksum ) THEN
    139          write(lunout,*) "integrd: negative surface pressure ",
    140      &                                                ps(stop_it)
     139           write(lunout,*) "integrd: ps = ", ps(stop_it)
    141140         write(lunout,*) " at node ij =", stop_it
    142141         ! since ij=j+(i-1)*jjp1 , we have
     
    145144         write(lunout,*) " lon = ",rlonv(i)*180./pi, " deg",
    146145     &                   " lat = ",rlatu(j)*180./pi, " deg"
     146         call abort_gcm("integrd_p", "negative surface pressure", 1)
    147147        ENDIF
    148148
  • LMDZ5/trunk/libf/phylmd/clesphys.h

    r2007 r2110  
    44!
    55!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
    6 !                 veillez à n'utiliser que des ! pour les commentaires
    7 !                 et à bien positionner les & des lignes de continuation
     6!                 veillez \`a n'utiliser que des ! pour les commentaires
     7!                 et \`a bien positionner les & des lignes de continuation
    88!                 (les placer en colonne 6 et en colonne 73)
    99!
  • LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

    r2075 r2110  
    16821682    !Config  Desc = activation de la version strato
    16831683    !Config  Def  = .FALSE.
    1684     !Config  Help = active la version stratosphérique de LMDZ de F. Lott
     1684    !Config  Help = active la version stratosph\'erique de LMDZ de F. Lott
    16851685
    16861686    ok_strato_omp=.FALSE.
  • LMDZ5/trunk/libf/phylmd/convect2.F90

    r1992 r2110  
    160160
    161161  INTEGER ncum
     162  INTEGER len
    162163  INTEGER idcum(len)
    163   INTEGER len
    164164  INTEGER nd
    165165  INTEGER ndp1
  • LMDZ5/trunk/libf/phylmd/cv3_routines.F90

    r2078 r2110  
    994994! Calcul de la quantité d'eau sous forme de glace
    995995! --------------------------------------------------------------------
     996  INTEGER nl, len
    996997  REAL qi(len, nl)
    997998  REAL t(len, nl), clw(len, nl)
    998999  REAL fracg
    999   INTEGER nl, len, k, i
     1000  INTEGER k, i
    10001001
    10011002  DO k = 3, nl
  • LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

    r2097 r2110  
    5959    include "iniprint.h"
    6060
    61     ! ug Nouveaux arguments nécessaires au histwrite_mod:
     61    ! ug Nouveaux arguments n\'ecessaires au histwrite_mod:
    6262    INTEGER, INTENT(IN)                         :: ivap
    6363    INTEGER, DIMENSION(klon), INTENT(IN)        :: lmax_th
     
    130130                                                                    90., 90., 90., 90., 90. /)                       
    131131#ifdef CPP_XIOS
    132     ! ug Variables utilisées pour récupérer le calendrier pour xios
     132    ! ug Variables utilis\'ees pour r\'ecup\'erer le calendrier pour xios
    133133    INTEGER :: x_an, x_mois, x_jour
    134134    REAL :: x_heure
     
    261261
    262262#ifdef CPP_XIOS
    263     ! ug Réglage du calendrier xios
     263    ! ug R\'eglage du calendrier xios
    264264    !Temps julian => an, mois, jour, heure
    265265    CALL ju2ymds(zjulian, x_an, x_mois, x_jour, x_heure)
Note: See TracChangeset for help on using the changeset viewer.