Changeset 1238 for trunk/LMDZ.COMMON


Ignore:
Timestamp:
May 7, 2014, 12:35:26 PM (11 years ago)
Author:
emillour
Message:

Mars GCM and common dynamics:

Common dynamics:

  • correction in inidissip (only matters in Martian case)
  • added correction in addfi on theta to account for surface pressure change.

Mars GCM:
Some fixes and updates:

  • addfi (dyn3d): Add correction on theta when surface pressure changes
  • vdif_cd (phymars): Correction for coefficients in stable nighttime case
  • jthermcalc (aeronomars): Fix for some pathological cases (further investigations on the origin of these is needed)

EM

Location:
trunk/LMDZ.COMMON/libf
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.COMMON/libf/dyn3d/addfi.F

    r1189 r1238  
    103103       ENDDO
    104104      ENDDO
    105 c
     105!***********************
     106! Correction on teta due to surface pressure changes
     107      DO k = 1,llm
     108        DO j = 1,ip1jmp1
     109           pteta(j,k)= pteta(j,k)*(1+pdpfi(j)*pdt/pps(j))**kappa
     110        ENDDO
     111      ENDDO
     112!***********************
    106113
    107114      DO k = 1,llm

  • trunk/LMDZ.COMMON/libf/dyn3d/inidissip.F90

    r1021 r1238  
    268268    DO l=1,llm
    269269     zz      = 1. - preff/presnivs(l)
    270      zvert(l)= dissip_fac_mid -( dissip_fac_mid-1.)/( 1.+zz*zz )
     270!     zvert(l)= dissip_fac_mid -( dissip_fac_mid-1.)/( 1.+zz*zz )
     271     zvert(l)= fac_mid -( fac_mid-1.)/( 1.+zz*zz )
    271272     
    272273     zvert(l)= zvert(l)*(1.0+((fac_up/fac_mid-1)*    &
    … …  
    275276                ))
    276277    ENDDO
    277 
     278    write(lunout,*) "inidissip: vert_prof_disip=1, scaleheight=",scaleheight
     279    write(lunout,*) "           fac_mid=",fac_mid,", fac_up=",fac_up
     280   
    278281   else
    279282     write(lunout,*) 'wrong value for vert_prof_dissip:',vert_prof_dissip

  • trunk/LMDZ.COMMON/libf/dyn3dpar/addfi_p.F

    r1189 r1238  
    128128      endif
    129129c
     130!***********************
     131! Correction on teta due to surface pressure changes
     132c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     133      DO k = 1,llm
     134        DO j = ijb,ije
     135           pteta(j,k)= pteta(j,k)*(1+pdpfi(j)*pdt/pps(j))**kappa
     136        ENDDO
     137      ENDDO
     138c$OMP END DO NOWAIT
     139!***********************
    130140
    131141      ijb=ij_begin

  • trunk/LMDZ.COMMON/libf/dyn3dpar/inidissip.F90

    r1021 r1238  
    268268    DO l=1,llm
    269269     zz      = 1. - preff/presnivs(l)
    270      zvert(l)= dissip_fac_mid -( dissip_fac_mid-1.)/( 1.+zz*zz )
     270!     zvert(l)= dissip_fac_mid -( dissip_fac_mid-1.)/( 1.+zz*zz )
     271     zvert(l)= fac_mid -( fac_mid-1.)/( 1.+zz*zz )
    271272     
    272273     zvert(l)= zvert(l)*(1.0+((fac_up/fac_mid-1)*    &
    … …  
    275276                ))
    276277    ENDDO
    277 
     278    write(lunout,*) "inidissip: vert_prof_disip=1, scaleheight=",scaleheight
     279    write(lunout,*) "           fac_mid=",fac_mid,", fac_up=",fac_up
     280   
    278281   else
    279282     write(lunout,*) 'wrong value for vert_prof_dissip:',vert_prof_dissip
Note: See TracChangeset for help on using the changeset viewer.