[Dart-dev] DART/branches Revision: 12386

dart at ucar.edu dart at ucar.edu
Thu Feb 1 15:59:15 MST 2018


hendric at ucar.edu
2018-02-01 15:59:14 -0700 (Thu, 01 Feb 2018)
44

had the inverse scale height computation.




Modified: DART/branches/recam/models/cam-fv/model_mod.f90
===================================================================
--- DART/branches/recam/models/cam-fv/model_mod.f90	2018-02-01 21:45:43 UTC (rev 12385)
+++ DART/branches/recam/models/cam-fv/model_mod.f90	2018-02-01 22:59:14 UTC (rev 12386)
@@ -804,15 +804,15 @@
 
 !print*, 'STATUS MI EVAL', status_array(1)
 
-! print*, 'lon_ind_below ', four_lons(1)
-! print*, 'lon_ind_above ', four_lons(2)
-! print*, 'lat_ind_below ', four_lats(2)
-! print*, 'lat_ind_above ', four_lats(3)
-! print*, 'lon_fract     ', lon_fract
-! print*, 'lat_fract     ', lat_fract
-! print*, 'quad_vals(:,1) ', quad_vals(:,1)
-! print*, 'quad_vals(:,2) ', quad_vals(:,2)
-! print*, 'quad_vals(:,3) ', quad_vals(:,3)
+! print*, 'lon_ind_below    ', four_lons(1)
+! print*, 'lon_ind_above    ', four_lons(2)
+! print*, 'lat_ind_below    ', four_lats(2)
+! print*, 'lat_ind_above    ', four_lats(3)
+! print*, 'lon_fract        ', lon_fract
+! print*, 'lat_fract        ', lat_fract
+! print*, 'quad_vals(:,1)   ', quad_vals(:,1)
+! print*, 'quad_vals(:,2)   ', quad_vals(:,2)
+! print*, 'quad_vals(:,3)   ', quad_vals(:,3)
 ! print*, 'interp_vals(:,1) ', interp_vals(1)
 ! print*, 'interp_vals(:,2) ', interp_vals(2)
 ! print*, 'interp_vals(:,3) ', interp_vals(3)
@@ -3958,10 +3958,12 @@
 
 else
    ! normal computation - should be safe now
-   scale_height = -log(p_surface / p_above)
+   scale_height = -log(p_above / p_surface )
 
 endif
 
+print*, 'p_surface, p_above', p_surface, p_above, diff
+
 end function scale_height
 
 !--------------------------------------------------------------------


More information about the Dart-dev mailing list