◆ vg6d_c2a()

subroutine, public volgrid6d_class::vg6d_c2a ( type(volgrid6d), dimension(:), intent(inout) this )

Convert grids type C to type A.

Use this to interpolate data from grid type C to grid type A Grids type are defined by Arakawa.

We need to find these types of area in a vg6d array T area of points with temterature etc. U area of points with u components of winds V area of points with v components of winds

this method works if it finds two volumes: 1) volume U and volume V: compute H and traslate on it 2) volume U/V and volume H : translate U/V on H three volumes: translate U and V on H

try to work well on more datasets at once

Parametri

[in,out] this vettor of volumes volgrid6d to elaborate

Definizione alla linea 3205 del file volgrid6d_class.F90.

     step_lat=(ymax-ymin)/dble(this%griddim%dim%ny-1)
     ymin=ymin-step_lat/2.d0
     ymax=ymax-step_lat/2.d0
     call set_val(this%griddim, ymin=ymin, ymax=ymax)
 ! improve grid definition precision
     CALL griddim_setsteps(this%griddim)
     
   case default
 
     call l4f_category_log(this%category,l4f_fatal,"C grid type not known")
     call raise_fatal_error ()
 
   end select
 
 end if
 
 
 call griddim_unproj(this%griddim)
 
 
 end subroutine vg6d_c2a_grid
 
 ! Convert C grid to A grid
 subroutine vg6d_c2a_mat(this,cgrid)
 
 type(volgrid6d),intent(inout) :: this ! object containing fields to be translated
 integer,intent(in) :: cgrid ! in C grid (Arakawa) we have 0=T,1=U,2=V points
 
 INTEGER :: i, j, k, iv, stallo
 REAL,ALLOCATABLE :: tmp_arr(:,:)
 REAL,POINTER :: voldatiuv(:,:)
 
 
 IF (cgrid == 0) RETURN ! nothing to do
 IF (cgrid /= 1 .AND. cgrid /= 2) THEN ! wrong cgrid
   CALL l4f_category_log(this%category,l4f_fatal,"C grid type "// &
    trim(to_char(cgrid))//" not known")
   CALL raise_error()
   RETURN
 ENDIF
 
 ! Temporary workspace
 ALLOCATE(tmp_arr(this%griddim%dim%nx, this%griddim%dim%ny),stat=stallo)
 if (stallo /=0)then
   call l4f_log(l4f_fatal,"allocating memory")
   call raise_fatal_error()
 end if
 
 ! allocate once for speed
 IF (.NOT.ASSOCIATED(this%voldati)) THEN
   ALLOCATE(voldatiuv(this%griddim%dim%nx, this%griddim%dim%ny), stat=stallo)
   IF (stallo /= 0) THEN
     CALL l4f_log(l4f_fatal,"allocating memory")
     CALL raise_fatal_error()
   ENDIF
 ENDIF
 
 IF (cgrid == 1) THEN ! U points to H points
   DO iv = 1, SIZE(this%var)
     DO k = 1, SIZE(this%timerange)
       DO j = 1, SIZE(this%time)
         DO i = 1, SIZE(this%level)
           tmp_arr(:,:) = rmiss
           CALL volgrid_get_vol_2d(this, i, j, k, iv, voldatiuv)
 
 ! West boundary extrapolation
           WHERE(voldatiuv(1,:) /= rmiss .AND. voldatiuv(2,:) /= rmiss)
             tmp_arr(1,:) = voldatiuv(1,:) - (voldatiuv(2,:) - voldatiuv(1,:)) / 2.
           END WHERE
 
 ! Rest of the field
           WHERE(voldatiuv(1:this%griddim%dim%nx-1,:) /= rmiss .AND. &
            voldatiuv(2:this%griddim%dim%nx,:) /= rmiss)
             tmp_arr(2:this%griddim%dim%nx,:) = &
              (voldatiuv(1:this%griddim%dim%nx-1,:) + &
              voldatiuv(2:this%griddim%dim%nx,:)) / 2.
           END WHERE
 
           voldatiuv(:,:) = tmp_arr
           CALL volgrid_set_vol_2d(this, i, j, k, iv, voldatiuv)
         ENDDO
       ENDDO
     ENDDO
   ENDDO
 
 ELSE IF (cgrid == 2) THEN ! V points to H points
   DO iv = 1, SIZE(this%var)
     DO k = 1, SIZE(this%timerange)
       DO j = 1, SIZE(this%time)
         DO i = 1, SIZE(this%level)
           tmp_arr(:,:) = rmiss
           CALL volgrid_get_vol_2d(this, i, j, k, iv, voldatiuv)
 
 ! South boundary extrapolation
           WHERE(voldatiuv(:,1) /= rmiss .AND. voldatiuv(:,2) /= rmiss)
             tmp_arr(:,1) = voldatiuv(:,1) - (voldatiuv(:,2) - voldatiuv(:,1)) / 2.
           END WHERE
 
 ! Rest of the field
           WHERE(voldatiuv(:,1:this%griddim%dim%ny-1) /= rmiss .AND. &
            voldatiuv(:,2:this%griddim%dim%ny) /= rmiss)
             tmp_arr(:,2:this%griddim%dim%ny) = &
              (voldatiuv(:,1:this%griddim%dim%ny-1) + &
              voldatiuv(:,2:this%griddim%dim%ny)) / 2.
           END WHERE
 
           voldatiuv(:,:) = tmp_arr
           CALL volgrid_set_vol_2d(this, i, j, k, iv, voldatiuv)
         ENDDO
       ENDDO
     ENDDO
   ENDDO
 ENDIF ! tertium non datur
 
 IF (.NOT.ASSOCIATED(this%voldati)) THEN
   DEALLOCATE(voldatiuv)
 ENDIF
 DEALLOCATE (tmp_arr)
 
 end subroutine vg6d_c2a_mat
 
 