Ignore:
Timestamp:
Apr 6, 2016 3:44:20 PM (5 years ago)
Author:
maronga
Message:

further modularization of land surface model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/header.f90

    r1809 r1817  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Moved land_surface_model header output to the respective module.
    2222!
    2323! Former revisions:
     
    258258       
    259259    USE kinds
    260   
     260 
    261261    USE land_surface_model_mod,                                                &
    262         ONLY:  conserve_water_content, land_surface, nzb_soil,                 &
    263                nzt_soil, root_fraction, soil_moisture, soil_temperature,       &
    264                soil_type, soil_type_name, veg_type, veg_type_name, zs
     262        ONLY: land_surface, lsm_header
    265263 
    266264    USE model_1d,                                                              &
     
    969967    ENDIF
    970968
    971 
    972     IF ( land_surface )  THEN
    973 
    974        temperatures = ''
    975        gradients    = '' ! use for humidity here
    976        coordinates  = '' ! use for height
    977        roots        = '' ! use for root fraction
    978        slices       = '' ! use for index
    979 
    980        i = 1
    981        DO i = nzb_soil, nzt_soil
    982           WRITE (coor_chr,'(F10.2,7X)') soil_temperature(i)
    983           temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr )
    984 
    985           WRITE (coor_chr,'(F10.2,7X)') soil_moisture(i)
    986           gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr )
    987 
    988           WRITE (coor_chr,'(F10.2,7X)')  - zs(i)
    989           coordinates = TRIM( coordinates ) // ' '  // TRIM( coor_chr )
    990 
    991           WRITE (coor_chr,'(F10.2,7X)')  root_fraction(i)
    992           roots = TRIM( roots ) // ' '  // TRIM( coor_chr )
    993 
    994           WRITE (coor_chr,'(I10,7X)')  i
    995           slices = TRIM( slices ) // ' '  // TRIM( coor_chr )
    996 
    997 
    998        ENDDO
    999 
    1000 !
    1001 !--    Write land surface model header
    1002        WRITE( io, 419 )
    1003        IF ( conserve_water_content )  THEN
    1004           WRITE( io, 440 )
    1005        ELSE
    1006           WRITE( io, 441 )
    1007        ENDIF
    1008 
    1009        WRITE( io, 438 ) TRIM( veg_type_name(veg_type) ),                       &
    1010                         TRIM (soil_type_name(soil_type) )
    1011        WRITE( io, 439 ) TRIM( coordinates ), TRIM( temperatures ),             &
    1012                         TRIM( gradients ), TRIM( roots ), TRIM( slices )
    1013 
    1014 
    1015     ENDIF
     969    IF ( land_surface )  CALL lsm_header ( io )
    1016970
    1017971    IF ( radiation )  THEN
     
    23082262418 FORMAT (/'    Day of the year at model start :   day_init      =     ',I3 &
    23092263            /'    UTC time at model start        :   time_utc_init = ',F7.1' s')
    2310 419 FORMAT (//' Land surface model information:'/ &
    2311               ' ------------------------------'/)
    23122264420 FORMAT (/'    Characteristic levels of the initial temperature profile:'// &
    23132265            '       Height:        ',A,'  m'/ &
     
    23682320                       '[0,1000] cm**2/s**3')
    23692321437 FORMAT ('    Droplet collision is switched off')
    2370 438 FORMAT (' --> Land surface type  : ',A,/ &
    2371             ' --> Soil porosity type : ',A)
    2372 439 FORMAT (/'    Initial soil temperature and moisture profile:'// &
    2373             '       Height:        ',A,'  m'/ &
    2374             '       Temperature:   ',A,'  K'/ &
    2375             '       Moisture:      ',A,'  m**3/m**3'/ &
    2376             '       Root fraction: ',A,'  '/ &
    2377             '       Gridpoint:     ',A)
    2378 440 FORMAT (' --> Soil bottom is closed (water content is conserved, default)')
    2379 441 FORMAT (' --> Soil bottom is open (water content is not conserved)')
    23802322444 FORMAT (//' Radiation model information:'/                                 &
    23812323              ' ----------------------------'/)
Note: See TracChangeset for help on using the changeset viewer.