[준비사항 및 참고사항]
Kato상의 일일예보 홈페이지 (ID:model, PW:result)
http://www.rain.hyarc.nagoya-u.ac.jp/~kato/resultJMA의 MSM 자료 생성 주기
00Z |-|-|-|—————–| 15시간 예측, 1시간간격 저장
03Z |——————————————-| 33시간 예측, 화일 2개로저장
06Z |——————–| 15시간 예측
09Z |—————————————| 33시간 예측참고 : map proj. = Rambert
[소스코드 풀기]
%cd ~dhkim/models/CReSS
%tar xvzf src/…..
%cd CReSS_20111208
%setenv CReSSHOME
pwd#CReSS의 홈디렉토리 지정을 하기위해 사용[모형 컴파일]
%cd $CReSSHOME/CReSS3.4.1_20111110/CReSS3.4.1m_20111110
%vi compile.conf
FFLAGS = -Mfree -byteswapio -O3 -fastsse
FC = ftn%./auto_compile.csh # log_compile에 에러 확인
# Src 디렉토리에 다음 실행화일들이 생성되어야 함
# asldata.exe, check.exe, gridata.exe, radata.exe,
# rstruct.exe, solver.exe, surface.exe, terrain.exe, unite.exe%mkdir $CReSSHOME/Work/YYYYMMDD
%cp -l *.exe $CReSSHOME/Work/YYYYMMDD/.
%cp run.csh $CReSSHOME/Work/YYYYMMDD/.
[IC & BC 만들기]
SST 자료 변환
자료취득 : NEAR-GOOS에서 MGDSST(해상도:0.25도) 자료를 받을 수 있음
http://goos.kishou.go.jp 에 온라인 등록하고
“JMA Products”->“MGDSST"에서 받을 수 있음%cd $CReSSHOME/mgdsst2cdf
%vi configure.make #시스템에 맞추어 compiler와 옵션 등을 수정
#module을 사용할 경우 netcdf 관련 라이브러리를 loading
%make
%vi sstlist.dat #SST 자료가 있는 디렉토리 지정
%./mgdsst2cdf
%cp out/mgdsst_YYYYMMDD.dat $CReSSHOME/Work/YYYYMMDD/data.sstYYYYMMDDhhmm.bin
# out 디렉토리에 *.ctl 이 생성되므로 GrADS 등으로 그림 확인 가능IC & BC 자료 변환
자료취득 : JMA는 GSM(global)과 MSM(mesoscale) 자료를
http://database.rish.kyoto-u.ac.jp/arch/jmadata/data/gpv/original/
에서 제공하고 있다. 자료의 상세 사항은
http://database.rish.kyoto-u.ac.jp/arch/jmadata/data/gpv/original/README
를 읽어보면 되며, 화일명에서 Rjp는 일본 지역, Rgl은 전지구를 의미한다.%cd $CReSSHOME/GRIB2
%make config # 현재 디렉토리에 configure.make 가 없으면 원본을 복사함.
%vi configure.make #시스템에 맞추어 compiler와 옵션 등을 수정
# CFLAGS = -DDEC 는 64bit 시스템에서 사용%vi gpv2dat.namelist
# config_file = ‘msm.conf’ 또는 ‘gsm.conf’, ‘rsm.conf’ 로 수정
# 여기에서는 ‘msm.conf’로 사용#msm.conf에 지정되어 있는 grib2_surface.dat(2D)와 grib2_upper.dat(3D)의 내용을 알맞게 수정
#msm.conf의 “InitialOnly = .true.“의 경우는 여러가지 초기장을 읽어서 하나의 초기장으로 만듦.
%ls ../Data/Lsurf > grib2_surface.dat
# 지정날짜의 “…Rjp_Lsurf…” 위치와 파일명 지정
%ls ../Data/L-pall > grib2_upper.dat
# 지정날짜의 “…Rjp_L-pall…” 위치와 파일명 지정%make
%./gpv2dat #GRIB2에서 MSM변환
%mv data.gpv*.bin ../Work/YYYYMMDD/.
#자료 확인은 template/MSM_GPV.ctl 을 이용하여 GrADS로 확인 가능
[모형 모의]
%cd $CReSSHOME/Work/YYYYMMDD
%cp $CReSSHOME/User_Conf_HIGE/user_conf_msm33-2km .
%vi user_conf_msm33-2km
&dirname
crsdir = ‘./’
datdir = ‘./’
&dimset
numpe = 36
xsub = 6
ysub = 6
&flength
sfcast = ‘YYYY/MM/DD hh:mm’
stime = 0.e0
etime = 108000.e0 # 30시간
&gpvpram
gpvitv = 10800.e0 # 3시간 간격 저장
&sfcphys
dtgrd = 7.5e0
&outfomat
dmpitv = 3600.e0
monitv = 3600.e0%qsub -q @haeon run.csh # 또는 다음의 순서대로 실행
%./check.exe < user.conf_msm33-2km # 마지막에 normally 나오면 o.k.
%ln -s $CReSSHOME/Data/data.*.bin .
%./terrain.exe < user.conf_msm33-2km >& log_msm33-2km
# 지형자료를 프로세서 갯수에 맞게 나눔
# 마지막 메세지에 normally가 나와야 함. 에러의 경우는 abnormally 표시됨.%./surface.exe < user.conf_msm33-2km »& log_msm33-2km
# 표층 경계치(SST)를 프로세서 갯수에 맞게 나눔
# 마지막에 normally 표시 확인%./gridata.exe < user.conf_msm33-2km »& log_msm33-2km
# 3차원장 경계치를 프로세서 갯수에 맞게 나눔
# 마지막에 normally 표시 확인%qsub -q @haeon run.csh ### 모형 모의 시작 ###
실제명령 : aprun -n 36 ./solver.exe < user.conf_msm33-2km >& log.solver%./unite.exe < user.conf_msm33-2km »& log_msm33-2km
끝.
