Home


About us
PALM
PALM-4U
Installation
Release notes
Documentation
Getting Help
Gallery
Developers area

Ignore:
Timestamp:
Apr 30, 2014 12:15:41 PM (10 years ago)
Author:
boeske
Message:

minor changes in profile data output of lsf tendencies, variables renamed

File:
1 edited

Legend:

Unmodified
Added
Removed

  • palm/trunk/SOURCE/nudging.f90

    r1381 r1382  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Changed the weighting factor that is used in the summation of nudging
     23! tendencies for profile data output from weight_pres to weight_substep,
     24! added Neumann boundary conditions for profile data output of nudging terms at
     25! nzt+1
    2326!
    2427! Former revisions:
    … …  
    299302
    300303       USE statistics,                                                         &
    301            ONLY:  hom, sums_ls_l, weight_pres
     304           ONLY:  hom, sums_ls_l, weight_substep
    302305
    303306       IMPLICIT NONE
    … …  
    333336             DO  i = nxl, nxr
    334337                DO  j = nys, nyn
     338
    335339                   DO  k = nzb_u_inner(j,i)+1, nzt
    336340
    … …  
    341345
    342346                      sums_ls_l(k,6) = sums_ls_l(k,6) + tmp_tend *             &
    343                                        weight_pres(intermediate_timestep_count)
     347                                     weight_substep(intermediate_timestep_count)
    344348                   ENDDO
     349                 
     350                   sums_ls_l(nzt+1,6) = sums_ls_l(nzt,6)
     351 
    345352                ENDDO
    346353            ENDDO
    … …  
    350357             DO  i = nxl, nxr
    351358                DO  j = nys, nyn
     359
    352360                   DO  k = nzb_u_inner(j,i)+1, nzt
    353361
    … …  
    358366
    359367                      sums_ls_l(k,7) = sums_ls_l(k,7) + tmp_tend *             &
    360                                        weight_pres(intermediate_timestep_count)
     368                                     weight_substep(intermediate_timestep_count)
    361369                   ENDDO
     370                 
     371                   sums_ls_l(nzt+1,7) = sums_ls_l(nzt,7)
     372
    362373                ENDDO
    363374            ENDDO
    … …  
    367378             DO  i = nxl, nxr
    368379                DO  j = nys, nyn
     380
    369381                   DO  k = nzb_u_inner(j,i)+1, nzt
    370382
    … …  
    375387
    376388                      sums_ls_l(k,4) = sums_ls_l(k,4) + tmp_tend *             &
    377                                        weight_pres(intermediate_timestep_count)
     389                                     weight_substep(intermediate_timestep_count)
    378390                   ENDDO
     391
     392                   sums_ls_l(nzt+1,4) = sums_ls_l(nzt,4)
     393
    379394                ENDDO
    380395            ENDDO
    … …  
    384399             DO  i = nxl, nxr
    385400                DO  j = nys, nyn
     401
    386402                   DO  k = nzb_u_inner(j,i)+1, nzt
    387403
    … …  
    392408
    393409                      sums_ls_l(k,5) = sums_ls_l(k,5) + tmp_tend *             &
    394                                        weight_pres(intermediate_timestep_count)
     410                                     weight_substep(intermediate_timestep_count)
    395411                   ENDDO
     412                 
     413                   sums_ls_l(nzt+1,5) = sums_ls_l(nzt,5)
     414
    396415                ENDDO
    397416            ENDDO
    … …  
    425444
    426445       USE statistics,                                                         &
    427            ONLY:  hom, sums_ls_l, weight_pres
     446           ONLY:  hom, sums_ls_l, weight_substep
    428447
    429448       IMPLICIT NONE
    … …  
    466485
    467486                sums_ls_l(k,6) = sums_ls_l(k,6) + tmp_tend                     &
    468                                  * weight_pres(intermediate_timestep_count)
     487                                 * weight_substep(intermediate_timestep_count)
    469488             ENDDO
     489
     490             sums_ls_l(nzt+1,6) = sums_ls_l(nzt,6)
    470491
    471492          CASE ( 'v' )
    … …  
    479500
    480501                sums_ls_l(k,7) = sums_ls_l(k,7) + tmp_tend                     &
    481                                  * weight_pres(intermediate_timestep_count)
     502                                 * weight_substep(intermediate_timestep_count)
    482503             ENDDO
    483504
     505             sums_ls_l(nzt+1,7) = sums_ls_l(nzt,7)
    484506
    485507          CASE ( 'pt' )
    … …  
    493515
    494516                sums_ls_l(k,4) = sums_ls_l(k,4) + tmp_tend                     &
    495                                  * weight_pres(intermediate_timestep_count)
     517                                 * weight_substep(intermediate_timestep_count)
    496518             ENDDO
     519
     520             sums_ls_l(nzt+1,4) = sums_ls_l(nzt,4)
    497521
    498522
    … …  
    507531
    508532                sums_ls_l(k,5) = sums_ls_l(k,5) + tmp_tend                     &
    509                                  * weight_pres(intermediate_timestep_count)
     533                                 * weight_substep(intermediate_timestep_count)
    510534             ENDDO
     535
     536             sums_ls_l(nzt+1,5) = sums_ls_l(nzt,5)
    511537
    512538          CASE DEFAULT
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |