source: trunk/fsource/mcsasubs/DOLLASE.FOR @ 971

Last change on this file since 971 was 971, checked in by vondreele, 9 years ago

update binwin2.7 with new G77 pyd files.
add mcsasubs to fsource directory
some mods to GSASIImath.py to accommodate them

File size: 723 bytes
Line 
1      SUBROUTINE DOLLASE(CA,RAT,CORR)
2
3!PURPOSE: Calculate MArch-Dollase function
4
5      INCLUDE       '../INCLDS/COPYRIGT.FOR' 
6
7!PSEUDOCODE:
8
9!CALLING ARGUMENTS:
10
11      REAL*4        CA                  !Cos(alpha)**2
12      REAL*4        RAT                 !PO ellipsoid axis ratio
13      REAL*4        CORR                !March-Dollase function
14
15!INCLUDE STATEMENTS:
16
17!LOCAL VARIABLES:
18
19      REAL*4        SA                  !Sin(alpha)**2
20      REAL*4        A                   !Intermediate value
21
22!SUBROUTINES CALLED:
23
24!FUNCTION DEFINITIONS:
25
26!DATA STATEMENTS:
27
28!CODE:
29
30      SA = 1.0-CA
31      A = CA*RAT**2+SA/RAT
32      A = 1.0/SQRT(A)
33      CORR = A**3
34      RETURN
35      END
Note: See TracBrowser for help on using the repository browser.