global.c Source File

00001 /*---------------------------------------------------------------------------*/
00002 /* SCCS Information: %W%    %G% */
00003 /*---------------------------------------------------------------------------*/
00004 #include <stdlib.h>
00005 #include "utilf.h"
00006 #include "markers.h"
00007 #include "surfaces.h"
00008 #include "extra.h"
00009 #include "global.h"
00010 #include "mymalloc.h"
00011 
00012 /*------------------------------------------------------------------*/
00013 /*                 Computational arrays                             */
00014 /*------------------------------------------------------------------*/
00015 real2D p;
00016 real2D u, v;
00017 real2D S11, S22, S12;
00018 real2D cc;
00019 real2D sxp, syp, sxu, syu, sxv, syv, ap, au, av;
00020 real2D rz1dr, rz2dr, rr1dz, rr2dz;
00021 real2D a1, a2, a3, divs;
00022 real2D work, work1;
00023 
00024 /*------------------------------------------------------------------*/
00025 /*                 Physical variables                               */
00026 /*------------------------------------------------------------------*/
00027 PhysStruct phys;
00028  
00029 /*------------------------------------------------------------------*/
00030 /*                 Numerical variables                              */
00031 /*------------------------------------------------------------------*/
00032 NumStruct num;
00033 
00034 /*------------------------------------------------------------------*/
00035 /*                 Output variables                                 */
00036 /*------------------------------------------------------------------*/
00037 OutStruct out;
00038 
00039 /*------------------------------------------------------------------*/
00040 /*                 Markers variables                                */
00041 /*------------------------------------------------------------------*/
00042 interface *interfaces;
00043 int Nint;
00044 interface_cut cutp, cutu, cutv;
00045 
00046 
00047 void initglobal()
00048 {
00049   p = (real2D) mymalloc(sizeof(real), num.nx, num.ny);
00050   u = (real2D) mymalloc(sizeof(real), num.nx + 1, num.ny); 
00051   v = (real2D) mymalloc(sizeof(real), num.nx, num.ny + 1); 
00052   S11 = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00053   S22 = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00054   S12 = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00055 
00056   cc = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00057   sxp = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00058   syp = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00059   sxu = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00060   syu = (real2D) mymalloc(sizeof(real), num.nx + 1, num.ny); 
00061   sxv = (real2D) mymalloc(sizeof(real), num.nx, num.ny + 1);
00062   syv = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00063   ap = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00064   au = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00065   av = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00066 
00067   rz1dr = (real2D) mymalloc(sizeof(real), num.nx, num.ny);
00068   rz2dr = (real2D) mymalloc(sizeof(real), num.nx, num.ny);
00069   rr1dz = (real2D) mymalloc(sizeof(real), num.nx, num.ny);
00070   rr2dz = (real2D) mymalloc(sizeof(real), num.nx, num.ny);
00071   
00072   a1 = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00073   a2 = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00074   a3 = (real2D) mymalloc(sizeof(real), num.nx, num.ny); 
00075   divs = (real2D) mymalloc(sizeof(real), num.nx, num.ny);
00076   work = (real2D) mymalloc(sizeof(real), num.nx, num.ny);
00077   work1 = (real2D) mymalloc(sizeof(real), num.nx, num.ny);
00078 }
00079 
00080 
00081 
00082 
00083 
00084