Page 1 Source Listing XVVEST_ALLNORMS 2013-05-14 19:54 xvvest_allnorms.f 1 program xvvest_allnorms 2 3 allocatable datav(:,:,:,:) 4 allocatable sistr(:) 5 allocatable sigtop(:) 6 allocatable dswkv(:,:) 7 allocatable dwv(:,:) 8 double precision,allocatable::xvvest(:,:) 9 10 real uvsq,uv850sq,prcpsq,mslpsq,z500sq 11 real uvsd,uv850sd,prcpsd,mslpsd,z500sd 12 13 read(5,*)idim,jdim,mem,ne1,jstr,casecode 14 15 allocate(datav(idim,jdim,mem,12)) 16 allocate(sistr(jstr)) 17 allocate(sigtop(jstr)) 18 allocate(dswkv(idim,jdim)) 19 allocate(dwv(idim,jstr)) 20 allocate(xvvest(ne1,mem)) 21 22 open(99,form='unformatted') 23 read(99) datav 24 c test 25 do 520 j=15,17 26 do 520 i=100,102 27 do 520 k=1,mem 28 do 520 kkkk=1,12 29 write(25,*) datav(i,j,k,kkkk) 30 520 continue 31 32 close(99) 33 34 c *** Calculate standard deviation of variables *** 35 uvsq = 0.0 36 uv850sq = 0.0 37 prcpsq = 0.0 38 mslpsq = 0.0 39 z500sq = 0.0 40 41 do 15 j=1,jdim 42 do 15 i=1,idim 43 do 15 nm=1,mem 44 uvsq = uvsq + datav(i,j,nm,1)**2 + datav(i,j,nm,2)**2 + 45 & datav(i,j,nm,3)**2 + datav(i,j,nm,4)**2 + 46 & datav(i,j,nm,5)**2 + datav(i,j,nm,6)**2 47 uv850sq = uv850sq + datav(i,j,nm,1)**2 + datav(i,j,nm,2)**2 48 prcpsq = prcpsq + datav(i,j,nm,10)**2 49 mslpsq = mslpsq + datav(i,j,nm,11)**2 50 z500sq = z500sq + datav(i,j,nm,12)**2 51 15 continue 52 53 uvsd = sqrt(uvsq/float(6*idim*jdim*mem)) 54 uv850sd = sqrt(uv850sq/float(2*idim*jdim*mem)) 55 prcpsd = sqrt(prcpsq/float(idim*jdim*mem)) 56 mslpsd = sqrt(mslpsq/float(idim*jdim*mem)) 57 z500sd = sqrt(z500sq/float(idim*jdim*mem)) Page 2 Source Listing XVVEST_ALLNORMS 2013-05-14 19:54 xvvest_allnorms.f 58 59 print*, 'u,v standard deviation = ',uvsd 60 print*, 'u,v 850mb standard deviation = ',uv850sd 61 print*, 'prcp standard deviation = ',prcpsd 62 print*, 'mslp standard deviation = ',mslpsd 63 print*, 'z500 standard deviation = ',z500sd 64 65 knum0=1234 66 open(knum0,form='formatted') 67 write(knum0,*) uvsd 68 write(knum0,*) uv850sd 69 write(knum0,*) prcpsd 70 write(knum0,*) mslpsd 71 write(knum0,*) z500sd 72 close(knum0) 73 74 c *** Interpolate ensemble perturbations onto stretched grid *** 75 do 35 nm=1,mem 76 do 28 iv=1,12 77 do 27 j=1,jdim 78 do 27 i=1,idim 79 dswkv(i,j)=datav(i,j,nm,iv) 80 27 continue 81 call interp(dswkv,dwv,sistr,sigtop,jstr,idim,jdim) 82 do 30 i=1,idim 83 do 30 j=1,jstr 84 datav(i,j,nm,iv)=dwv(i,j) 85 30 continue 86 28 continue 87 35 continue 88 89 do 1300 nm=1,mem 90 ie1=0 91 do 1000 j=1,jstr 92 do 1000 i=1,idim 93 do 1100 iv=1,12 94 ie1=ie1+1 95 xvvest(ie1,nm)=datav(i,j,nm,iv) 96 1100 continue 97 1000 continue 98 1300 continue 99 100 knum1=8600+casecode 101 write(knum1) xvvest 102 103 STOP 104 END Page 3 Source Listing XVVEST_ALLNORMS 2013-05-14 19:54 Entry Points xvvest_allnorms.f ENTRY POINTS Name MAIN__ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 1000 Label 97 91,92 1100 Label 96 93 1300 Label 98 89 15 Label 51 41,42,43 27 Label 80 77,78 28 Label 86 76 30 Label 85 82,83 35 Label 87 75 520 Label 30 25,26,27,28 CASECODE Local 13 R(4) 4 scalar 13,100 DATAV Local 3 R(4) 4 4 1 ALC 15,23,29,44,45,46,47,48,49,50,79,8 4,95 DSWKV Local 6 R(4) 4 2 1 ALC 18,79,81 DWV Local 7 R(4) 4 2 1 ALC 19,81,84 FLOAT Func 53 scalar 53,54,55,56,57 I Local 26 I(4) 4 scalar 26,29,42,44,45,46,47,48,49,50,78,7 9,82,84,92,95 IDIM Local 13 I(4) 4 scalar 13,15,18,19,42,53,54,55,56,57,78,8 1,82,92 IE1 Local 90 I(4) 4 scalar 90,94,95 INTERP Subr 81 81 IV Local 76 I(4) 4 scalar 76,79,84,93,95 J Local 25 I(4) 4 scalar 25,29,41,44,45,46,47,48,49,50,77,7 9,83,84,91,95 JDIM Local 13 I(4) 4 scalar 13,15,18,41,53,54,55,56,57,77,81 JSTR Local 13 I(4) 4 scalar 13,16,17,19,81,83,91 K Local 27 I(4) 4 scalar 27,29 KKKK Local 28 I(4) 4 scalar 28,29 KNUM0 Local 65 I(4) 4 scalar 65,66,67,68,69,70,71,72 KNUM1 Local 100 I(4) 4 scalar 100,101 MEM Local 13 I(4) 4 scalar 13,15,20,27,43,53,54,55,56,57,75,8 9 MSLPSD Local 11 R(4) 4 scalar 56,62,70 MSLPSQ Local 10 R(4) 4 scalar 38,49,56 NE1 Local 13 I(4) 4 scalar 13,20 NM Local 43 I(4) 4 scalar 43,44,45,46,47,48,49,50,75,79,84,8 9,95 PRCPSD Local 11 R(4) 4 scalar 55,61,69 PRCPSQ Local 10 R(4) 4 scalar 37,48,55 SIGTOP Local 5 R(4) 4 1 1 ALC 17,81 SISTR Local 4 R(4) 4 1 1 ALC 16,81 SQRT Func 53 scalar 53,54,55,56,57 UV850SD Local 11 R(4) 4 scalar 54,60,68 UV850SQ Local 10 R(4) 4 scalar 36,47,54 Page 4 Source Listing XVVEST_ALLNORMS 2013-05-14 19:54 Symbol Table xvvest_allnorms.f Name Object Declared Type Bytes Dimen Elements Attributes References UVSD Local 11 R(4) 4 scalar 53,59,67 UVSQ Local 10 R(4) 4 scalar 35,44,53 XVVEST Local 8 R(8) 8 2 1 ALC 20,95,101 XVVEST_ALLNORMS Prog 1 Z500SD Local 11 R(4) 4 scalar 57,63,71 Z500SQ Local 10 R(4) 4 scalar 39,50,57 Page 5 Source Listing XVVEST_ALLNORMS 2013-05-14 19:54 xvvest_allnorms.f 105 Page 6 Source Listing XVVEST_ALLNORMS 2013-05-14 19:54 Subprograms/Common Blocks xvvest_allnorms.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References XVVEST_ALLNORMS Prog 1 COMPILER OPTIONS BEING USED -align nocommons -align nodcommons -align noqcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume byterecl -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_unit_star -assume old_ldout_format -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume nostd_mod_proc_name -assume norealloc_lhs -assume underscore -assume no2underscores no -auto -auto_scalar no -bintext -ccdefault default -check noargs -check noarg_temp_created -check nobounds -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check power -check noshape -check nounderflow -check nouninitialized -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest Page 7 Source Listing XVVEST_ALLNORMS 2013-05-14 19:54 xvvest_allnorms.f no -ftz -fp_model noprecise -fp_model fast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -heap_arrays 0 no -threadprivate_compat -g0 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude -O2 no -pad_source -real_size 32 no -recursive -reentrancy none no -sharable_localsaves -vec=simd -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w noargument_checking -w nodeclarations -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage -includepath : /gpfs/gp2/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/,.f,./.f,/usrx/local/intel/composerxe/mkl/include/.f, /usrx/local/intel/composerxe/tbb/include/.f,/gpfs/gp2/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/.f, /gpfs/gp2/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/.f,/usr/local/include/.f,/usr/lib/gcc/x86_64-redhat-linux/4.4.6/include/.f, /usr/include/.f,/usr/include/.f -list filename : xvvest_allnorms.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100