Ignore:
Timestamp:
Jan 24, 2013, 3:00:03 PM (12 years ago)
Author:
aslmd
Message:

LMDZ.GENERIC. Added patches for large domains. Added C2H2 and C2H6 which must be named 2H2 and 2H6 in gases.def

Location:
trunk/LMDZ.GENERIC/libf/phystd
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/libf/phystd/calc_cpp_mugaz.F90

    r716 r868  
    6767               cpp_c   = cpp_c   + 2.175*gfrac(igas)
    6868               mugaz_c = mugaz_c + 17.03*gfrac(igas)
    69                print*,'WARNING, cpp for NH3 may be for liquid'
     69               print*,'WARNING, cpp for NH3 may be for liquid'
     70            elseif(gnom(igas).eq.'2H6')then
     71               ! C2H6 http://encyclopedia.airliquide.com/Encyclopedia.asp?GasID=28
     72               cpp_c   = cpp_c   + 1.763*gfrac(igas)
     73               mugaz_c = mugaz_c + 30.07*gfrac(igas)
     74            elseif(gnom(igas).eq.'2H2')then
     75               ! C2H2 http://encyclopedia.airliquide.com/Encyclopedia.asp?GasID=1
     76               cpp_c   = cpp_c   + 1.575*gfrac(igas)
     77               mugaz_c = mugaz_c + 26.04*gfrac(igas)
    7078            else
    7179               print*,'Error in calc_cpp_mugaz: Gas species not recognised!'

  • trunk/LMDZ.GENERIC/libf/phystd/gases_h.F90

    r716 r868  
    2929      integer :: igas_CH4
    3030      integer :: igas_NH3
     31      integer :: igas_2H2
     32      integer :: igas_2H6
    3133
    3234      end module gases_h

  • trunk/LMDZ.GENERIC/libf/phystd/su_gases.F90

    r716 r868  
    102102           igas_NH3=igas
    103103           count=count+1
     104        elseif (gnom(igas).eq."2H6") then
     105           ! this is ethane
     106           igas_2H6=igas
     107           count=count+1
     108        elseif (gnom(igas).eq."2H2") then
     109           ! this is acetylene
     110           igas_2H2=igas
     111           count=count+1
    104112        endif
    105113     enddo

  • trunk/LMDZ.GENERIC/libf/phystd/sugas_corrk.F90

    r810 r868  
    124124                 'match that in gases.def (',gnom(igas),'), exiting. You should compare ', &
    125125                 'gases.def with Q.dat in your radiative transfer directory.'
    126             call abort
     126            !call abort
    127127         endif
    128128      enddo
Note: See TracChangeset for help on using the changeset viewer.