From 42d3aa5b76375bb60d304a029cbb69b756f89c12 Mon Sep 17 00:00:00 2001 From: Scott Rabenhorst Date: Sat, 4 May 2024 11:32:46 -0400 Subject: [PATCH 1/2] Add moist exports for mass fraction of suspended rain, snow, and graupel --- .../GEOS_BACM_1M_InterfaceMod.F90 | 9 ++ .../GEOS_GFDL_1M_InterfaceMod.F90 | 9 ++ .../GEOS_MGB2_2M_InterfaceMod.F90 | 9 ++ .../GEOSmoist_GridComp/GEOS_MoistGridComp.F90 | 82 ++++++++++++------- .../GEOS_THOM_1M_InterfaceMod.F90 | 9 ++ 5 files changed, 89 insertions(+), 29 deletions(-) diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_BACM_1M_InterfaceMod.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_BACM_1M_InterfaceMod.F90 index cb842e5f6..e0e1e43ff 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_BACM_1M_InterfaceMod.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_BACM_1M_InterfaceMod.F90 @@ -821,6 +821,15 @@ subroutine BACM_1M_Run (GC, IMPORT, EXPORT, CLOCK, RC) endif endif + call MAPL_GetPointer(EXPORT, PTR3D, 'QRTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = RAD_QR*RAD_CF + + call MAPL_GetPointer(EXPORT, PTR3D, 'QSTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = RAD_QS*RAD_CF + + call MAPL_GetPointer(EXPORT, PTR3D, 'QGTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = RAD_QG*RAD_CF + call MAPL_TimerOff (MAPL,"--BACM_1M") end subroutine BACM_1M_Run diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_GFDL_1M_InterfaceMod.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_GFDL_1M_InterfaceMod.F90 index e23df955c..4e90dbdd1 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_GFDL_1M_InterfaceMod.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_GFDL_1M_InterfaceMod.F90 @@ -962,6 +962,15 @@ subroutine GFDL_1M_Run (GC, IMPORT, EXPORT, CLOCK, RC) endif + call MAPL_GetPointer(EXPORT, PTR3D, 'QRTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = QRAIN + + call MAPL_GetPointer(EXPORT, PTR3D, 'QSTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = QSNOW + + call MAPL_GetPointer(EXPORT, PTR3D, 'QGTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = QGRAUPEL + call MAPL_TimerOff(MAPL,"--GFDL_1M",RC=STATUS) end subroutine GFDL_1M_Run diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MGB2_2M_InterfaceMod.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MGB2_2M_InterfaceMod.F90 index 7d2070ebf..abc6ece7e 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MGB2_2M_InterfaceMod.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MGB2_2M_InterfaceMod.F90 @@ -2494,6 +2494,15 @@ subroutine MGB2_2M_Run (GC, IMPORT, EXPORT, CLOCK, RC) endif endif + call MAPL_GetPointer(EXPORT, PTR3D, 'QRTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = QRAIN + + call MAPL_GetPointer(EXPORT, PTR3D, 'QSTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = QSNOW + + call MAPL_GetPointer(EXPORT, PTR3D, 'QGTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = QGRAUPEL + call MAPL_TimerOff(MAPL,"--MGB2_2M",__RC__) end subroutine MGB2_2M_Run diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 index 46c83c19d..48deccc3d 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 @@ -913,29 +913,53 @@ subroutine SetServices ( GC, RC ) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'QCTOT', & - LONG_NAME = 'mass_fraction_of_total_cloud_water', & - UNITS = 'kg kg-1', & + SHORT_NAME = 'QCTOT', & + LONG_NAME = 'mass_fraction_of_total_cloud_water', & + UNITS = 'kg kg-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'QLTOT', & - LONG_NAME = 'grid_box_mass_fraction_of_cloud_liquid_water', & - UNITS = 'kg kg-1', & - DIMS = MAPL_DimsHorzVert, & + SHORT_NAME = 'QLTOT', & + LONG_NAME = 'grid_box_mass_fraction_of_cloud_liquid_water', & + UNITS = 'kg kg-1', & + DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'QITOT', & - LONG_NAME = 'grid_box_mass_fraction_of_cloud_ice_water', & - UNITS = 'kg kg-1', & + SHORT_NAME = 'QITOT', & + LONG_NAME = 'grid_box_mass_fraction_of_cloud_ice_water', & + UNITS = 'kg kg-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) + call MAPL_AddExportSpec(GC, & + SHORT_NAME = 'QRTOT', & + LONG_NAME = 'grid_box_mass_fraction_of_suspended_rain', & + UNITS = 'kg kg-1', & + DIMS = MAPL_DimsHorzVert, & + VLOCATION = MAPL_VLocationCenter, RC=STATUS ) + VERIFY_(STATUS) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME = 'QSTOT', & + LONG_NAME = 'grid_box_mass_fraction_of_suspended_snow', & + UNITS = 'kg kg-1', & + DIMS = MAPL_DimsHorzVert, & + VLOCATION = MAPL_VLocationCenter, RC=STATUS ) + VERIFY_(STATUS) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME = 'QGTOT', & + LONG_NAME = 'grid_box_mass_fraction_of_suspended_graupel', & + UNITS = 'kg kg-1', & + DIMS = MAPL_DimsHorzVert, & + VLOCATION = MAPL_VLocationCenter, RC=STATUS ) + VERIFY_(STATUS) + call MAPL_AddExportSpec(GC, & SHORT_NAME = 'QPTOTLS', & LONG_NAME = 'mass_fraction_of_large_scale_falling_precip', & @@ -1880,15 +1904,15 @@ subroutine SetServices ( GC, RC ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'FCLD' , & - LONG_NAME = 'cloud_fraction_for_radiation', & - UNITS = '1', & + SHORT_NAME = 'FCLD' , & + LONG_NAME = 'cloud_fraction_for_radiation', & + UNITS = '1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME ='QV', & + SHORT_NAME ='QV', & LONG_NAME ='water_vapor_for_radiation', & UNITS ='kg kg-1', & DIMS = MAPL_DimsHorzVert, & @@ -1896,41 +1920,41 @@ subroutine SetServices ( GC, RC ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'QL', & - LONG_NAME = 'in_cloud_cloud_liquid_for_radiation', & - UNITS = 'kg kg-1', & + SHORT_NAME = 'QL', & + LONG_NAME = 'in_cloud_cloud_liquid_for_radiation', & + UNITS = 'kg kg-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'QI', & - LONG_NAME = 'in_cloud_cloud_ice_for_radiation', & - UNITS = 'kg kg-1', & + SHORT_NAME = 'QI', & + LONG_NAME = 'in_cloud_cloud_ice_for_radiation', & + UNITS = 'kg kg-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'QR', & - LONG_NAME = 'Falling_rain_for_radiation', & - UNITS = 'kg kg-1', & + SHORT_NAME = 'QR', & + LONG_NAME = 'in_cloud_ralling_rain_for_radiation', & + UNITS = 'kg kg-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'QS', & - LONG_NAME = 'Falling_snow_for_radiation', & - UNITS = 'kg kg-1', & + SHORT_NAME = 'QS', & + LONG_NAME = 'in_cloud_falling_snow_for_radiation', & + UNITS = 'kg kg-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'QG', & - LONG_NAME = 'Falling_graupel_for_radiation', & - UNITS = 'kg kg-1', & + SHORT_NAME = 'QG', & + LONG_NAME = 'in_cloud_falling_graupel_for_radiation', & + UNITS = 'kg kg-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_THOM_1M_InterfaceMod.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_THOM_1M_InterfaceMod.F90 index 9ff2269a0..4a84fcbe1 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_THOM_1M_InterfaceMod.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_THOM_1M_InterfaceMod.F90 @@ -1132,6 +1132,15 @@ subroutine THOM_1M_Run (GC, IMPORT, EXPORT, CLOCK, RC) endif + call MAPL_GetPointer(EXPORT, PTR3D, 'QRTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = QRAIN + + call MAPL_GetPointer(EXPORT, PTR3D, 'QSTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = QSNOW + + call MAPL_GetPointer(EXPORT, PTR3D, 'QGTOT', RC=STATUS); VERIFY_(STATUS) + if (associated(PTR3D)) PTR3D = QGRAUPEL + call MAPL_TimerOff(MAPL,"--THOM_1M",RC=STATUS) end subroutine THOM_1M_Run From db55c301840d98b788b0e17045510af726c0f555 Mon Sep 17 00:00:00 2001 From: Matthew Thompson Date: Mon, 6 May 2024 10:06:55 -0400 Subject: [PATCH 2/2] Fix misspelling --- .../GEOSmoist_GridComp/GEOS_MoistGridComp.F90 | 80 +++++++++---------- 1 file changed, 40 insertions(+), 40 deletions(-) diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 index 48deccc3d..c43eb1e79 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 @@ -524,7 +524,7 @@ subroutine SetServices ( GC, RC ) SHORT_NAME = 'TAUOROX', & LONG_NAME = 'surface_eastward_orographic_gravity_wave_stress', & UNITS = 'N m-2', & - RESTART = MAPL_RestartSkip, & + RESTART = MAPL_RestartSkip, & DIMS = MAPL_DimsHorzOnly, & VLOCATION = MAPL_VLocationNone, RC=STATUS ) VERIFY_(STATUS) @@ -533,7 +533,7 @@ subroutine SetServices ( GC, RC ) SHORT_NAME = 'TAUOROY', & LONG_NAME = 'surface_northward_orographic_gravity_wave_stress', & UNITS = 'N m-2', & - RESTART = MAPL_RestartSkip, & + RESTART = MAPL_RestartSkip, & DIMS = MAPL_DimsHorzOnly, & VLOCATION = MAPL_VLocationNone, RC=STATUS ) VERIFY_(STATUS) @@ -543,7 +543,7 @@ subroutine SetServices ( GC, RC ) SHORT_NAME = 'OMEGA', & LONG_NAME = 'vertical_pressure_velocity', & UNITS = 'Pa s-1', & - RESTART = MAPL_RestartSkip, & + RESTART = MAPL_RestartSkip, & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) VERIFY_(STATUS) @@ -940,25 +940,25 @@ subroutine SetServices ( GC, RC ) SHORT_NAME = 'QRTOT', & LONG_NAME = 'grid_box_mass_fraction_of_suspended_rain', & UNITS = 'kg kg-1', & - DIMS = MAPL_DimsHorzVert, & - VLOCATION = MAPL_VLocationCenter, RC=STATUS ) - VERIFY_(STATUS) + DIMS = MAPL_DimsHorzVert, & + VLOCATION = MAPL_VLocationCenter, RC=STATUS ) + VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & SHORT_NAME = 'QSTOT', & LONG_NAME = 'grid_box_mass_fraction_of_suspended_snow', & UNITS = 'kg kg-1', & - DIMS = MAPL_DimsHorzVert, & - VLOCATION = MAPL_VLocationCenter, RC=STATUS ) - VERIFY_(STATUS) + DIMS = MAPL_DimsHorzVert, & + VLOCATION = MAPL_VLocationCenter, RC=STATUS ) + VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & SHORT_NAME = 'QGTOT', & LONG_NAME = 'grid_box_mass_fraction_of_suspended_graupel', & UNITS = 'kg kg-1', & - DIMS = MAPL_DimsHorzVert, & - VLOCATION = MAPL_VLocationCenter, RC=STATUS ) - VERIFY_(STATUS) + DIMS = MAPL_DimsHorzVert, & + VLOCATION = MAPL_VLocationCenter, RC=STATUS ) + VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & SHORT_NAME = 'QPTOTLS', & @@ -1146,18 +1146,18 @@ subroutine SetServices ( GC, RC ) RC=STATUS ) VERIFY_(STATUS) - call MAPL_AddExportSpec(GC, & + call MAPL_AddExportSpec(GC, & SHORT_NAME = 'UMF_DC', & - LONG_NAME = 'Deep_updraft_mass_flux_at_interfaces', & + LONG_NAME = 'Deep_updraft_mass_flux_at_interfaces', & UNITS = 'kg m-2 s-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationEdge, & RC=STATUS ) VERIFY_(STATUS) - - call MAPL_AddExportSpec(GC, & + + call MAPL_AddExportSpec(GC, & SHORT_NAME = 'MFD_DC', & - LONG_NAME = 'Deep_updraft_detrained_mass_flux', & + LONG_NAME = 'Deep_updraft_detrained_mass_flux', & UNITS = 'kg m-2 s-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, & @@ -1192,7 +1192,7 @@ subroutine SetServices ( GC, RC ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'CNV_TOPP_DP', & + SHORT_NAME = 'CNV_TOPP_DP', & LONG_NAME = 'pressure_at_deep_convective_cloud_top', & UNITS = 'Pa', & DIMS = MAPL_DimsHorzOnly, & @@ -1201,7 +1201,7 @@ subroutine SetServices ( GC, RC ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'CNV_TOPP_MD', & + SHORT_NAME = 'CNV_TOPP_MD', & LONG_NAME = 'pressure_at_congestus_convective_cloud_top', & UNITS = 'Pa', & DIMS = MAPL_DimsHorzOnly, & @@ -1210,7 +1210,7 @@ subroutine SetServices ( GC, RC ) VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & - SHORT_NAME = 'CNV_TOPP_SH', & + SHORT_NAME = 'CNV_TOPP_SH', & LONG_NAME = 'pressure_at_shallow_convective_cloud_top', & UNITS = 'Pa', & DIMS = MAPL_DimsHorzOnly, & @@ -1240,23 +1240,23 @@ subroutine SetServices ( GC, RC ) call MAPL_AddExportSpec(GC, & SHORT_NAME = 'ENTR_MD', & - LONG_NAME = 'Lateral_entrainment_rate_in_gf_mid_convection_updraft', & + LONG_NAME = 'Lateral_entrainment_rate_in_gf_mid_convection_updraft', & UNITS = 'm-1', & DIMS = MAPL_DimsHorzVert, & - VLOCATION = MAPL_VLocationCenter, & + VLOCATION = MAPL_VLocationCenter, & RC=STATUS ) - - VERIFY_(STATUS) + + VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & SHORT_NAME = 'ENTR_SH', & - LONG_NAME = 'Lateral_entrainment_rate_in_gf_shallow_convection_updraft', & + LONG_NAME = 'Lateral_entrainment_rate_in_gf_shallow_convection_updraft', & UNITS = 'm-1', & DIMS = MAPL_DimsHorzVert, & - VLOCATION = MAPL_VLocationCenter, & + VLOCATION = MAPL_VLocationCenter, & RC=STATUS ) - - VERIFY_(STATUS) + + VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & SHORT_NAME = 'CNV_UPDF', & @@ -1937,7 +1937,7 @@ subroutine SetServices ( GC, RC ) call MAPL_AddExportSpec(GC, & SHORT_NAME = 'QR', & - LONG_NAME = 'in_cloud_ralling_rain_for_radiation', & + LONG_NAME = 'in_cloud_falling_rain_for_radiation', & UNITS = 'kg kg-1', & DIMS = MAPL_DimsHorzVert, & VLOCATION = MAPL_VLocationCenter, RC=STATUS ) @@ -2031,13 +2031,13 @@ subroutine SetServices ( GC, RC ) VLOCATION = MAPL_VLocationNone, RC=STATUS ) VERIFY_(STATUS) - call MAPL_AddExportSpec(GC, & + call MAPL_AddExportSpec(GC, & SHORT_NAME = 'SNOW_RATIO', & LONG_NAME = 'ratio_of_snow_to_total_precip',& UNITS = '1', & DIMS = MAPL_DimsHorzOnly, & VLOCATION = MAPL_VLocationNone, RC=STATUS ) - VERIFY_(STATUS) + VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & SHORT_NAME ='LS_PRCP', & @@ -2571,13 +2571,13 @@ subroutine SetServices ( GC, RC ) VLOCATION = MAPL_VLocationNone, RC=STATUS ) VERIFY_(STATUS) - call MAPL_AddExportSpec(GC, & + call MAPL_AddExportSpec(GC, & SHORT_NAME='RKFRE', & - LONG_NAME ='fraction_of_tke_associated_with_vertical_velocity', & - UNITS ='' , & + LONG_NAME ='fraction_of_tke_associated_with_vertical_velocity', & + UNITS ='' , & DIMS = MAPL_DimsHorzOnly, & VLOCATION = MAPL_VLocationNone, RC=STATUS ) - VERIFY_(STATUS) + VERIFY_(STATUS) call MAPL_AddExportSpec(GC, & SHORT_NAME='STOCH_CNV', & @@ -5389,19 +5389,19 @@ subroutine RUN ( GC, IMPORT, EXPORT, CLOCK, RC ) ! Trajectory for the moist TLM/ADJ !--------------------------------------------------------- call MAPL_GetPointer(EXPORT, PTR3D, 'TH_moist' , RC=STATUS); VERIFY_(STATUS) - if (associated(PTR3D)) PTR3D = T/PK + if (associated(PTR3D)) PTR3D = T/PK call MAPL_GetPointer(EXPORT, PTR3D, 'Q_moist' , RC=STATUS); VERIFY_(STATUS) if (associated(PTR3D)) PTR3D = Q call MAPL_GetPointer(EXPORT, PTR2D, 'TS_moist' , RC=STATUS); VERIFY_(STATUS) - if (associated(PTR2D)) PTR2D = TS + if (associated(PTR2D)) PTR2D = TS call MAPL_GetPointer(EXPORT, PTR2D, 'KHu_moist', RC=STATUS); VERIFY_(STATUS) if (associated(PTR2D)) then - PTR2D = -1 + PTR2D = -1 do l = 0,LM where( (PTR2D == -1) .AND. (KH(:,:,l) > 2.0) ) PTR2D = l * 1.0 end where - enddo + enddo endif call MAPL_GetPointer(EXPORT, PTR2D, 'KHl_moist', RC=STATUS); VERIFY_(STATUS) if (associated(PTR2D)) then @@ -5427,7 +5427,7 @@ subroutine RUN ( GC, IMPORT, EXPORT, CLOCK, RC ) endif PTR2D(I,J) = MAX( PTR2D(I,J), KCBLMIN ) enddo - enddo + enddo endif ! These may be used by children