Procesing.pptVIP

* * * Accessing LAPACK (2 of 2) external DGELS ./dgels_interface.so” ; NAG example: /lapack-ex/node45.html ; These are transposed from the fortran example A = (/ (/ -0.57, -1.93, 2.30, -1.93, 0.15, -0.02 /), \ ; (4,6) (/ -1.28, 1.08, 0.24, 0.64, 0.30, 1.03 /), \ (/ -0.39, -0.31, 0.40, -0.66, 0.15, -1.43 /), \ (/ 0.25, -2.14,-0.35, 0.08,-2.13, 0.50 /) /)*1d0 ; must be double dimA = dimsizes(A) N = dimA(0) ; 4 M = dimA(1) ; 6 B = (/-2.67 ,-0.55 ,3.34, -0.77, 0.48, 4.10/)*1d0 ; must be double ; LAPACK wants 2D nrhs = 1 B2 = conform_dims ( (/nrhs,M/), B, 1 ) ; (1,6) B2(0,:) = B lwork = 500 ; allocate space work = new ( lwork, double, No_FillValue) ; must be double DGELS::dgelsiw(M, N, nrhs, A , B2, lwork, work ) print(B2(0,0:N-1)) 楷鹃盂笑赖瘤汇盐姓吊到澎知想崔朔腮线杯闷球柞共蘑创划度丢圃嘘癣甄ProcessingProcessing Combining NCL and Fortran in C-shell #!/usr/bin/csh # =========== NCL ============ cat ! main.ncl END_NCL load ”$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl load ”$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl“ load $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl“ external SUB ./sub.so begin ... end END_NCL # ===========FORTRAN ======== cat ! sub.f END_SUBF C NCLFORTSTART ... C NCLEND END_SUBF # =========== WRAPIT========== WRAPIT sub.f # =========== EXECUTE ======== ncl main.ncl ! main.out 挪祭员拎虫己假崎氰树坑巢吾诚尖浦崭捆影躇乌踩歹宵驻法劲君传贫躇灸ProcessingProcessing Global Variables and Scope [1 of 2] load “dummy_1.ncl” ; not aware of constants below GRAVITY = 9.8 RGAS = 204 load “dummy_2.ncl” ; can use GRAVITY and RGAS REARTH = load “dummy_3.ncl” ; can use GRAVITY, RGAS, REARTH begin ; can use GRAVITY, RGAS, REARTH : end Global Variable(s) by definition: can