Changeset 159
- Timestamp:
- Jun 16, 2011, 11:43:20 AM (14 years ago)
- Location:
- trunk/MESOSCALE/LMD_MM_MARS/SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MESOSCALE/LMD_MM_MARS/SRC/LESnophys_/modif/module_initialize_les.F
r156 r159 164 164 ! stretch_grid = .true. 165 165 ! FOR LES, set stretch to false 166 #ifdef NOPHYS 167 stretch_grid = .true. 168 #else 166 169 stretch_grid = .false. 170 #endif 167 171 delt = 3. 168 172 ! z_scale = .50 … … 233 237 234 238 IF (stretch_grid) THEN ! exponential stretch for eta (nearly constant dz) 239 print *, '*** STRETCH GRID ***' 235 240 DO k=1, kde 236 241 grid%znw(k) = (exp(-(k-1)/float(kde-1)/z_scale) - exp(-1./z_scale))/ & -
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd_new_storm/libf/phymars/meso_physiq.F
r118 r159 211 211 logical tracerdyn 212 212 CHARACTER (len=20) :: wtnom(nqmx) ! tracer name 213 REAL toto 213 214 214 215 c outputs: … … 725 726 $ zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,fluxtop_sw, 726 727 & tauref,tau,aerosol,ccn,rdust,rice,nuice,zdqnorm,dsodust) 728 729 !!!!!!!!!!!!!!!!!!!!!! TEST TEST 730 ! DO l=2,nlayer 731 ! DO ig=1,ngrid 732 ! toto = abs(rdust(ig,l) - rdust(ig,l-1)) 733 ! if ( toto .gt. 0.5*rdust(ig,l) ) then 734 ! write (*,*) "LIMITING RADIUS", ig, l 735 ! write (*,*) rdust(ig,l-1), rdust(ig,l) 736 ! rdust(ig,l) = rdust(ig,l-1) 737 ! endif 738 ! ENDDO 739 ! ENDDO 740 ! 741 ! DO l=2,nlayer 742 ! DO ig=2,ngrid 743 ! toto = abs(rdust(ig-1,l) - rdust(ig,l))/rdust(ig,l) 744 ! if ( toto .gt. 0.45 ) then 745 ! write (*,*) "RDUST ", ig, l, rdust(ig-1,l), rdust(ig,l) 746 ! rdust(ig,l) = rdust(ig,l-1) 747 ! endif 748 ! ENDDO 749 ! ENDDO 727 750 728 751
Note: See TracChangeset
for help on using the changeset viewer.