Ignore:
Timestamp:
Mar 9, 2010, 3:42:32 PM (15 years ago)
Author:
Laurent Fairhead
Message:

The ability to read annual aerosol forcing files that had disappeared in
recent versions of the code is put back in


On ré-intégre la possibilité de lire des fichiers de forçage aérosols de
type annuel

File:
1 edited

Legend:

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

    r1319 r1321  
    6666     cyear='.nat'
    6767     ! get_aero_fromfile returns pt_out allocated and initialized with data for 12 month
     68     ! pt_out has dimensions (klon, klev_src, 12)
     69     CALL get_aero_fromfile(name_aero, cyear, klev_src, pt_ap, pt_b, p0, pt_out, psurf, load)
     70     
     71  ELSE IF (type == 'annuel') THEN
     72! Read and return data from scenario annual files
     73!****************************************************************************************     
     74     WRITE(cyear,'(I4)') iyr_in
     75     WRITE(lunout,*) 'get_aero 3 iyr_in=', iyr_in,'   ',cyear
     76     ! get_aero_fromfile returns pt_out allocated and initialized with data for nbr_tsteps month
    6877     ! pt_out has dimensions (klon, klev_src, 12)
    6978     CALL get_aero_fromfile(name_aero, cyear, klev_src, pt_ap, pt_b, p0, pt_out, psurf, load)
     
    431440
    432441       IF ((pt_b(1) < pt_b(klev_src)) .OR. .NOT. new_file) THEN
    433           WRITE(lunout,*) 'Vertical axis in file ',TRIM(fname), ' needs to be inverted'
    434           WRITE(lunout,*) 'before pt_ap = ', pt_ap
    435           WRITE(lunout,*) 'before pt_b = ', pt_b
     442!          WRITE(lunout,*) 'Vertical axis in file ',TRIM(fname), ' needs to be inverted'
     443!          WRITE(lunout,*) 'before pt_ap = ', pt_ap
     444!          WRITE(lunout,*) 'before pt_b = ', pt_b
    436445         
    437446          ! Inverse vertical levels for varyear
Note: See TracChangeset for help on using the changeset viewer.