Changeset 3526 for trunk/LMDZ.COMMON/libf/evolution/NS_soilthprop_mars.F90

Timestamp:

Nov 19, 2024, 5:54:18 PM (4 days ago)

Author:

jbclement

Message:

PEM:
Removing unused or redundant Norbert Schorghofer's subroutines (follow-up of r3493) + cleaning and some modifications of related subroutines.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/NS_soilthprop_mars.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/evolution/NS_soilthprop_mars.F90

@@ -76 +76 @@
   
 end subroutine soilthprop
-
-
-      
-subroutine smartgrid(nz,z,zdepth,thIn,rhoc,porosity,ti,rhocv,layertype,icefrac)
-!***********************************************************************
-! smartgrid: returns intelligently spaced grid and appropriate 
-!            values of thermal inertia ti and rho*c in icy layer
-!                  
-!     INPUTS: 
-!             nz = number of grid points
-!             z = grid spacing for dry regolith 
-!                 (will be partially overwritten)
-!             zdepth = depth where ice table starts
-!                      negative values indicate no ice
-!             rhoc = heat capacity per volume of dry regolith [J/m^3]
-!             thIn = thermal inertia of dry regolith [SI-units]
-!             porosity = void space / total volume
-!             layertypes are explained below  
-!             icefrac = fraction of ice in icelayer
-!
-!     OUTPUTS: z = grid coordinates
-!              vectors ti and rhocv
-!***********************************************************************
-  implicit none
-  integer, intent(IN) :: nz, layertype
-  real(8), intent(IN) :: zdepth, thIn, rhoc, porosity, icefrac
-  real(8), intent(INOUT) :: z(nz)
-  real(8), intent(OUT) :: ti(nz), rhocv(nz)
-  integer j, b
-  real(8) stretch, newrhoc, newti
-  real(8), parameter :: NULL=0.
-  
-  if (zdepth>0 .and. zdepth<z(nz)) then
-
-     select case (layertype)
-     case (1)  ! interstitial ice
-        call soilthprop(porosity,1.d0,rhoc,thIn,1,newrhoc,newti,NULL)
-     case (2)  ! mostly ice (massive ice)
-        call soilthprop(porosity,NULL,rhoc,thIn,2,newrhoc,newti,icefrac)
-     case (3)  ! all ice (pure ice)
-        call soilthprop(NULL,NULL,NULL,NULL,3,newrhoc,newti,NULL)
-     case (4)  ! ice + rock + nothing else (ice-cemented soil)
-        call soilthprop(porosity,NULL,rhoc,NULL,4,newrhoc,newti,NULL)
-     case default
-        error stop 'invalid layer type'
-     end select
-
-     ! Thermal skin depth is proportional to sqrt(kappa)
-     ! thermal diffusivity kappa = k/(rho*c) = I^2/(rho*c)**2
-     stretch = (newti/thIn)*(rhoc/newrhoc) ! ratio of sqrt(thermal diffusivity)
-     
-     b = 1
-     do j=1,nz
-        if (z(j)<=zdepth) then 
-           b = j+1
-           rhocv(j) = rhoc
-           ti(j) = thIn
-        else
-           rhocv(j) = newrhoc
-           ti(j) = newti
-        endif
-        ! print *,j,z(j),ti(j),rhocv(j)
-     enddo
-     do j=b+1,nz
-        z(j) = z(b) + stretch*(z(j)-z(b))
-     enddo
-     
-     ! print *,'zmax=',z(nz),' b=',b,' stretch=',stretch
-     ! print *,'depth at b-1, b ',z(b-1),z(b)
-     ! print *,'ti1=',thIn,' ti2=',newti
-     ! print *,'rhoc1=',rhoc,' rhoc2=',newrhoc 
-  endif
-  
-end subroutine smartgrid
-