Ignore:
Timestamp:
Jul 24, 2024, 4:23:34 PM (6 months ago)
Author:
abarral
Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

Location:
LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat
Files:
13 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.f90

    r5116 r5117  
    299299  !   convergence dynamique pour les traceurs "EAU"
    300300  ! Earth-specific treatment of first 2 tracers (water)
    301    if (planet_type=="earth") THEN
     301   IF (planet_type=="earth") THEN
    302302    DO iq=1,2
    303303     DO l=1,llm
    … …  
    479479
    480480     jH_cur_split=jH_cur+(isplit-1) * dtvr / (daysec *nsplit_phys)
    481      debut_split=debut.and.isplit==1
    482      lafin_split=lafin.and.isplit==nsplit_phys
     481     debut_split=debut.AND.isplit==1
     482     lafin_split=lafin.AND.isplit==nsplit_phys
    483483
    484484   ! if (planet_type=="earth") THEN
    … …  
    493493          zdufi,zdvfi,zdtfi,zdqfi,zdpsrf)
    494494
    495    ! else if ( planet_type=="generic" ) THEN
     495   ! ELSE IF ( planet_type=="generic" ) THEN
    496496   !    CALL physiq (ngridmx,     !! ngrid
    497497  ! .             llm,            !! nlayer
    … …  
    518518  ! .             tracerdyn)      !! tracerdyn <-- utilite ???
    519519
    520   !  endif ! of if (planet_type=="earth")
     520  !  ENDIF ! of if (planet_type=="earth")
    521521
    522522     zufi(:,:)=zufi(:,:)+zdufi(:,:)*zdt_split

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/grid_atob_m.f90

    r5113 r5117  
    244244  WHERE(cham2tmp<0.001) cham2tmp=0.001
    245245  CALL fine2coarse(xtmp,ytmp,xmod,ymod,out,REAL(LOG(cham2tmp)))
    246   out=EXP(out)
    247   amin=MINVAL(out); amax=MAXVAL(out)
     246  out=EXP(OUT)
     247  amin=MINVAL(OUT); amax=MAXVAL(OUT)
    248248  PRINT*, 'Ecart-type du modele:', amin, amax
    249249  out=out/amax*20.0
    250   amin=MINVAL(out); amax=MAXVAL(out)
     250  amin=MINVAL(OUT); amax=MAXVAL(OUT)
    251251  PRINT*, 'Longueur de rugosite du modele:', amin, amax
    252   rugs=REAL(out)
     252  rugs=REAL(OUT)
    253253
    254254END SUBROUTINE rugsoro

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/inigeomphy_mod.F90

    r5116 r5117  
    2323                                 south_west, south_east
    2424  USE mod_interface_dyn_phys, ONLY:  init_interface_dyn_phys
    25   USE nrtype, ONLY: pi
     25  USE lmdz_physical_constants, ONLY: pi
    2626  USE comvert_mod, ONLY: preff, ap, bp, aps, bps, presnivs, &
    2727                         scaleheight, pseudoalt, presinter
    … …  
    144144        boundslonfi_glo(k,north_east)=rlonu(i)
    145145        boundslatfi_glo(k,north_east)=rlatv(j-1)
    146         if (i==1) THEN
     146        IF (i==1) THEN
    147147          ! special case for the first longitude's west bound
    148148          boundslonfi_glo(k,north_west)=rlonu(iim)-2*PI

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_calfis_loc.F90

    r5116 r5117  
    4040    USE parallel_lmdz,ONLY:omp_chunk,using_mpi,jjb_u,jje_u,jjb_v,jje_v &
    4141          ,jj_begin_dyn=>jj_begin,jj_end_dyn=>jj_end
    42     USE Write_Field
     42    USE lmdz_write_field
    4343    Use Write_field_p
    4444    USE Times
    … …  
    188188  ! J'ai été surpris au début que les tableaux zufi_omp, zdufi_omp n'co soitent
    189189  ! en SAVE. Je crois comprendre que c'est parce que tu voulais qu'il
    190   ! soit allocatable (plutot par exemple que de passer une dimension
     190  ! soit ALLOCATABLE (plutot par exemple que de passer une dimension
    191191  ! dépendant du process en argument des routines) et que, du coup,
    192192  ! le SAVE évite d'avoir à refaire l'allocation à chaque appel.
    … …  
    222222!$OMP THREADPRIVATE(firstcal,debut)
    223223
    224   REAL,SAVE,dimension(1:iim,1:llm):: du_send,du_recv,dv_send,dv_recv
     224  REAL,SAVE,DIMENSION(1:iim,1:llm):: du_send,du_recv,dv_send,dv_recv
    225225  INTEGER :: ierr
    226   INTEGER,dimension(MPI_STATUS_SIZE,4) :: Status
    227   INTEGER, dimension(4) :: Req
     226  INTEGER,DIMENSION(MPI_STATUS_SIZE,4) :: Status
     227  INTEGER, DIMENSION(4) :: Req
    228228  REAL,ALLOCATABLE,SAVE:: zdufi2(:,:),zdvfi2(:,:)
    229229  INTEGER :: k,kstart,kend
    … …  
    412412  kend=klon
    413413
    414   if (is_north_pole_dyn) kstart=2
    415   if (is_south_pole_dyn) kend=klon-1
     414  IF (is_north_pole_dyn) kstart=2
     415  IF (is_south_pole_dyn) kend=klon-1
    416416
    417417!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    … …  
    423423      i=index_i(ig0)
    424424      j=index_j(ig0)
    425       if (i==1) THEN
     425      IF (i==1) THEN
    426426        zufi(ig0,l)= 0.5 *(  pucov(iim,j,l)/cu(iim,j) &
    427427              + pucov(1,j,l)/cu(1,j) )
    … …  
    440440  jjb=jj_begin_dyn-1
    441441  jje=jj_end_dyn+1
    442   if (is_north_pole_dyn) jjb=1
    443   if (is_south_pole_dyn) jje=jjm
     442  IF (is_north_pole_dyn) jjb=1
     443  IF (is_south_pole_dyn) jje=jjm
    444444
    445445!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    … …  
    469469      zvfi(ig0,l)= 0.5 *(  pvcov(i,j-1,l)/cv(i,j-1) &
    470470            + pvcov(i,j,l)/cv(i,j) )
    471       if (j==1 .OR. j==jjp1) then !  AdlC MAY 2014
     471      IF (j==1 .OR. j==jjp1) then !  AdlC MAY 2014
    472472        zrfi(ig0,l) = 0 !  AdlC MAY 2014
    473473      else
    … …  
    491491     ! V = 1 / pi  *  integrale [ v * sin(long) * d long ]
    492492
    493   if (is_north_pole_dyn) THEN
     493  IF (is_north_pole_dyn) THEN
    494494!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    495495    DO l=1,llm
    … …  
    511511    ENDDO
    512512!$OMP END DO NOWAIT
    513   endif
     513  ENDIF
    514514
    515515
    … …  
    519519     ! V = 1 / pi  *  integrale [ v * sin(long) * d long ]
    520520
    521   if (is_south_pole_dyn) THEN
     521  IF (is_south_pole_dyn) THEN
    522522!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    523523    DO l=1,llm
    … …  
    538538    ENDDO
    539539!$OMP END DO NOWAIT
    540   endif
     540  ENDIF
    541541
    542542  ! On change de grille, dynamique vers physiq, pour le flux de masse verticale
    … …  
    561561
    562562!$OMP BARRIER
    563   if (first_omp) THEN
     563  IF (first_omp) THEN
    564564    klon=klon_omp
    565565
    … …  
    586586    allocate(flxwfi_omp(klon,llm))
    587587    first_omp=.FALSE.
    588   endif
     588  ENDIF
    589589
    590590
    … …  
    708708
    709709     jH_cur_split=jH_cur+(isplit-1) * dtvr / (daysec *nsplit_phys)
    710      debut_split=debut.and.isplit==1
    711      lafin_split=lafin.and.isplit==nsplit_phys
     710     debut_split=debut.AND.isplit==1
     711     lafin_split=lafin.AND.isplit==nsplit_phys
    712712
    713713    CALL call_physiq(klon,llm,nqtot,tracers(:)%name, &
    … …  
    837837  IF (using_mpi) THEN
    838838
    839   if (MPI_rank>0) THEN
     839  IF (MPI_rank>0) THEN
    840840!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    841841   DO l=1,llm
    … …  
    858858!$OMP BARRIER
    859859
    860   endif
    861 
    862   if (MPI_rank<MPI_Size-1) THEN
     860  ENDIF
     861
     862  IF (MPI_rank<MPI_Size-1) THEN
    863863!$OMP BARRIER
    864864
    … …  
    872872!$OMP END MASTER
    873873
    874   endif
     874  ENDIF
    875875
    876876!$OMP BARRIER
    … …  
    879879!$OMP MASTER
    880880!$OMP CRITICAL (MPI)
    881   if (MPI_rank>0 .and. MPI_rank< MPI_Size-1) THEN
     881  IF (MPI_rank>0 .AND. MPI_rank< MPI_Size-1) THEN
    882882    CALL MPI_WAITALL(4,Req(1),Status,ierr)
    883   else if (MPI_rank>0) THEN
     883  ELSE IF (MPI_rank>0) THEN
    884884    CALL MPI_WAITALL(2,Req(1),Status,ierr)
    885   else if (MPI_rank <MPI_Size-1) THEN
     885  ELSE IF (MPI_rank <MPI_Size-1) THEN
    886886    CALL MPI_WAITALL(2,Req(3),Status,ierr)
    887   endif
     887  ENDIF
    888888!$OMP END CRITICAL (MPI)
    889889!$OMP END MASTER
    … …  
    908908    pdvfi(:,jj_begin,l)=0
    909909
    910     if (.not. is_south_pole_dyn) THEN
     910    IF (.NOT. is_south_pole_dyn) THEN
    911911      pdhfi(:,jj_end:jj_end+1,l)=0
    912912      pdqfi(:,jj_end:jj_end+1,l,:)=0
    … …  
    921921    pdpsfi(:,jj_begin)=0
    922922
    923    if (.not. is_south_pole_dyn) THEN
     923   IF (.NOT. is_south_pole_dyn) THEN
    924924     pdpsfi(:,jj_end:jj_end+1)=0
    925925   endif
    … …  
    937937  kend=klon
    938938
    939   if (is_north_pole_dyn) kstart=2
    940   if (is_south_pole_dyn)  kend=klon-1
     939  IF (is_north_pole_dyn) kstart=2
     940  IF (is_south_pole_dyn)  kend=klon-1
    941941
    942942  !CDIR ON_ADB(index_i)
    … …  
    947947      j=index_j(ig0)
    948948      pdpsfi(i,j) = zdpsrf(ig0)
    949       if (i==1) pdpsfi(iip1,j) =  zdpsrf(ig0)
     949      IF (i==1) pdpsfi(iip1,j) =  zdpsrf(ig0)
    950950     enddo
    951951
    952     if (is_north_pole_dyn) THEN
     952    IF (is_north_pole_dyn) THEN
    953953        DO i=1,iip1
    954954          pdpsfi(i,1)    = zdpsrf(1)
    … …  
    956956    endif
    957957
    958     if (is_south_pole_dyn) THEN
     958    IF (is_south_pole_dyn) THEN
    959959        DO i=1,iip1
    960960          pdpsfi(i,jjp1) = zdpsrf(klon)
    … …  
    970970  kend=klon
    971971
    972   if (is_north_pole_dyn) kstart=2
    973   if (is_south_pole_dyn)  kend=klon-1
     972  IF (is_north_pole_dyn) kstart=2
     973  IF (is_south_pole_dyn)  kend=klon-1
    974974
    975975!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    … …  
    983983      j=index_j(ig0)
    984984      pdhfi(i,j,l) = cpp * zdtfi(ig0,l) / ppk(i,j,l)
    985       if (i==1) pdhfi(iip1,j,l) =  cpp * zdtfi(ig0,l) / ppk(i,j,l)
     985      IF (i==1) pdhfi(iip1,j,l) =  cpp * zdtfi(ig0,l) / ppk(i,j,l)
    986986     enddo
    987987
    988     if (is_north_pole_dyn) THEN
     988    IF (is_north_pole_dyn) THEN
    989989        DO i=1,iip1
    990990          pdhfi(i,1,l)    = cpp *  zdtfi(1,l)      / ppk(i, 1  ,l)
    … …  
    992992    endif
    993993
    994     if (is_south_pole_dyn) THEN
     994    IF (is_south_pole_dyn) THEN
    995995        DO i=1,iip1
    996996          pdhfi(i,jjp1,l) = cpp *  zdtfi(klon,l)/ ppk(i,jjp1,l)
    … …  
    10531053          j=index_j(ig0)
    10541054          pdqfi(i,j,l,iq) = zdqfi(ig0,l,itr)
    1055           if (i==1) pdqfi(iip1,j,l,iq) = zdqfi(ig0,l,itr)
     1055          IF (i==1) pdqfi(iip1,j,l,iq) = zdqfi(ig0,l,itr)
    10561056        ENDDO
    10571057
    … …  
    10831083       j=index_j(ig0)
    10841084
    1085        if (i/=iim) THEN
     1085       IF (i/=iim) THEN
    10861086         pdufi(i,j,l)=0.5*(zdufi2(ig0,l)+zdufi2(ig0+1,l))*cu(i,j)
    10871087       endif
    10881088
    1089        if (i==1) THEN
     1089       IF (i==1) THEN
    10901090          pdufi(iim,j,l)=0.5*(  zdufi2(ig0,l) &
    10911091                + zdufi2(ig0+iim-1,l))*cu(iim,j)
    … …  
    10951095     enddo
    10961096
    1097      if (is_north_pole_dyn) THEN
     1097     IF (is_north_pole_dyn) THEN
    10981098       DO i=1,iip1
    10991099        pdufi(i,1,l)    = 0.
    … …  
    11011101     endif
    11021102
    1103      if (is_south_pole_dyn) THEN
     1103     IF (is_south_pole_dyn) THEN
    11041104       DO i=1,iip1
    11051105        pdufi(i,jjp1,l) = 0.
    … …  
    11161116  kend=klon
    11171117
    1118   if (is_north_pole_dyn) kstart=2
    1119   if (is_south_pole_dyn)  kend=klon-1-iim
     1118  IF (is_north_pole_dyn) kstart=2
     1119  IF (is_south_pole_dyn)  kend=klon-1-iim
    11201120
    11211121!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    … …  
    11281128       j=index_j(ig0)
    11291129       pdvfi(i,j,l)=0.5*(zdvfi2(ig0,l)+zdvfi2(ig0+iim,l))*cv(i,j)
    1130        if (i==1) pdvfi(iip1,j,l) = 0.5*(zdvfi2(ig0,l)+ &
     1130       IF (i==1) pdvfi(iip1,j,l) = 0.5*(zdvfi2(ig0,l)+ &
    11311131             zdvfi2(ig0+iim,l)) &
    11321132             *cv(i,j)
    … …  
    11411141   ! v = U * cos(long) + V * SIN(long)
    11421142
    1143   if (is_north_pole_dyn) THEN
     1143  IF (is_north_pole_dyn) THEN
    11441144!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    11451145    DO l=1,llm
    … …  
    11581158!$OMP END DO NOWAIT
    11591159
    1160   endif
    1161 
    1162   if (is_south_pole_dyn) THEN
     1160  ENDIF
     1161
     1162  IF (is_south_pole_dyn) THEN
    11631163!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    11641164     DO l=1,llm
    … …  
    11771177!$OMP END DO NOWAIT
    11781178
    1179   endif
     1179  ENDIF
    11801180  !-----------------------------------------------------------------------
    11811181

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/lmdz_gr_fi_dyn_p.F90

    r5066 r5117  
    2727        j=index_j(ig)
    2828        pdyn(i,j,ifield)=pfi(ig,ifield)
    29         if (i==1) pdyn(im,j,ifield)=pdyn(i,j,ifield)
     29        IF (i==1) pdyn(im,j,ifield)=pdyn(i,j,ifield)
    3030      ENDDO
    3131

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/mod_interface_dyn_phys.F90

    r5110 r5117  
    33
    44MODULE mod_interface_dyn_phys
    5   INTEGER,SAVE,dimension(:),allocatable :: index_i
    6   INTEGER,SAVE,dimension(:),allocatable :: index_j
     5  INTEGER,SAVE,DIMENSION(:),ALLOCATABLE :: index_i
     6  INTEGER,SAVE,DIMENSION(:),ALLOCATABLE :: index_j
    77 
    88 

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phydev/iniphysiq_mod.F90

    r5113 r5117  
    1818  USE inifis_mod, ONLY: inifis
    1919  USE phyaqua_mod, ONLY: iniaqua
    20   USE nrtype, ONLY: pi
     20  USE lmdz_physical_constants, ONLY: pi
    2121  IMPLICIT NONE
    2222

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/callphysiq_mod.F90

    r5116 r5117  
    6969
    7070  !$OMP MASTER
    71   if (ok_dyn_xios) THEN
     71  IF (ok_dyn_xios) THEN
    7272     CALL xios_get_current_context(dyn3d_ctx_handle)
    73   endif
     73  ENDIF
    7474  !$OMP END MASTER
    7575
    … …  
    9999! switching back to LMDZDYN context
    100100!$OMP MASTER
    101   if (ok_dyn_xios) THEN
     101  IF (ok_dyn_xios) THEN
    102102     CALL xios_set_current_context(dyn3d_ctx_handle)
    103   endif
     103  ENDIF
    104104!$OMP END MASTER
    105105

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/ce0l.F90

    r5116 r5117  
    7373
    7474#ifdef CPP_PARA
    75   integer ierr
     75  INTEGER ierr
    7676#else
    7777! for iniphysiq in serial mode

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90

    r5116 r5117  
    3939  USE comvert_mod, ONLY: ap, bp, preff, pressure_exner
    4040  USE temps_mod, ONLY: annee_ref, day_ref, itau_dyn, itau_phy, start_time
    41   USE strings_mod, ONLY: strLower
     41  USE lmdz_strings, ONLY: strLower
    4242
    4343  IMPLICIT NONE
    … …  
    7878    USE lmdz_filtreg
    7979    USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
     80    USE lmdz_q_sat, ONLY: q_sat
    8081    IMPLICIT NONE
    8182    !-------------------------------------------------------------------------------

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90

    r5116 r5117  
    539539  do j=2,jmp1-1
    540540     PRINT*,'avant if ',cos(rlatu(j)),coslat0
    541      if (cos(rlatu(j))<coslat0) THEN
     541     IF (cos(rlatu(j))<coslat0) THEN
    542542         ! nb de pts affectes par le filtrage de part et d'autre du pt
    543543         ifiltre=(coslat0/cos(rlatu(j))-1.)/2.
    … …  
    548548         wwf(ifiltre+1)=(coslat0/cos(rlatu(j))-1.)/2.-ifiltre
    549549         do i=1,imp1-1
    550             if (masque(i,j)>0.9) THEN
     550            IF (masque(i,j)>0.9) THEN
    551551               ssz=phis(i,j)
    552552               do ifi=1,ifiltre+1
    553553                  ii=i+ifi
    554                   if (ii>imp1-1) ii=ii-imp1+1
     554                  IF (ii>imp1-1) ii=ii-imp1+1
    555555                  ssz=ssz+wwf(ifi)*phis(ii,j)
    556556                  ii=i-ifi
    557                   if (ii<1) ii=ii+imp1-1
     557                  IF (ii<1) ii=ii+imp1-1
    558558                  ssz=ssz+wwf(ifi)*phis(ii,j)
    559559               enddo

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/limit_netcdf.f90

    r5116 r5117  
    378378      REAL :: chmin, chmax, timeday, al
    379379      INTEGER ierr, idx
    380       integer n_extrap ! number of extrapolated points
    381       logical skip
     380      INTEGER n_extrap ! number of extrapolated points
     381      LOGICAL skip
    382382
    383383      !------------------------------------------------------------------------------
    … …  
    599599            CALL pchfe_95(timeyear, champtime(i, j, :), yder, skip, &
    600600                    arth(0.5, real(ndays_in) / ndays, ndays), champan(i, j, :), ierr)
    601             if (ierr < 0) CALL abort_physic("get_2Dfield", "", 1)
     601            IF (ierr < 0) CALL abort_physic("get_2Dfield", "", 1)
    602602            n_extrap = n_extrap + ierr
    603603          END DO

  • LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/test_disvert_m.F90

    r5116 r5117  
    1313    ! the surface pressure, which sample possible values on Earth.
    1414
    15     use exner_hyb_m, ONLY: exner_hyb
    16     use lmdz_vertical_layers, ONLY: ap,bp,preff
    17     use comconst_mod, ONLY: kappa, cpp
     15    USE exner_hyb_m, ONLY: exner_hyb
     16    USE lmdz_vertical_layers, ONLY: ap,bp,preff
     17    USE comconst_mod, ONLY: kappa, cpp
    1818    USE lmdz_abort_physic, ONLY: abort_physic
    1919
    … …  
    2222
    2323    ! Local:
    24     integer l, i
    25     integer, parameter:: ngrid = 7
    26     real p(ngrid, llm + 1) ! pressure at half-level, in Pa
    27     real pks(ngrid) ! exner function at the surface, in J K-1 kg-1
    28     real pk(ngrid, llm) ! exner function at full level, in J K-1 kg-1
    29     real ps(ngrid) ! surface pressure, in Pa
    30     real p_lay(ngrid, llm) ! pressure at full level, in Pa
    31     real delta_ps ! in Pa
     24    INTEGER l, i
     25    INTEGER, parameter:: ngrid = 7
     26    REAL p(ngrid, llm + 1) ! pressure at half-level, in Pa
     27    REAL pks(ngrid) ! exner function at the surface, in J K-1 kg-1
     28    REAL pk(ngrid, llm) ! exner function at full level, in J K-1 kg-1
     29    REAL ps(ngrid) ! surface pressure, in Pa
     30    REAL p_lay(ngrid, llm) ! pressure at full level, in Pa
     31    REAL delta_ps ! in Pa
    3232
    3333    !---------------------
    … …  
    4242
    4343    ! Are pressure values in the right order?
    44     if (any(p(:, :llm) <= p_lay .or. p_lay <= p(:, 2:))) THEN
     44    IF (any(p(:, :llm) <= p_lay .OR. p_lay <= p(:, 2:))) THEN
    4545       ! List details and stop:
    4646       do l = 1, llm
    4747          do i = 1, ngrid
    48              if (p(i, l) <= p_lay(i, l)) THEN
     48             IF (p(i, l) <= p_lay(i, l)) THEN
    4949                print 1000, "ps = ", ps(i) / 100., "hPa, p(level ",  l, &
    5050                     ") = ", p(i, l) / 100., " hPa <= p(layer ", l, ") = ", &
    5151                     p_lay(i, l) / 100., " hPa"
    5252             end if
    53              if (p_lay(i, l) <= p(i, l + 1)) THEN
     53             IF (p_lay(i, l) <= p(i, l + 1)) THEN
    5454                print 1000, "ps = ", ps(i) / 100., "hPa, p(layer ", l, ") = ", &
    5555                     p_lay(i, l) / 100., " hPa <= p(level ", l + 1, ") = ", &
Note: See TracChangeset for help on using the changeset viewer.