Changeset 952 for trunk/LMDZ.GENERIC/libf

Timestamp:

May 1, 2013, 2:01:43 AM (12 years ago)

Author:

aslmd

Message:

LMDZ.GENERIC. added the possibility to remove surface (nosurf option) and to add an internal heat flux (intheat value, default is 0.)

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callcorrk.F90 (modified) (1 diff)
• callkeys.h (modified) (3 diffs)
• inifis.F (modified) (1 diff)
• physiq.F90 (modified) (3 diffs)
• turbdiff.F90 (modified) (1 diff)
• vdifc.F (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

@@ -450 +450 @@
          albv=albedo(ig)        ! shortwave 
 
-      if(noradsurf.and.(albv.gt.0.0))then
+      if(nosurf.and.(albv.gt.0.0))then
          print*,'For open lower boundary in callcorrk must'
          print*,'have surface albedo set to zero!'

trunk/LMDZ.GENERIC/libf/phystd/callkeys.h

@@ -11 +11 @@
      &   , Nmix_co2, radfixed, dusttau                                  &
      &   , meanOLR, specOLR                                             &
-     &   , kastprof, noradsurf, Tstrat                                  &
+     &   , kastprof, Tstrat                                             &
+     &   , nosurf, intheat                                              &     
      &   , newtonian, tau_relax, testradtimes                           &
      &   , check_cpp_match, force_cpp                                   &
@@ -62 +63 @@
       logical sourceevol
       logical CLFvarying
-      logical noradsurf
+      logical nosurf
 
       integer iddist
@@ -87 +88 @@
       real cloudlvl
       real icetstep
-
+      real intheat

trunk/LMDZ.GENERIC/libf/phystd/inifis.F

@@ -307 +307 @@
 
          write(*,*)"Remove lower boundary?"
-         noradsurf=.false.
-         call getin("noradsurf",noradsurf)
-         write(*,*)" noradsurf = ",noradsurf
+         nosurf=.false.
+         call getin("nosurf",nosurf)
+         write(*,*)" nosurf = ",nosurf
 
 ! Tests of incompatibility:
-         if (noradsurf.and.callsoil) then
-           print*,'noradsurf not compatible with soil scheme!'
+         if (nosurf.and.callsoil) then
+           print*,'nosurf not compatible with soil scheme!'
+           print*,'... got to make a choice!'
            call abort
          endif
-         !if (noradsurf.and.calldifv) then
-         !  print*,'noradsurf not compatible with a boundary layer!'
-         !  call abort
-         !endif
+
+         write(*,*)"Add an internal heat flux?",
+     .             "... matters only if callsoil=F"
+         intheat=0.
+         call getin("intheat",intheat)
+         write(*,*)" intheat = ",intheat
 
          write(*,*)"Use Newtonian cooling for radiative transfer?"

trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

@@ -520 +520 @@
             print*,'WARNING! Thermal conduction in the soil turned off'
             capcal(:)=1.e6
-            fluxgrd(:)=0.
-            if(noradsurf)then
-                  fluxgrd(:)=10.0
-            endif
-            print*,'Flux from ground = ',fluxgrd,' W m^-2'
+            fluxgrd(:)=intheat
+            print*,'Flux from ground = ',intheat,' W m^-2'
          endif
          icount=1
@@ -760 +757 @@
               fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)+fluxsurf_sw(1:ngrid)*(1.-albedo(1:ngrid))
 
-              if(noradsurf)then ! no lower surface; SW flux just disappears
-                 GSR = SUM(fluxsurf_sw(1:ngrid)*area(1:ngrid))/totarea
-                 fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)
-                 print*,'SW lost in deep atmosphere = ',GSR,' W m^-2'
-              endif
+              !if(noradsurf)then ! no lower surface; SW flux just disappears
+              !   GSR = SUM(fluxsurf_sw(1:ngrid)*area(1:ngrid))/totarea
+              !   fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)
+              !   print*,'SW lost in deep atmosphere = ',GSR,' W m^-2'
+              !endif
 
 !             Net atmospheric radiative heating rate (K.s-1)
@@ -1324 +1321 @@
 !   9. Surface and sub-surface soil temperature
 !-----------------------------------------------------------------------
-
 
 !     Increment surface temperature

trunk/LMDZ.GENERIC/libf/phystd/turbdiff.F90

@@ -226 +226 @@
          zcdv_true(ig) = cd0
          zcdh_true(ig) = cd0
-!        zcdv_true(ig)=0.!JL12 uncomment to disable atm/surface momentum flux
-!        zcdh_true(ig)=0.!JL12 uncomment to disable sensible heat flux
+       if(nosurf)then
+         zcdv_true(ig)=0.D+0 !JL12 disable atm/surface momentum flux
+         zcdh_true(ig)=0.D+0 !JL12 disable sensible heat flux 
+       endif
       ENDDO
 

trunk/LMDZ.GENERIC/libf/phystd/vdifc.F

@@ -234 +234 @@
          zcdv_true(ig) = cd0
          zcdh_true(ig) = cd0
+         if (nosurf) then
+             zcdv_true(ig) = 0.   !! disable sensible momentum flux
+             zcdh_true(ig) = 0.   !! disable sensible heat flux
+         endif
       ENDDO