Ignore:
Timestamp:
Oct 21, 2012 1:03:21 PM (9 years ago)
Author:
letzel
Message:
  • mask locations determined based on scalar positions (init_masks)
  • save memory by not allocating pt_2 in case of neutral = .T. (init_3d_model, swap_timelevel)
  • minor reformatting (check_for_restart)
File:
1 edited

Legend:

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

    r1031 r1032  
    44! Current revisions:
    55! -----------------
     6! hotfix: mask locations determined based on scalar positions
    67! netCDF4 without parallel file support implemented
    78!
     
    464465             count = count + 1
    465466             IF ( dim == 1 .OR. dim == 2 )  THEN
    466                 m = NINT( mask(mid,dim,count) * mask_scale(dim) * ddxyz )
     467                m = NINT( mask(mid,dim,count) * mask_scale(dim) * ddxyz - 0.5 )
    467468             ELSEIF ( dim == 3 )  THEN
    468469                ind_array =  &
     
    515516                CALL message( 'init_masks', 'PA0332', 1, 2, 0, 6, 0 )
    516517             ENDIF
    517              loop_begin  = NINT( mask_loop(mid,dim,1) * mask_scale(dim) * ddxyz )
    518              loop_end    = NINT( mask_loop(mid,dim,2) * mask_scale(dim) * ddxyz )
    519              loop_stride = NINT( mask_loop(mid,dim,3) * mask_scale(dim) * ddxyz )
     518             loop_begin  = NINT( mask_loop(mid,dim,1) * mask_scale(dim) &
     519                  * ddxyz - 0.5 )
     520             loop_end    = NINT( mask_loop(mid,dim,2) * mask_scale(dim) &
     521                  * ddxyz - 0.5 )
     522             loop_stride = NINT( mask_loop(mid,dim,3) * mask_scale(dim) &
     523                  * ddxyz )
    520524          ELSEIF ( dim == 3 )  THEN
    521525             IF ( mask_loop(mid,dim,2) < 0.0 )  THEN
Note: See TracChangeset for help on using the changeset viewer.