clm_out.f90

CLM集合卡曼滤波数据同化算法

       SUBROUTINE clm_out (lun_out ,lon ,lat,loc_yea,loc_day,loc_hou) 
!=======================================================================
!      Source file: clm_out.f90
! Original version: Yongjiu Dai, September 15, 1999
!=======================================================================

       Use clm2d_module
       Use clmtv_module
       Use clmtc_module

       IMPLICIT NONE

       integer, INTENT(in) :: & 
                lun_out,      &! logical unit number of output file 
                lon,          &! atm number of longitudes
                lat            ! atm number of latitudes
!      Local variables
       integer i, j
       real, dimension(lon,lat) :: &
                wind,         &! wind velocity (m/s)
                solar,        &! incident solar radiation (W/m2)
                prcp,         &! precipitation (mm/s)
                net_rad        ! net radiation (w/m2)
       real calday             ! calendar day + fraction (1.xx -> 365.xx)
       real tir
       integer::loc_yea,loc_day,loc_hou
! ----------------------------------------------------------------------
      loc_yea=loc_yea
      loc_day=loc_day
      loc_hou=loc_hou

       calday = float(jday) + (float(mcsec)-dtime)/86400.
       if((float(mcsec)-dtime) < 0.)then
           calday = float(jday-1) + (86400.+float(mcsec)-dtime)/86400.
       else
           calday = float(jday) + (float(mcsec)-dtime)/86400.
       endif

       do j = 1, lat
       do i = 1, lon
       wind(i,j) = sqrt(ugcmxy(i,j)**2 + vgcmxy(i,j)**2)
       prcp(i,j) = prcxy(i,j) + prlxy(i,j)
       solar(i,j) = solsxy(i,j) + sollxy(i,j) + solsdxy(i,j) + solldxy(i,j)
       net_rad(i,j)= flwdsxy(i,j) - lwupxy(i,j) + sabvxy(i,j) + sabgxy(i,j)  
      ! plot:  solar(i,j), flwdsxy(i,j), lwupxy(i,j), sabvxy(i,j), sabgxy(i,j)
       print*,flwdsxy 
       enddo
       enddo

       write(lun_out)                &!

       calday                       ,&! [1] 
! ------------------------------------------------------
! Fluxes
       solar                        ,&! [2]  [W/m2]
       sabvxy                       ,&! [3]  [W/m2]
       sabgxy                       ,&! [4]  [W/m2]
       flwdsxy                      ,&! [5]  [W/m2]
       lwupxy                       ,&! [6]  [W/m2]

       shxy                         ,&! [7]  [w/m2] 
       fsenlxy                      ,&! [8]  [W/m2]
       fsengxy                      ,&! [9]  [W/m2]

       lhxy                         ,&! [10] [W/m2] 
       fevplxy                      ,&! [11] [mm/s]
       etrxy                        ,&! [12] [mm/s]
       fevpgxy                      ,&! [13] [mm/s]

       fgrndxy                      ,&! [14] [W/m2]
       tauxxy                       ,&! [15] (kg/m/s2)
       tauyxy                       ,&! [16] (kg/m/s2)
! ------------------------------------------------------
! water
       prcp                         ,&! [17] (mm h2o/s)
       rsurxy                       ,&! [18] (mm h2o/s)
       rnofxy                       ,&! [19] (mm h2o/s)
       ldewxy                       ,&! [20] [mm]
       scv2xy                       ,&! [21] (mm)
       snowdpxy                     ,&! [22] [meter]
! ------------------------------------------------------
! assimilation of CO2
       rstxy                        ,&! [23]
       assimxy                      ,&! [24]
       respcxy                      ,&! [25]
! ------------------------------------------------------
! temperatures and some model state variables
       tgcmxy                       ,&! [26] (K)
       tsxy                         ,&! [27] (kelvin)
       trefxy                       ,&! [28] (kelvin)
       tlsunxy                      ,&! [29] [K]
       tlshaxy                      ,&! [30]
       tgxy                         ,&! [31] [K]

       wind                         ,&! [32] (m/s)
       qgcmxy                       ,&! [33] (kg/kg)

       fsnoxy                       ,&! [34] -
       sigfxy                       ,&! [35] -

       soitemxy                     ,&! [36-45] [K]
       soiliqxy                     ,&! [46-55] [kg/m2]
       soiicexy                       ! [56-65] [kg/m2]
! ------------------------------------------------------
!       print*,'dadadada......'
!       print*,'lalalala......'  
       tir = (lwupxy(1,1)/5.6697E-8)**(0.25)
       write(8,250)soitemxy 
   250 format(i4,i7,i7,25f30.20)     
       END SUBROUTINE clm_out