Changeset 1468


Ignore:
Timestamp:
Aug 19, 2015, 1:45:31 PM (9 years ago)
Author:
slebonnois
Message:

SL: bug corrections

Location:
trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.TITAN/Tools/tmc.F90

    r1454 r1468  
    281281  do ilev=1,altlength
    282282! Need to be consistent with GCM computations
    283 !    if (za(ilon,ilat,ilev,itim).ne.miss_val) then
     283    if (vitu(ilon,ilat,ilev,itim).ne.miss_val) then
    284284     rayon(ilon,ilat,ilev,itim) = a0
    285285!     rayon(ilon,ilat,ilev,itim) = za(ilon,ilat,ilev,itim) + a0
    286286      grav(ilon,ilat,ilev,itim) = g0*a0*a0 &
    287287                 /(rayon(ilon,ilat,ilev,itim)*rayon(ilon,ilat,ilev,itim))
    288 !    else
    289 !     rayon(ilon,ilat,ilev,itim) = miss_val
    290 !      grav(ilon,ilat,ilev,itim) = miss_val
    291 !    endif
     288    else
     289     rayon(ilon,ilat,ilev,itim) = miss_val
     290      grav(ilon,ilat,ilev,itim) = miss_val
     291    endif
    292292  enddo
    293293 enddo

  • trunk/LMDZ.VENUS/Tools/tmc.F90

    r1454 r1468  
    281281  do ilev=1,altlength
    282282! Need to be consistent with GCM computations
    283 !    if (za(ilon,ilat,ilev,itim).ne.miss_val) then
     283    if (vitu(ilon,ilat,ilev,itim).ne.miss_val) then
    284284     rayon(ilon,ilat,ilev,itim) = a0
    285285!     rayon(ilon,ilat,ilev,itim) = za(ilon,ilat,ilev,itim) + a0
    286286      grav(ilon,ilat,ilev,itim) = g0*a0*a0 &
    287287                 /(rayon(ilon,ilat,ilev,itim)*rayon(ilon,ilat,ilev,itim))
    288 !    else
    289 !     rayon(ilon,ilat,ilev,itim) = miss_val
    290 !      grav(ilon,ilat,ilev,itim) = miss_val
    291 !    endif
     288    else
     289     rayon(ilon,ilat,ilev,itim) = miss_val
     290      grav(ilon,ilat,ilev,itim) = miss_val
     291    endif
    292292  enddo
    293293 enddo

  • trunk/LMDZ.VENUS/libf/phyvenus/chemparam_mod.F90

    r1442 r1468  
    694694                i_h2=i
    695695                PRINT*,'h2',i_h2
    696                 M_tr(i_h2)= 18.0153
     696                M_tr(i_h2)= 2.01588
    697697                CASE('h2o')
    698698                i_h2o=i

  • trunk/LMDZ.VENUS/libf/phyvenus/ini_histins.h

    r1453 r1468  
    103103         CALL histdef(nid_ins, "mmean", "Mean molecular mass", "g/mol",
    104104     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    105      .                "ave(X)", zsto,zout)
     105     .                "ins(X)", zsto,zout)
    106106
    107107         CALL histdef(nid_ins, "rho", "Air density [mass/Vol]", "kg/m3",
    108108     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    109      .                "ave(X)", zsto,zout)
     109     .                "ins(X)", zsto,zout)
    110110
    111111
    … …  
    128128          CALL histdef(nid_ins, "d_qmoldif CO2", "Dif molec" , "kg/kg",
    129129     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    130      .                "ave(X)", zsto,zout)
     130     .                "ins(X)", zsto,zout)
    131131          CALL histdef(nid_ins, "d_qmoldif O3p", "Dif molec" , "kg/kg",
    132132     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    133      .                "ave(X)", zsto,zout)
     133     .                "ins(X)", zsto,zout)
    134134          CALL histdef(nid_ins, "d_qmoldif N2", "Dif molec" , "kg/kg",
    135135     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
    136      .                "ave(X)", zsto,zout)
     136     .                "ins(X)", zsto,zout)
    137137       ENDIF
    138138c
Note: See TracChangeset for help on using the changeset viewer.