Changeset 138 for palm/trunk/SOURCE/header.f90

Timestamp:

Nov 28, 2007 10:03:58 AM (16 years ago)

Author:

letzel

Message:

Plant canopy model of Watanabe (2004,BLM 112,307-341) added.

File:

• palm/trunk/SOURCE/header.f90 (modified) (3 diffs)

palm/trunk/SOURCE/header.f90

@@ -78 +78 @@
     CHARACTER (LEN=70) ::  char1, char2, coordinates, gradients, dopr_chr, &
                            do2d_xy, do2d_xz, do2d_yz, do3d_chr, &
-                           run_classification, slices, temperatures, &
+                           run_classification, slices, temperatures, learde, &
                            ugcomponent, vgcomponent
     CHARACTER (LEN=85) ::  roben, runten
@@ -364 +364 @@
 
     END SELECT
+
+    IF ( plant_canopy ) THEN
+
+       WRITE ( io, 280 ) canopy_mode, pch_index, drag_coefficient
+
+!
+!--    Leaf area density profile
+!--    Building output strings, starting with surface value
+       WRITE ( learde, '(F6.2)' )  lad_surface
+       gradients = '------'
+       slices = '     0'
+       coordinates = '   0.0'
+       i = 1
+       DO  WHILE ( lad_vertical_gradient_level_ind(i) /= -9999 )
+
+          WRITE (coor_chr,'(F7.2)')  lad(lad_vertical_gradient_level_ind(i))
+          learde = TRIM( learde ) // ' ' // TRIM( coor_chr )
+
+          WRITE (coor_chr,'(F7.2)')  lad_vertical_gradient(i)
+          gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
+
+          WRITE (coor_chr,'(I7)')  lad_vertical_gradient_level_ind(i)
+          slices = TRIM( slices ) // ' ' // TRIM( coor_chr )
+
+          WRITE (coor_chr,'(F7.1)')  lad_vertical_gradient_level(i)
+          coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
+
+          i = i + 1
+       ENDDO
+
+       WRITE ( io, 281 )  TRIM( coordinates ), TRIM( learde ), &
+                          TRIM( gradients ), TRIM( slices )
+
+    ENDIF
 
 !
@@ -1333 +1367 @@
               ' Horizontal index bounds (l/r/s/n): ',I4,' / ',I4,' / ',I4, &
                 ' / ',I4)
+280 FORMAT (//' Vegetation canopy (drag) model:'/ &
+              ' ------------------------------'// &
+              ' Canopy mode: ', A / &
+              ' Canopy top: ',I4 / &
+              ' Leaf drag coefficient: ',F6.2 /)
+281 FORMAT (/ ' Characteristic levels of the leaf area density:'// &
+              ' Height:              ',A,'  m'/ &
+              ' Leaf area density:   ',A,'  m**2/m**3'/ &
+              ' Gradient:            ',A,'  m**2/m**4'/ &
+              ' Gridpoint:           ',A)
+               
 300 FORMAT (//' Boundary conditions:'/ &
              ' -------------------'// &