Add TSOIL1

for the new soil NOx extension option
noaa-oar-arl · Nov 25, 2024 · f9ef39c · f9ef39c
1 parent 6dee43f
commit f9ef39c
Showing 1 changed file with 16 additions and 1 deletion.
diff --git a/src/nuopc/cap.F90 b/src/nuopc/cap.F90
@@ -1710,7 +1710,7 @@ subroutine hco_ext_set_fields ( HcoState, ExtState, RC )
       end if
     end if
 
-    !%%%%% Air and skin temperature %%%%%
+    !%%%%% Air temperature %%%%%
     if ( ExtState%T2M%DoUse ) then
       Name = 'T2M'
       call ExtDat_Set( HcoState,     ExtState%T2M,                          &
@@ -1724,6 +1724,7 @@ subroutine hco_ext_set_fields ( HcoState, ExtState, RC )
       end if
     end if
 
+    !%%%%% Skin temperature %%%%%
     if ( ExtState%TSKIN%DoUse ) then
       Name = 'TS'
       call ExtDat_Set( HcoState,     ExtState%TSKIN,                        &
@@ -1737,6 +1738,20 @@ subroutine hco_ext_set_fields ( HcoState, ExtState, RC )
       end if
     end if
 
+    !%%%%% Soil temperature %%%%%
+    IF ( ExtState%TSOIL1%DoUse ) THEN
+      Name = 'TSOIL1'
+      CALL ExtDat_Set( HcoState,     ExtState%TSOIL1,                       &
+      trim( Name ), RC,       FIRST=FIRST                 )
+      if ( RC /= HCO_SUCCESS ) then
+        ErrMsg = 'Could not find quantity "' // trim( Name )            // &
+          '" for the HEMCO standalone simulation!'
+        call HCO_Error( HcoConfig%Err, ErrMsg, RC, ThisLoc )
+        call HCO_Leave( HcoState%Config%Err, RC )
+        return
+      end if
+    end if
+
     !%%%%% Soil moisture %%%%%
     if ( ExtState%GWETROOT%DoUse ) then
       Name = 'GWETROOT'