Changeset 3834


Ignore:
Timestamp:
Jul 7, 2025, 1:39:14 PM (3 days ago)
Author:
afalco
Message:

Pluto: fix for zrecast, reading phisfi when phisinit unavailable.
AF

Location:
trunk/LMDZ.PLUTO/util
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.PLUTO/util/compile

    r3815 r3834  
    1111
    1212FLAGS=
    13 # FLAGS=-g -O0
     13# FLAGS="-g -O0"
    1414
    1515#pgf90 $SOURCE.F90 \
  • trunk/LMDZ.PLUTO/util/zrecast.F90

    r3815 r3834  
    596596    write(*,*) 'OK, got phisinit or phsifi'
    597597    ierr=NF_CLOSE(infid2)
     598  else
     599    ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
     600    if (ierr.ne.NF_NOERR) then
     601        stop "Error: Failed reading phisinit/phisfi"
     602    endif
    598603  endif
    599604else
     
    603608  endif
    604609endif
     610print*,"phisinit range:", minval(phisinit), maxval(phisinit)
    605611
    606612!===============================================================================
     
    953959    call crude_gcm_za(lonlength,latlength,altlength,timelength, &
    954960                      phisinit,ps,press,temp,za_gcm)
     961    print*, "GCM altitude range:", minval(za_gcm),maxval(za_gcm)
    955962  else if (ztype.eq.3) then ! above local surface altitudes
    956963    allocate(zs_gcm(lonlength,latlength,altlength,timelength),stat=ierr)
Note: See TracChangeset for help on using the changeset viewer.