<p><b>laura@ucar.edu</b> 2013-03-07 16:27:25 -0700 (Thu, 07 Mar 2013)</p><p>Added the local variables ctopo_p,ctopo2_p, and regime_p. These variables are declared as optional in module_bl_ysu.F, but need to be passed in the argument list. Results are unchanged relative to previous revision (tested using pgi DEBUG=true on yellowstone. If not in the argument lsit, MPAS crashes when compiled with gfortran (but not with pgi).<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_pbl.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_pbl.F        2013-03-07 23:24:55 UTC (rev 2568)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_pbl.F        2013-03-07 23:27:25 UTC (rev 2569)
@@ -24,11 +24,14 @@
 
 !from surface-layer model:
  if(.not.allocated(br_p)      ) allocate(br_p(ims:ime,jms:jme)               )
+ if(.not.allocated(ctopo_p)   ) allocate(ctopo_p(ims:ime,jms:jme)            )
+ if(.not.allocated(ctopo2_p)  ) allocate(ctopo2_p(ims:ime,jms:jme)           )
  if(.not.allocated(gz1oz0_p)  ) allocate(gz1oz0_p(ims:ime,jms:jme)           )
  if(.not.allocated(hfx_p)     ) allocate(hfx_p(ims:ime,jms:jme)              )
  if(.not.allocated(psih_p)    ) allocate(psih_p(ims:ime,jms:jme)             )
  if(.not.allocated(psim_p)    ) allocate(psim_p(ims:ime,jms:jme)             )
  if(.not.allocated(qfx_p)     ) allocate(qfx_p(ims:ime,jms:jme)              )
+ if(.not.allocated(regime_p)  ) allocate(regime_p(ims:ime,jms:jme)           )
  if(.not.allocated(u10_p)     ) allocate(u10_p(ims:ime,jms:jme)              )
  if(.not.allocated(ust_p)     ) allocate(ust_p(ims:ime,jms:jme)              )
  if(.not.allocated(v10_p)     ) allocate(v10_p(ims:ime,jms:jme)              )
@@ -59,11 +62,14 @@
 
 !from surface-layer model:
  if(allocated(br_p)      ) deallocate(br_p       )
+ if(allocated(ctopo_p)   ) deallocate(ctopo_p    )
+ if(allocated(ctopo2_p)  ) deallocate(ctopo2_p   )
  if(allocated(gz1oz0_p)  ) deallocate(gz1oz0_p   )
  if(allocated(hfx_p)     ) deallocate(hfx_p      )
  if(allocated(psih_p)    ) deallocate(psih_p     )
  if(allocated(psim_p)    ) deallocate(psim_p     )
  if(allocated(qfx_p)     ) deallocate(qfx_p      )
+ if(allocated(regime_p)  ) deallocate(regime_p   )
  if(allocated(u10_p)     ) deallocate(u10_p      )
  if(allocated(ust_p)     ) deallocate(ust_p      )
  if(allocated(v10_p)     ) deallocate(v10_p      )
@@ -108,6 +114,7 @@
     psim_p(i,j)   = diag_physics % psim   % array(i)
     psih_p(i,j)   = diag_physics % psih   % array(i)
     qfx_p(i,j)    = diag_physics % qfx    % array(i)
+    regime_p(i,j) = diag_physics % regime % array(i)
     u10_p(i,j)    = diag_physics % u10    % array(i)
     ust_p(i,j)    = diag_physics % ust    % array(i)
     v10_p(i,j)    = diag_physics % v10    % array(i)
@@ -116,6 +123,8 @@
 
     xland_p(i,j)  = sfc_input % xland  % array(i)       
 !initialization for YSU PBL scheme:
+    ctopo_p(i,j)  = 1._RKIND
+    ctopo2_p(i,j) = 1._RKIND
     kpbl_p(i,j)   = 1
  enddo
  enddo
@@ -214,8 +223,9 @@
                  psih     = psih_p     , xland    = xland_p     , hfx      = hfx_p      , &amp;
                  qfx      = qfx_p      , wspd     = wspd_p      , br       = br_p       , &amp;
                  dt       = dt_pbl     , kpbl2d   = kpbl_p      , exch_h   = exch_p     , &amp;
-                 u10      = u10_p       , v10      = v10_p      , kzhout   = kzh_p      , &amp;
-                 kzmout   = kzm_p       , kzqout   = kzq_p      , rho      = rho_p      , &amp;
+                 u10      = u10_p      , v10      = v10_p       , ctopo    = ctopo_p    , &amp; 
+                 ctopo2   = ctopo2_p   , regime   = regime_p    , rho      = rho_p      , &amp;
+                 kzhout   = kzh_p      , kzmout   = kzm_p       , kzqout   = kzq_p      , &amp;
                  ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde  , &amp;
                  ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme  , &amp;
                  its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte    &amp;
@@ -237,8 +247,9 @@
                  psih     = psih_p     , xland    = xland_p     , hfx      = hfx_p      , &amp;
                  qfx      = qfx_p      , wspd     = wspd_p      , br       = br_p       , &amp;
                  dt       = dt_pbl     , kpbl2d   = kpbl_p      , exch_h   = exch_p     , &amp;
-                 u10      = u10_p       , v10      = v10_p      , kzhout   = kzh_p      , &amp;
-                 kzmout   = kzm_p       , kzqout   = kzq_p      , rho      = rho_p      , &amp;
+                 u10      = u10_p      , v10      = v10_p       , ctopo    = ctopo_p    , &amp; 
+                 ctopo2   = ctopo2_p   , regime   = regime_p    , rho      = rho_p      , &amp;
+                 kzhout   = kzh_p      , kzmout   = kzm_p       , kzqout   = kzq_p      , &amp;
                  ids = ids , ide = ide , jds = jds , jde = jde , kds = kds , kde = kde  , &amp;
                  ims = ims , ime = ime , jms = jms , jme = jme , kms = kms , kme = kme  , &amp;
                  its = its , ite = ite , jts = jts , jte = jte , kts = kts , kte = kte    &amp;

</font>
</pre>