! Last change: SS 19 Sep 100 1:31 am PROGRAM Setup_RCM_winds USE Rcm_mod_subs IMPLICIT NONE INTEGER (iprec) :: iflag, i, j CHARACTER (LEN=1) :: answer_1char CHARACTER (LEN=80) :: form_string LOGICAL :: logic_flag ! INQUIRE (FILE = 'rcmwind', EXIST = logic_flag) IF (logic_flag) THEN WRITE (*,'(T2,A)') 'WINDS INPUT FILE ALREADY EXIST' WRITE (*,'(T2,A)') 'IF YOU PROCEED, IT WILL BE OVERWRITTEN' WRITE (*,'(T2,A)',ADVANCE='NO') 'PROCEED ANYWAY ? (Y/N):______' READ (*,*) answer_1char IF (answer_1char /= 'y' .AND. answer_1char /= 'Y') THEN STOP END IF END IF ! ! pwe = zero pwn = zero hwe = zero hwn = zero sw = zero ! ! OPEN (UNIT = LUN, FILE = 'rcmwind', STATUS = 'REPLACE', FORM = 'FORMATTED') ! form_string = '(4(TR2,ES23.15))' WRITE (LUN, '(I10.10)') SIZE(pwe) WRITE (LUN,'(A80)') form_string DO j = 1, jsize DO i = 1, isize WRITE (LUN,form_string) pwe(i,j), hwn(i,j), hwe(i,j), pwn(i,j) END DO END DO ! ! form_string = '(1(TR2,ES23.15))' WRITE (LUN, '(I10.10)') SIZE(sw) WRITE (LUN,'(A80)') form_string DO j = 1, jsize WRITE (LUN,form_string) sw(j) END DO CLOSE (LUN) STOP END PROGRAM Setup_RCM_winds