From 17d1507ac61d377dd4251f0b240a72db5f11d694 Mon Sep 17 00:00:00 2001 From: Arlindo da Silva Date: Tue, 2 Nov 2021 18:31:30 -0400 Subject: [PATCH] Merged in GOCART-2G updates --- .../GEOS_PhysicsGridComp.F90 | 4 +- .../GEOSmoist_GridComp/ConvPar_GF_GEOS5.F90 | 1 + .../GEOSmoist_GridComp/GEOS_MoistGridComp.F90 | 238 ++++++++++++++++-- .../GEOSirrad_GridComp/GEOS_IrradGridComp.F90 | 2 +- .../GEOSsolar_GridComp/GEOS_SolarGridComp.F90 | 3 +- 5 files changed, 228 insertions(+), 20 deletions(-) diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOS_PhysicsGridComp.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOS_PhysicsGridComp.F90 index 336ff0896..bc87a02b4 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOS_PhysicsGridComp.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOS_PhysicsGridComp.F90 @@ -1315,7 +1315,9 @@ subroutine SetServices ( GC, RC ) !Aerosol call MAPL_AddConnectivity ( GC, & - SHORT_NAME = (/'AERO_ACI'/), & + SHORT_NAME = (/'AERO'/), & +! SHORT_NAME = (/'AERO_ACI '/), & +! SHORT_NAME = (/'AERO_ACI ', 'AERO2G_ACI'/), & DST_ID = MOIST, & SRC_ID = CHEM, & RC=STATUS ) diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/ConvPar_GF_GEOS5.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/ConvPar_GF_GEOS5.F90 index 0e4ea9207..a7cb89de7 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/ConvPar_GF_GEOS5.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/ConvPar_GF_GEOS5.F90 @@ -9939,6 +9939,7 @@ SUBROUTINE get_incloud_sc_chem_up(cumulus,fscav,mtp,se,se_cup,sc_up,pw_up,tot_pw if( TRIM(CHEM_name (ispc)) == 'sulfur' .or. & TRIM(CHEM_name (ispc)(1:len_trim('ss') )) == 'ss' .or. & ! 'seasalt' + TRIM(CHEM_name (ispc)(1:len_trim('SS') )) == 'SS' .or. & ! 'seasalt 2G' TRIM(CHEM_name (ispc)) == 'SO2' .or. & TRIM(CHEM_name (ispc)) == 'SO4' .or. & diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 index e31dcf4fd..301a7ef0c 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSmoist_GridComp/GEOS_MoistGridComp.F90 @@ -730,13 +730,13 @@ subroutine SetServices ( GC, RC ) RC=STATUS ) VERIFY_(STATUS) - call MAPL_AddImportSpec(GC, & - LONG_NAME = 'aerosol_cloud_interaction', & - UNITS = '1', & - SHORT_NAME = 'AERO_ACI', & - DIMS = MAPL_DimsHorzVert, & - VLOCATION = MAPL_VLocationCenter, & - DATATYPE = MAPL_StateItem, & + call MAPL_AddImportSpec(GC, & + LONG_NAME = 'aerosols', & + UNITS = '1', & + SHORT_NAME = 'AERO', & + DIMS = MAPL_DimsHorzVert, & + VLOCATION = MAPL_VLocationCenter, & + DATATYPE = MAPL_StateItem, & RC=STATUS ) VERIFY_(STATUS) @@ -3560,6 +3560,73 @@ subroutine SetServices ( GC, RC ) ! CAR 12/5/08 ! Aerosol Scavenging diagnostics/export states ! ------------------------------ + call MAPL_AddExportSpec(GC, & + SHORT_NAME='DDU2gDT ', & + LONG_NAME ='dust_tendency_due_to_conv_scav', & + UNITS ='kg m-2 s-1', & + DIMS = MAPL_DimsHorzOnly, & + __RC__ ) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME='DSS2gDT ', & + LONG_NAME ='sea_salt_tendency_due_to_conv_scav', & + UNITS ='kg m-2 s-1', & + DIMS = MAPL_DimsHorzOnly, & + __RC__ ) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME='DOC2gDT ', & + LONG_NAME ='organic_carbon_tendency_due_to_conv_scav', & + UNITS ='kg m-2 s-1', & + DIMS = MAPL_DimsHorzOnly, & + __RC__ ) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME='DBC2gDT ', & + LONG_NAME ='black_carbon_tendency_due_to_conv_scav', & + UNITS ='kg m-2 s-1', & + DIMS = MAPL_DimsHorzOnly, & + __RC__ ) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME='DSU2gDT ', & + LONG_NAME ='sulfate_tendency_due_to_conv_scav', & + UNITS ='kg m-2 s-1', & + DIMS = MAPL_DimsHorzOnly, & + __RC__ ) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME='DNI2gDT', & + LONG_NAME ='nitrate_tendency_due_to_conv_scav', & + UNITS ='kg m-2 s-1', & + DIMS = MAPL_DimsHorzOnly, & + __RC__ ) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME='DNH4A2gDT', & + LONG_NAME ='ammonium_aerosol_tendency_due_to_conv_scav', & + UNITS ='kg m-2 s-1', & + DIMS = MAPL_DimsHorzOnly, & + __RC__ ) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME='DNH32gDT', & + LONG_NAME ='ammonia_tendency_due_to_conv_scav', & + UNITS ='kg m-2 s-1', & + DIMS = MAPL_DimsHorzOnly, & + __RC__ ) + + call MAPL_AddExportSpec(GC, & + SHORT_NAME='DBRC2gDT', & + LONG_NAME ='brown_carbon_tendency_due_to_conv_scav', & + UNITS ='kg m-2 s-1', & + DIMS = MAPL_DimsHorzOnly, & + __RC__ ) + + + + + call MAPL_AddExportSpec(GC, & SHORT_NAME='DDUDT ', & LONG_NAME ='dust_tendency_due_to_conv_scav', & @@ -5301,6 +5368,8 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) ! CAR real, pointer, dimension(:,: ) :: DDUDT, & DSSDT, DOCDT, DBCDT, DSUDT, DNIDT, DNH4ADT, DNH3DT, DBRCDT, DDUDTcarma, DSSDTcarma + real, pointer, dimension(:,: ) :: DDU2gDT, & + DSS2gDT, DOC2gDT, DBC2gDT, DSU2gDT, DNI2gDT, DNH4A2gDT, DNH32gDT, DBRC2gDT character(len=ESMF_MAXSTR) :: QNAME, CNAME, ENAME character(len=ESMF_MAXSTR), pointer, dimension(:) :: QNAMES, CNAMES integer :: ind @@ -5562,6 +5631,7 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) character(len=ESMF_MAXSTR), allocatable, dimension(:) :: aero_aci_modes integer :: ACI_STATUS + type (AerProps), dimension (IM, JM, LM) :: AeroProps !Storages aerosol properties for activation type (AerProps) :: AeroAux, AeroAux_b @@ -5977,6 +6047,7 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) real :: cNN, cNN_OCEAN, cNN_LAND, CONVERT real , dimension(IM,JM) :: CMDU, CMSS, CMOC, CMBC, CMSU, CMNI, CMNH3, CMNH4A, CMBRC + real , dimension(IM,JM) :: CMDU2g, CMSS2g, CMOC2g, CMBC2g, CMSU2g, CMNI2g, CMNH32g, CMNH4A2g, CMBRC2g real , dimension(IM,JM) :: CMDUcarma, CMSScarma real :: sigmaqt, qcn, cfn, qsatn, dqlls, dqils, qt @@ -6627,6 +6698,16 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) call MAPL_GetPointer(EXPORT, ENTLAM, 'ENTLAM' , RC=STATUS); VERIFY_(STATUS) ! Aerosol Scavenging ! CAR 12/5/08 + call MAPL_GetPointer(EXPORT, DDU2gDT, 'DDU2gDT' , RC=STATUS); VERIFY_(STATUS) + call MAPL_GetPointer(EXPORT, DSS2gDT, 'DSS2gDT' , RC=STATUS); VERIFY_(STATUS) + call MAPL_GetPointer(EXPORT, DBC2gDT, 'DBC2gDT' , RC=STATUS); VERIFY_(STATUS) + call MAPL_GetPointer(EXPORT, DOC2gDT, 'DOC2gDT' , RC=STATUS); VERIFY_(STATUS) + call MAPL_GetPointer(EXPORT, DSU2gDT, 'DSU2gDT' , RC=STATUS); VERIFY_(STATUS) + call MAPL_GetPointer(EXPORT, DNI2gDT, 'DNI2gDT' , RC=STATUS); VERIFY_(STATUS) + call MAPL_GetPointer(EXPORT, DNH4A2gDT, 'DNH4A2gDT' , RC=STATUS); VERIFY_(STATUS) + call MAPL_GetPointer(EXPORT, DNH32gDT, 'DNH32gDT' , RC=STATUS); VERIFY_(STATUS) + call MAPL_GetPointer(EXPORT, DBRC2gDT, 'DBRC2gDT' , RC=STATUS); VERIFY_(STATUS) + call MAPL_GetPointer(EXPORT, DDUDT, 'DDUDT' , RC=STATUS); VERIFY_(STATUS) call MAPL_GetPointer(EXPORT, DSSDT, 'DSSDT' , RC=STATUS); VERIFY_(STATUS) call MAPL_GetPointer(EXPORT, DBCDT, 'DBCDT' , RC=STATUS); VERIFY_(STATUS) @@ -7605,12 +7686,12 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) end do - call ESMF_StateGet(IMPORT, 'AERO_ACI', aero_aci, __RC__) + call ESMF_StateGet(IMPORT, 'AERO', aero_aci, __RC__) - call ESMF_AttributeGet(aero_aci, name='implements_aerosol_activation_properties_method', & - value=implements_aerosol_activation_properties, __RC__) +! call ESMF_AttributeGet(aero_aci, name='implements_aerosol_activation_properties_method', & +! value=implements_aerosol_activation_properties, __RC__) - if (implements_aerosol_activation_properties) then +! if (implements_aerosol_activation_properties) then call ESMF_AttributeGet(aero_aci, name='number_of_aerosol_modes', value=n_modes, __RC__) @@ -7618,8 +7699,8 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) allocate(aero_aci_modes(n_modes), __STAT__) call ESMF_AttributeGet(aero_aci, name='aerosol_modes', itemcount=n_modes, valuelist=aero_aci_modes, __RC__) - - call ESMF_AttributeGet(aero_aci, name='air_pressure', value=aci_field_name, __RC__) + + call ESMF_AttributeGet(aero_aci, name='air_pressure_for_aerosol_optics', value=aci_field_name, __RC__) if (aci_field_name /= '') then call MAPL_GetPointer(aero_aci, aci_ptr_3d, trim(aci_field_name), __RC__) aci_ptr_3d = PLE @@ -7675,7 +7756,7 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) #if (0) if (MAPL_AM_I_ROOT()) then print * - print *, 'AERO_ACI::' // trim(aero_aci_modes(n)) + print *, 'AERO::' // trim(aero_aci_modes(n)) print *, 'num : ', aci_num(1,1,LM) print *, 'dgn : ', aci_dgn(1,1,LM) @@ -7725,11 +7806,11 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) call MAPL_TimerOff(STATE,"--USE_AEROSOL_NN1") - else +! else ! options: ! *) set aerosol concentrations to 0.0, i.e., no aerosol ! *) raise an exception if aerosol is required! - end if +! end if call init_Aer(AeroAux) @@ -8001,6 +8082,16 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) ! Compute initial mass loading for aerosols; CAR 12/19/08 ! ------------------------------------------------------- !! First initialize everything to zero + if(associated(DDU2gDT)) DDU2gDT = 0.0 + if(associated(DSS2gDT)) DSS2gDT = 0.0 + if(associated(DBC2gDT)) DBC2gDT = 0.0 + if(associated(DOC2gDT)) DOC2gDT = 0.0 + if(associated(DSU2gDT)) DSU2gDT = 0.0 + if(associated(DNI2gDT)) DNI2gDT = 0.0 + if(associated(DNH4A2gDT)) DNH4A2gDT = 0.0 + if(associated(DNH32gDT)) DNH32gDT = 0.0 + if(associated(DBRC2gDT)) DBRC2gDT= 0.0 + if(associated(DDUDT)) DDUDT = 0.0 if(associated(DSSDT)) DSSDT = 0.0 if(associated(DBCDT)) DBCDT = 0.0 @@ -8013,6 +8104,16 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) if(associated(DDUDTcarma)) DDUDTcarma = 0.0 if(associated(DSSDTcarma)) DSSDTcarma = 0.0 + CMDU2g = 0.0 + CMSS2g = 0.0 + CMOC2g = 0.0 + CMBC2g = 0.0 + CMSU2g = 0.0 + CMNI2g = 0.0 + CMNH4A2g = 0.0 + CMNH32g = 0.0 + CMBRC2g = 0.0 + CMDU = 0.0 CMSS = 0.0 CMOC = 0.0 @@ -8027,13 +8128,59 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) !! Now loop over tracers and accumulate initial column loading !! tendency kg/m2/s CAR - +!if(mapl_am_i_root()) print*,'MOIST CNAMES = ',CNAMES KK=0 do K=1,KM if(IS_FRIENDLY(K)) then KK = KK + 1 QNAME = trim(QNAMES(K)) CNAME = trim(CNAMES(K)) + if((CNAME == 'DU') .or. (CNAME == 'SS') .or. (CNAME == 'NI') .or. (CNAME == 'SU') .or. & + (CNAME == 'CA.oc') .or. (CNAME == 'CA.bc') .or. (CNAME == 'CA.br')) then ! Diagnostics for GOCART2G tracers + SELECT CASE (QNAME(1:3)) + CASE ('DU0') + if(associated(DDU2gDT)) then + CMDU2g = CMDU2g + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('SS0') + if(associated(DSS2gDT)) then + CMSS2g = CMSS2g + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('SO4') + if(associated(DSU2gDT)) then + CMSU2g = CMSU2g + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('NO3') + if(associated(DNI2gDT)) then + CMNI2g = CMNI2g + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('NH3') + if(associated(DNH32gDT)) then + CMNH32g = CMNH32g + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('NH4') + if(associated(DNH4A2gDT)) then + CMNH4A2g = CMNH4A2g + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + END SELECT + + SELECT CASE (QNAME(1:13)) + CASE ('CAphilicCA.bc') + if(associated(DBC2gDT)) then + CMBC2g = CMBC2g + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if +! CASE ('OCp') + CASE ('CAphilicCA.oc') + if(associated(DOC2gDT)) then + CMOC2g = CMOC2g + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('CAphilicCA.br') + if(associated(DBRC2gDT)) then + CMBRC2g = CMBRC2g + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + END SELECT + endif + if(CNAME == 'GOCART') then ! Diagnostics for GOCART tracers SELECT CASE (QNAME(1:3)) CASE ('du0') @@ -8074,6 +8221,7 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) end if END SELECT endif + if(CNAME == 'CARMA') then ! Diagnostics for CARMA tracers ! Check name to see if it is a "pc" element ENAME = '' @@ -8534,6 +8682,52 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) KK = KK + 1 QNAME = trim(QNAMES(K)) CNAME = trim(CNAMES(K)) + if((CNAME == 'DU') .or. (CNAME == 'SS') .or. (CNAME == 'NI') .or. (CNAME == 'SU') .or. & + (CNAME == 'CA.oc') .or. (CNAME == 'CA.bc') .or. (CNAME == 'CA.br')) then ! Diagnostics for GOCART2G tracers + SELECT CASE (QNAME(1:3)) + CASE ('DU0') + if(associated(DDU2gDT)) then + DDU2gDT = DDU2gDT + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('SS0') + if(associated(DSS2gDT)) then + DSS2gDT = DSS2gDT + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('SO4') + if(associated(DSU2gDT)) then + DSU2gDT = DSU2gDT + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('NO3') + if(associated(DNI2gDT)) then + DNI2gDT = DNI2gDT + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('NH3') + if(associated(DNH32gDT)) then + DNH32gDT = DNH32gDT + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('NH4') + if(associated(DNH4A2gDT)) then + DNH4A2gDT = DNH4A2gDT + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + END SELECT + + SELECT CASE (QNAME(1:13)) + CASE ('CAphilicCA.bc') + if(associated(DBC2gDT)) then + DBC2gDT = DBC2gDT + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if +! CASE ('OCp') + CASE ('CAphilicCA.oc') + if(associated(DOC2gDT)) then + DOC2gDT = DOC2gDT + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + CASE ('CAphilicCA.br') + if(associated(DBRC2gDT)) then + DBRC2gDT = DBRC2gDT + sum(XHO(:,:,:,KK)*DP(:,:,:),dim=3) + end if + END SELECT + endif + if(CNAME == 'GOCART') then ! Diagnostics for GOCART tracers SELECT CASE (QNAME(1:3)) CASE ('du0') @@ -8597,6 +8791,16 @@ subroutine MOIST_DRIVER(IM,JM,LM, RC) endif end do + if (associated(DDU2gDT)) DDU2gDT = (DDU2gDT - CMDU2g) / (MAPL_GRAV*DT_MOIST) + if (associated(DSS2gDT)) DSS2gDT = (DSS2gDT - CMSS2g) / (MAPL_GRAV*DT_MOIST) + if (associated(DBC2gDT)) DBC2gDT = (DBC2gDT - CMBC2g) / (MAPL_GRAV*DT_MOIST) + if (associated(DOC2gDT)) DOC2gDT = (DOC2gDT - CMOC2g) / (MAPL_GRAV*DT_MOIST) + if (associated(DSU2gDT)) DSU2gDT = (DSU2gDT - CMSU2g) / (MAPL_GRAV*DT_MOIST) + if (associated(DNI2gDT)) DNI2gDT = (DNI2gDT - CMNI2g) / (MAPL_GRAV*DT_MOIST) + if (associated(DNH32gDT)) DNH32gDT = (DNH32gDT - CMNH32g) / (MAPL_GRAV*DT_MOIST) + if (associated(DNH4A2gDT)) DNH4A2gDT = (DNH4A2gDT - CMNH4A2g) / (MAPL_GRAV*DT_MOIST) + if (associated(DBRC2gDT)) DBRC2gDT= (DBRC2gDT- CMBRC2g)/ (MAPL_GRAV*DT_MOIST) + if (associated(DDUDT)) DDUDT = (DDUDT - CMDU) / (MAPL_GRAV*DT_MOIST) if (associated(DSSDT)) DSSDT = (DSSDT - CMSS) / (MAPL_GRAV*DT_MOIST) if (associated(DBCDT)) DBCDT = (DBCDT - CMBC) / (MAPL_GRAV*DT_MOIST) diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSradiation_GridComp/GEOSirrad_GridComp/GEOS_IrradGridComp.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSradiation_GridComp/GEOSirrad_GridComp/GEOS_IrradGridComp.F90 index e2fe430cf..c51a5d9c3 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSradiation_GridComp/GEOSirrad_GridComp/GEOS_IrradGridComp.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSradiation_GridComp/GEOSirrad_GridComp/GEOS_IrradGridComp.F90 @@ -1939,7 +1939,7 @@ subroutine Lw_Driver(IM,JM,LM,LATS,LONS,CoresPerNode,RC) VERIFY_(STATUS) ! execute the aero provider's optics method - call ESMF_MethodExecute(AERO, label="aerosol_optics", RC=STATUS) + call ESMF_MethodExecute(AERO, label="run_aerosol_optics", RC=STATUS) VERIFY_(STATUS) ! EXT from AERO_PROVIDER diff --git a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSradiation_GridComp/GEOSsolar_GridComp/GEOS_SolarGridComp.F90 b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSradiation_GridComp/GEOSsolar_GridComp/GEOS_SolarGridComp.F90 index 10dcb2538..02d07b194 100644 --- a/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSradiation_GridComp/GEOSsolar_GridComp/GEOS_SolarGridComp.F90 +++ b/GEOSagcm_GridComp/GEOSphysics_GridComp/GEOSradiation_GridComp/GEOSsolar_GridComp/GEOS_SolarGridComp.F90 @@ -2033,7 +2033,7 @@ subroutine RUN ( GC, IMPORT, EXPORT, CLOCK, RC ) VERIFY_(STATUS) ! execute the aero provider's optics method - call ESMF_MethodExecute(AERO, label="aerosol_optics", userRC=AS_STATUS, RC=STATUS) + call ESMF_MethodExecute(AERO, label="run_aerosol_optics", userRC=AS_STATUS, RC=STATUS) VERIFY_(AS_STATUS) VERIFY_(STATUS) @@ -2632,6 +2632,7 @@ subroutine SORADCORE(IM,JM,LM,NO_AERO,CURRTIME,LoadBalance,RC) if(SLICESimp(K) > 0) then ! Skip PREF if (NAMESimp(K)=="AERO") then +! if (NAMESimp(K)=="AERO_RAD") then _ASSERT(size(AEROSOL_EXT,3)==LM,'needs informative message') _ASSERT(size(AEROSOL_SSA,3)==LM,'needs informative message')