Ignore:
Timestamp:
Jul 12, 2016 4:34:24 PM (5 years ago)
Author:
suehring
Message:

Separate balance equations for humidity and passive_scalar

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/read_3d_binary.f90

    r1852 r1960  
    107107        ONLY:  e, kh, km, ol, p, pt, q, ql, qc, nr, nrs, nrsws, nrswst,        &
    108108               prr, precipitation_amount, qr,                                  &
    109                qrs, qrsws, qrswst, qs, qsws, qswst, sa, saswsb, saswst,        &
    110                rif_wall, shf, ts, tswst, u, u_m_l, u_m_n, u_m_r, u_m_s, us,    &
    111                usws, uswst, v, v_m_l, v_m_n, v_m_r, v_m_s, vpt, vsws, vswst,   &
    112                w, w_m_l, w_m_n, w_m_r, w_m_s, z0, z0h, z0q
     109               qrs, qrsws, qrswst, qs, qsws, qswst, s, sa, saswsb, saswst,     &
     110               ss, ssws, sswst, rif_wall, shf, ss, ts, tswst, u, u_m_l, u_m_n, &
     111               u_m_r, u_m_s, us, usws, uswst, v, v_m_l, v_m_n, v_m_r, v_m_s,   &
     112               vpt, vsws, vswst, w, w_m_l, w_m_n, w_m_r, w_m_s, z0, z0h, z0q
    113113
    114114    USE averaging
     
    985985                   rif_wall(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp,:) = &
    986986                            tmp_4d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp,:)
     987                           
     988                CASE ( 's' )
     989                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     990                   s(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     991                                    tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    987992
    988993                CASE ( 's_av' )
     
    10671072                      ENDIF
    10681073                   ENDIF
    1069 
     1074                   
     1075                CASE ( 'ss' )
     1076                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     1077                   ss(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     1078                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     1079
     1080                CASE ( 'ssws' )
     1081                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     1082                   ssws(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     1083                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     1084
     1085                CASE ( 'ssws_av' )
     1086                   IF ( .NOT. ALLOCATED( ssws_av ) )  THEN
     1087                      ALLOCATE( ssws_av(nysg:nyng,nxlg:nxrg) )
     1088                   ENDIF 
     1089                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     1090                   ssws_av(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     1091                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     1092
     1093                CASE ( 'sswst' )
     1094                   IF ( k == 1 )  READ ( 13 )  tmp_2d
     1095                   sswst(nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp)  = &
     1096                                          tmp_2d(nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
     1097                 
    10701098                CASE ( 'ts' )
    10711099                   IF ( k == 1 )  READ ( 13 )  tmp_2d
Note: See TracChangeset for help on using the changeset viewer.