source: palm/trunk/UTIL/chemistry/gasphase_preproc/kpp/util/Mex_Fun.f90 @ 4446

Last change on this file since 4446 was 2696, checked in by kanani, 7 years ago

Merge of branch palm4u into trunk

File size: 1.4 KB
Line 
1
2!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
3 SUBROUTINE mexFunction( nlhs, plhs, nrhs, prhs )
4!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
5!                 Matlab Gateway for the Derivative Function Fun
6!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
7
8 USE KPP_ROOT_Model
9
10      INTEGER nlhs, nrhs
11      INTEGER plhs(*), prhs(*)
12      INTEGER mxGetPr, mxCreateFull, mxGetM, mxgetN
13      INTEGER VPtr, FPtr, RPtr, VdotPtr
14      KPP_REAL V(KPP_NVAR), F(KPP_NFIX), RCT(KPP_NREACT)
15      KPP_REAL Vdot(KPP_NVAR)
16
17! Check for the right number of input arguments
18      IF ( nrhs .ne. 3 ) THEN
19         CALL mexErrMsgTxt('Fun requires 3 input vectors: &
20     &V(KPP_NVAR), F(KPP_NFIX), RCT(KPP_NREACT)')
21      END IF 
22! Check for the right number of output arguments
23      IF ( nlhs .ne. 1 ) THEN
24         CALL mexErrMsgTxt('Fun requires 1 output vector: &
25     &Vdot(KPP_NVAR)')
26      END IF
27
28      plhs(1) = mxCreateDoubleMatrix(KPP_NVAR,1,0)
29
30      VPtr = mxGetPr(prhs(1))
31      CALL mxCopyPtrToReal8(VPtr,V,KPP_NVAR)
32     
33      FPtr = mxGetPr(prhs(2))
34      CALL mxCopyPtrToReal8(FPtr,F,KPP_NFIX)
35     
36      RPtr = mxGetPr(prhs(3))
37      CALL mxCopyPtrToReal8(RPtr,RCT,KPP_NREACT)
38
39      VdotPtr = mxGetPr(plhs(1))
40
41      CALL Fun( V, F, RCT, Vdot )
42
43      CALL mxCopyReal8ToPtr(Vdot, VdotPtr, KPP_NVAR)
44
45 END SUBROUTINE mexFunction
Note: See TracBrowser for help on using the repository browser.