TOPAZ Column physcis에서 fed 바둑판 error

[HCJung-jbnu]

[HCJung-jbnu]

• ~/TOPAZ/generic_tracer_utils.F90 : SUB. g_tracer_set_common

• mask 데이터를 활용해서 coast mask 데이터를 생성

• 그러므로 isc:iec, jsc:jec 외의 halo 포함 영역(isd:ied, jsd:jed)의 mask 자료가 필요함. !Determine the coast line. !In order to that grid_tmask must have the proper value on the data domain boundaries isd,ied,jsd,jed !so that we can decide if the coast line coinsides with a point on the compute domain boundary

  do j =jsc, jec ; do i = isc, iec if (g_tracer_com%grid_tmask(i,j,1) .gt. 0) then if (g_tracer_com%grid_tmask(i-1,j,1) .lt. 1 .or. g_tracer_com%grid_tmask(i,j-1,1) .lt. 1 .or. & g_tracer_com%grid_tmask(i+1,j,1) .lt. 1 .or. g_tracer_com%grid_tmask(i,j+1,1) .lt. 1) then !{ g_tracer_com%grid_mask_coast(i,j) = 1 endif endif enddo; enddo

[HCJung-jbnu]

• ~/Nonlinear/initial.F : SUB. initial
• isd:ied, jsd:jed의 mask 자료 받도록 수정

      !Determine the coast line.
      !In order to that grid_tmask must have the proper value on the data domain boundaries isd,ied,jsd,jed
      !so that we can decide if the coast line coinsides with a point on the compute domain boundary
      DO nk = 1, N(ng)
        TP(ng)%tmask(isd:ied,jsd:jed,nk)=GRID(ng)%rmask(isd:ied,jsd:jed) !rho mask
  #ifdef WET_DRY
        TP(ng)%tmask(isd:ied,jsd:jed,nk)=GRID(ng)%rmask_wet(isd:ied,jsd:jed) !rho mask
  #endif
      END DO
      TP(ng)%mbathy(is:ie,js:je)=N(ng) ! sigma coord
      TP(ng)%area(is:ie,js:je)=GRID(ng)%omn(is:ie,js:je) ! area(m2)
    END DO
    END DO
  #endif

• ~/TOPAZ/generic_TOPAZ_model.F90 : SUB. generic_TOPAZ_init

    !Determine the coast line.
    !In order to that grid_tmask must have the proper value on the data domain boundaries isd,ied,jsd,jed
    !so that we can decide if the coast line coinsides with a point on the compute domain boundary
    do i = isd, ied; do j = jsd, jed; do k = 1, Grids%nk
       INTP(ng)%grid_tmask(i,j,k)=TP(ng)%tmask(i,j,k)
    end do; end do; end do