在進行數值實驗時,常常會藉由修改地形高度來檢驗一個天氣系統的變化,大概每一篇颱風模擬的Paper都會出現一個關於地形敏感度實驗的章節XD,以下內容出自我寫給學弟妹的MM5使用手冊:
MM5/TERRAIN的地形製作流程:
! =========================================================================== ! = MM5 Terrain modification, ! = Insert me between "continue" (CRTER.171) and "return" (CRTER.174) ! = by Dennis and Cypresslin - http://cypresslin.web.fc2.com/ ! =========================================================================== | |
iActivate = 1 factor = 0.0 | ! 1 - yes; 0 - no ! multiplier, 0.0 - no terrain; 50.0 - 50% terrain ... etc |
if ( iActivate .eq. 1 ) then do j = 1,JX do i = 1,IX ! Specify the area you want to modify if ( XLAT(i,j) .gt. 21.5 .and. XLAT(i,j) .lt. 25.5 .and. & XLON(i,j) .gt. 119.95 .and. XLON(i,j) .lt. 122.02) then HT(i,j) = HT(i,j) * factor endif enddo enddo endif |
! =========================================================================== ! = MM5 Land-Water mask modification, ! = Insert me between "ENDIF" (CRLWMSK.175) and "ENDDO" (CRLWMSK.176) ! = by Cypresslin - http://cypresslin.web.fc2.com/ ! =========================================================================== iActivate = 1 if ( iActivate .eq. 1 ) then ! Specify the area you want to modify if ( XLAT(i,j) .gt. 21.5 .and. XLAT(i,j) .lt. 25.5 .and. & XLON(i,j) .gt. 119.95 .and. XLON(i,j) .lt. 122.02) then FORCETER(i,j) = 0.001 endif endif |
製作出TERRAIN_DOMAIN#之後,不要忘記先看看他的TER.PLT,或者是用MM5toGrADS之類的轉檔程式把TERRAIN_DOMAIN#轉出來修改後的檢查正確性。
Special thanks to: Dennis (NCU)
<-Memo