fortran+grads
- 格式：docx
- 大小：29.25 KB
- 文档页数：30
fortran+grads
WRF
Could not chdir to home directory /home/walilai: Stale NFS
file handle
/usr/X11R6/bin/xauth: error in locking authority file
/home/walilai/.Xauthority
Bug
INF ! 1/0=無限大(可用程式寫寫看)
NaN ! 1/0=無限大
read(unit,*) ! 注意unit數字使用限制
Segmentation fault!矩陣沒定義到的空間(I or j…)或從第6行開始寫(須從7)
do
i=i+1 ! 寫在這會出現Segmentation fault
station(i)=' ' ! 可能是if沒被定義到I,station
read(10,21,end=31) st,lon,lat,mm,dd,hh,rain
if(mm.eq.7.and.dd.eq.13.and.hh.eq.24)then
i=i+1 !OK
station(i)=' ' !OK
do j=1,i
do
i=i+1
c station(i)=' ' !若沒c掉，會出現Segmentation fault
read(11,110,end=97) st,yy,mm,dd,hh,c,rain
do_ud: off end of record:
k=0 !a.out OK!!
c do j=1,a !c處出現do_ud: off en
d of record:
c k=j-1 do i=1,a write(11,rec=3*k+1)station(i)
write(11,rec=3*k+2)st_lon(i)
write(11,rec=3*k+3)st_lat(i)
k=k+1
enddo
open: null file name!沒有讀到檔案，或是原本要輸入end才可結束程式的地方案到其他按鍵
list in: end of file!檔案讀不完全，可能status=unknown
open: 'new' file exists:status!因檔案存在，不可用new
do_ud: end of file!字串字數沒寫對
apparent state: unit 10 named f060200_112.5_35.dat
fmt: read unexpected character !找以下檔案會發現錯誤
apparent state: unit 11 named c-cwb200807.txt
last format: (A6,I4,I2,I2,I2,A47,F5.1)!所讀數字出現文字T(雨跡)
!format格式讀錯
!可能計算出現nan而失敗
invalid integer: read unexpected character
apparent state: unit 87 named autop2010-07.txt
last format: list io
lately reading sequential formatted external IO
Aborted
! read(87,100,end=920)st,yy,mm,dd,hh,c,rain
將格式由*改為100就可
fmt: end of file
apparent state: unit 92 named remove_2.txt
last format: (A6,1X,F6.1,1X,F5.1,2X,F7.1,2X,F7.1)
lately reading sequential formatted external IO
Aborted
!沒有close(??)或是無線迴圈的end=900，或是開檔時只有11個，卻讀到do i=1,12，或是1D矩陣寫成3D矩陣，或是已開啟資料寫兩次read()
list in: end of file
apparent state: unit 31 named fort.31
last format: list io
lately reading direct formatted external IO
Aborted
!txt_name(31)='xxx'之後，後面出現read(31,*) txt_name(31)
/tmp/ccI6mztx.o(.text+0x139d): In function `MAIN__':
: undefined reference to `filt25_'
副程式最後要加 include ‘filt25.f’
do_ud: off end of record !recl那nx和ny寫錯
apparent state: unit 85 named grads_5.dat
stnmap -i xxx.ctl
Name of binary data set:
/work4/tony1321/rain/200807/0710_daily.dat
Number of times in the data set: 24
Number of surface variables: 7
Number of level dependent variables: 0
Starting scan of station data binary file. Binary data file open:
/work4/tony1321/rain/200807/0710_daily.dat
Processing time = 1
Invalid station hdr found in station binary file
Possible causes: Invalid level count in hdr
Descriptor file mismatch
File not station data
Invalid relative time
levs = -1027080205 flag = -998639206 time = 1.4013e-45
Processing time = 24 !ok
Time = 24 has stn count = 415
Max reports per time: 415 reports at t = 1 Max data elements in largest report: 5
Version 2 Station map file created: rain.map
stnmap: WARNING!! This stnmap file can only be accessed
by GrADS Version 1.9b4
stnmap: WARNING!! However, GrADS Version 1.9b4 can read
both versions
stnmap: COMMENT -- use the -1 command line option to
create a map for older GrADS versions
grads
shaded會把contour蓋過去，因此要先寫shaded，再寫contour
q file 1.2……:不需開ctl即可看到各種變數
若ctl觀測為一段時間，沒再gs定義set t XX，他會自動跑下去
Constant field. Value = 2.60788e-09
'set font 0-5 調整字體
'set cstyle 0-7 調整線的類型
'set t 'tt 要空格(空格 = 等於)
'set digsize 數字'設定風標長度
'set mpdraw off' 將grads內建地圖關掉
Grads裡都用小寫較不會有問題('d Tem' ok但'd 'Tem不行)
會出現畫出來的圖亂讀，譬如讀溫度卻跑出壓力，那可能式氣壓層數寫錯，如有25層卻只寫到21層
Define error: Define not yet valid for station data
Default file is a station data file
! station data必放在第一個
Low Level I/O Error: Read error on data file
Data file name = grads.dat
Error reading 81 bytes at location 9882 !此處排序錯誤
Data Request Error: Error for variable 'inttem'
Error ocurred at column 1 Low Level I/O Error: Read error on data file
Data file name = 20091012.dat !只有一層但ctl卻寫21層
Error reading 141 bytes at location 417501
Data Request Error: Error for variable 'dbz'
Error ocurred at column 1
DISPLAY error: Invalid expression
Expression = dbz
Cannot plot color bar: No shading information
DISPLAY error: Invalid expression Expression = inttem
Data Request Warning: Request beyond file limits
Cannot plot data - all undefined values !set dfile順序要放對
Open Error: Can't open description file !開到不存在的ctl檔SET Error: No files open yet
DISPLAY error: no file open yet
Syntax Error: Invalid Operand !ctl檔名稱寫錯，沒讀到
'ccslp' not a variable or function name
Error ocurred at column 1
DISPLAY error: Invalid expression
Expression = ccslp
Low Level I/O Error: Read error on data file
Data file name = 20091012_RCWF.dat
Error reading 141 bytes at location 79524
Data Request Error: Error for variable 'dbz'
Error ocurred at column 1
DISPLAY error: Invalid expression
Expression = dbz
