Ignore:
Timestamp:
05/26/22 20:21:20 (2 years ago)
Author:
oboucher
Message:

Correct some errors in the writing of ss_w and ss_g

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TOOLS/CMIP6_FORCING/AER_OPTICS/mie_6bands_nitrate_v4.f

    r3385 r6150  
    145145c 
    146146       OPEN (unit=14, file='SEXT_nitrate_soluble_6bands.txt') 
    147        OPEN (unit=15, file='G_nitrate_soluble_6bands.txt') 
    148        OPEN (unit=16, file='SSA_nitrate_soluble_6bands.txt') 
     147       OPEN (unit=15, file='SSA_nitrate_soluble_6bands.txt') 
     148       OPEN (unit=16, file='G_nitrate_soluble_6bands.txt') 
    149149c 
    150150       OPEN (unit=34, file='SEXT_nitrate_soluble_5wave.txt') 
     
    154154c 
    155155       OPEN (unit=24, file='SEXT_nitrate_insoluble_6bands.txt') 
    156        OPEN (unit=25, file='G_nitrate_insoluble_6bands.txt') 
    157        OPEN (unit=26, file='SSA_nitrate_insoluble_6bands.txt') 
     156       OPEN (unit=25, file='SSA_nitrate_insoluble_6bands.txt') 
     157       OPEN (unit=26, file='G_nitrate_insoluble_6bands.txt') 
    158158c 
    159159c--initializing wavelengths 
     
    410410        WRITE(15,*) '  ! '//chmode(mode) 
    411411        DO k=1, NbandSW 
    412           WRITE(15,951) (ss_g(k,IRH),IRH=1,rh_int) 
     412          WRITE(15,951) (ss_w(k,IRH),IRH=1,rh_int) 
    413413        ENDDO 
    414414        WRITE(16,*) '  ! '//chmode(mode) 
    415415        DO k=1, NbandSW 
    416           WRITE(16,951) (ss_w(k,IRH),IRH=1,rh_int) 
     416          WRITE(16,951) (ss_g(k,IRH),IRH=1,rh_int) 
    417417        ENDDO 
    418418c 
     
    434434        WRITE(24,953) (ss_a(k,1),k=1,NbandSW) 
    435435        WRITE(25,*) '  ! Nitrate insoluble' 
    436         WRITE(25,953) (ss_g(k,1),k=1,NbandSW) 
     436        WRITE(25,953) (ss_w(k,1),k=1,NbandSW) 
    437437        WRITE(26,*) '  ! Nitrate insoluble' 
    438         WRITE(26,953) (ss_w(k,1),k=1,NbandSW) 
     438        WRITE(26,953) (ss_g(k,1),k=1,NbandSW) 
    439439c 
    440440        WRITE(35,*) '  ! extinction Nitrate insoluble' 
Note: See TracChangeset for help on using the changeset viewer.