Coupling support for GEOS-ESM

columnphysics/icepack_flux.F90

@@ -43,7 +43,9 @@ subroutine merge_fluxes (aicen,                &
                                fhocnn,   fswthrun,   &
                                fswthrun_vdr, fswthrun_vdf,&
                                fswthrun_idr, fswthrun_idf,&
-                               strairxT, strairyT,   &
+                               fswthrun_uvrdr, fswthrun_uvrdf,&
+                               fswthrun_pardr, fswthrun_pardf,&
+                               strairxT, strairyT,   &  
                                Cdn_atm_ratio,        &
                                fsurf,    fcondtop,   &
                                fcondbot,             &
@@ -56,6 +58,8 @@ subroutine merge_fluxes (aicen,                &
                                fhocn,    fswthru,    &
                                fswthru_vdr, fswthru_vdf,&
                                fswthru_idr, fswthru_idf,&
+                               fswthru_uvrdr, fswthru_uvrdf,&
+                               fswthru_pardr, fswthru_pardf,&
                                melttn, meltsn, meltbn, congeln, snoicen, &
                                meltt,  melts,        &
                                meltb,  dsnow, dsnown,&
@@ -89,6 +93,10 @@ subroutine merge_fluxes (aicen,                &
           fsaltn  , & ! salt flux to ocean              (kg/m2/s)
           fhocnn  , & ! actual ocn/ice heat flx         (W/m**2)
           fswthrun, & ! sw radiation through ice bot    (W/m**2)
+          fswthrun_uvrdr, & ! vis dir sw radiation through ice bot    (W/m**2)
+          fswthrun_uvrdf, & ! vis dir sw radiation through ice bot    (W/m**2)
+          fswthrun_pardr, & ! vis dir sw radiation through ice bot    (W/m**2)
+          fswthrun_pardf, & ! vis dir sw radiation through ice bot    (W/m**2)
           melttn  , & ! top ice melt                    (m)
           meltbn  , & ! bottom ice melt                 (m)
           meltsn  , & ! snow melt                       (m)
@@ -139,6 +147,12 @@ subroutine merge_fluxes (aicen,                &
           fswthru_idf     ! nir dif sw radiation through ice bot    (W/m**2)
 
       real (kind=dbl_kind), intent(inout), optional :: &
+          fswthru_uvrdr , & ! vis dir sw radiation through ice bot    (W/m**2)
+          fswthru_uvrdf , & ! vis dif sw radiation through ice bot    (W/m**2)
+          fswthru_pardr , & ! nir dir sw radiation through ice bot    (W/m**2)
+          fswthru_pardf     ! nir dif sw radiation through ice bot    (W/m**2)
+
+      real (kind=dbl_kind), optional, intent(inout):: &
           dsnow,    & ! change in snow depth            (m)
           Uref        ! air speed reference level       (m/s)
 
@@ -212,6 +226,15 @@ subroutine merge_fluxes (aicen,                &
       if (present(fswthru_idf)) &
          fswthru_idf   = fswthru_idf   + fswthrun_idf  * aicen
 
+      if (present(fswthru_uvrdr)) &
+         fswthru_uvrdr   = fswthru_uvrdr   + fswthrun_uvrdr  * aicen
+      if (present(fswthru_uvrdf)) &
+         fswthru_uvrdf   = fswthru_uvrdf   + fswthrun_uvrdf  * aicen
+      if (present(fswthru_pardr)) &
+         fswthru_pardr   = fswthru_pardr   + fswthrun_pardr  * aicen
+      if (present(fswthru_pardf)) &
+         fswthru_pardf   = fswthru_pardf   + fswthrun_pardf  * aicen
+
       ! ice/snow thickness
 
       meltt     = meltt     + melttn    * aicen